Answered
Issues Using Dirac Function Trying to add Impulse Force at specific time
Hi Michael, The ODE solvers aren't set up to handle Dirac impulse functions this way. As I read the code (make sure to double ...

2 years ago | 0

| accepted

Answered
Errors using ode45, state representation with time variant parameter matrices
%System Parameters a = 4.5; M = 2.5e3; I = 3.2e2; k = 5e5; c = 3.6e3; v0 = 18; g = 9.81; rho = 8e2; EI = 7e8; L = 400; %Para...

2 years ago | 0

| accepted

Answered
Obscure error in mvncdf()
Walking through the debugger for this code: %{ aaa1=[4.546541072119, -19.52465262]; aaa2=[4.548657475153, -17.361228432903]; ...

2 years ago | 0

| accepted

Answered
Extracting data from SimulationOutput from multiple simulations
signalname = 'PCC_Vabc'; PCC_V_ALL = cell2mat(arrayfun(@(out) TP_rms(get(out.logsout,signalname).Values.Data),out,'UniformOutpu...

2 years ago | 0

| accepted

Answered
Finding a definite integral
Hi Sharon, It seems like we want to compute the following integral syms x u(x) k(x) w(x) I = int((u(x)+k(x))*w(x),x,0,1,'Hold...

2 years ago | 0

Answered
How can i tune Q and R parameters of a LQI controller?
This Answer thread may be of interest. Among other things, it shows how to use the Symmetric Root Locus (SRL) technique to desig...

2 years ago | 0

Answered
Discrete filter bode shows wrong results
FIR filters are typically specified with coefficients in ascending powers of z^-1, but the the default for tf in the Control Sys...

2 years ago | 0

Question

Should simplify() Work Better on Symbolic Expressions Involving heavisde() with an Ineteger Argument?
I was playing around with what I thought was a simple problem but turned out to not be. Let x(n) be the input to a discrete-tim...

2 years ago | 1 answer | 0

1

answer

Answered
Why does freqz function make a phase difference?
Hi Mohammad, See code below for corrections %calcuating fourier transform of x by freqz function. A=1; x=[1 2 3 4]; The arg...

2 years ago | 0

Answered
Integration of real variable function on arbitrary interval.
Why would you expect the result to be 3/4? Plotting the function, assuming we want the real cube root plot(-1:.01:1,nthroot(-1...

2 years ago | 0

| accepted

Answered
Create FFT filter of small FFT amplitudes and replot orginal signal
Hi Coleman, If I understand your question correctly .... Once you do this: P2 = abs(Datafft/L); P1 = P2(1:L/2+1); P1(2:end-...

2 years ago | 0

Answered
How to convert symbolic transfer function to state space?
Hi Sam, I was hoping that the OP would show a little more effort before giving the answer. Anyway ... The code should be modi...

2 years ago | 1

Answered
How to find the frequency and amplitude of an oscillating signal?
Hi Rohan, For this signal, the output of fft without windowing gives a very good estimate of the amplitude and frequency of the...

2 years ago | 1

Answered
Simulink add Transferfunction as binominal equation
Hi Gunnar, You can use the LTI System block and enter as the block parameter: p(1)*tf(1,[p(2) 1])^p(3) assuming that only the ...

2 years ago | 0

Answered
finding the coordinate where imaginary value is zero in nyquist plot
The doc page for nyquist shows how to use it to return the real (Re) and imaginary (Im) parts along with the frequency vector (w...

2 years ago | 0

Answered
Getting trouble in EKF for pendulum (not tracking perfectly)
The state prediction equation needs to be modified, see below. I didn't verify the rest of the equations. g = 9.81; ...

2 years ago | 0

Answered
string match and replace string with matching number from another matrix
From the desription, it sounds like the StimNo file contains a cell array (assuming all file names are unique) stimno = {'air.w...

2 years ago | 0

| accepted

Answered
Observer gain matrix calculations give warning when replicated with Matlab
Hi NGiannis, Here's the plant model and desired pole locations Gp=tf(1, [1 0 1 0]); s1=-1+i; %Desired closed loop poles s2=-...

2 years ago | 0

Answered
Root Locus Plot Is Way Off
Hi Arpad, Let's take a look at that root locus for the uncompensated system. I asusme that the red x's indicate the desired clo...

2 years ago | 0

Answered
lla2eci conversion for geo sat
Hi Derrick, I think the effect we're seeing here is a slight angular offset betweent the ECI z-axis and the Earth's axis of rot...

2 years ago | 0

| accepted

Answered
Unable to find transfer functon using model linearizer
Hi S L, Can you provide the parameters for your sinestream input signal? I couldn't get the Model Linearizer to run to complet...

2 years ago | 0

Answered
stem a convolve signal
Because the first point in x corresponds to t = 0, so does the first point in y. t_conv has to have as many points as y. And the...

2 years ago | 0

| accepted

Answered
Error using .* Matrix dimensions must agree.
The anonymous function fun won't work unless the input arguments x and y are compatible for elementwise-multiplication with u an...

2 years ago | 0

Answered
unexpected discontinuity in graphic
Hi Konstantinos, Your plot, the plot I got on my local machine, and the plot generated here aren't quite the same,presumbably d...

2 years ago | 0

Answered
How to simply multiply a symbolic variable into a matrix?
I'm pretty sure that if u is to be symbolic then so must be A (unless you make A into a cell array, which is probably not the wa...

2 years ago | 0

Answered
How to use multiple input-output datasets in ssest function (MIMO System Identification Toolbox).
The article also describes how to use the tt input to describe MIMO data sets (and for multiple experiments if that's needed): l...

2 years ago | 0

| accepted

Answered
Update M-script parameter value Continuously based on Workspace Data fed by a Simulink Model
The To Workspace block doesn't update the output in the workspace "continuously" while the simulation is running. Also, Simulink...

2 years ago | 0

| accepted

Answered
What approach can I take to design a filter that fits (approximately) the average power spectrum of multiple signals.
Hi G, Do you have the amplitude and phase spectra for the red curve? If so, you can try using tfest, though that spike at arou...

2 years ago | 0

Answered
How to save the result of each loop separately
One option is to use a cell array to store dissimilar objects y = [3 4]; for i=1:length(y) x{i}=zeros(y(i),2) end Thoug...

2 years ago | 0

| accepted

Answered
How to enter vector coefficients to z-transform
g = [2,2.2,2.3]; a1 = [2, 5,8]; a2 =[ 1,0.3,0.4]; HZ = tf(zeros(3),'Ts',-1); for ii = 1:3 HZ(ii,ii) = tf(g(ii)*[1 2 1],...

2 years ago | 1