What frustrates you about MATLAB? #2
It would be nice if I could open a file in the editor as read-only, click somewhere to set the editor for read only after a file...

Most effective way to solve nonhomogeneous linear ODE problem
Hello Michal, Following up on this comment, I'd be concerned about the solution of the "reference" method if that's what's bein...

Discrete Transfer Function Truncating Signal?
If the sample time of the DTF block is 0.3, and the simulation runs for 20.1 - 20.39999999 seconds, there should be 20.1/0.3 = 6...

Simplify Transfer Function Matrix from State Space Form
Hi Diego, It looks like the code used the Symbolic Math Toolbox? If so, there's no need to go that route unless you really nee...

| accepted

nufft wrong result ?
Hello Kraka, Have you tried setting up your query points to exactly include f = 0.025? Also, it wasn't clear why the plot comma...

| accepted

How to find the eigenvectors for multiple degree of system?
Hello KG, Divide each mode shape by its first element so that the first element of each mode shape is unity. M=[3 2 1;2 2 1;1 ...

Question

Does ArrayValued Affect the Efficiency of Integral()
Suppose I want to numerically integrate a function of a scalar that returns a vector f = @(x) [1 ; sin(1000*x)]; To integrate ...

### 1

mirroring a magnitude response of a filter from -1 to 1
freqz() takes a third input that allows you to specify the frequency vector. n1= 15; %order Fs= 100; %sampling frequency fc= ...

How to reformat the display output precision of e.g. the TF() function?
The function poly2str in the CST may be of interest, I think that the numerical formatting is governed by a single line, so I g...

| accepted

partial sum of a series
Hi SSBGH, To plot the function, we need a set of x-values to plot over. As you've done, define the function and the CFS coeffi...

Fourier coefficients of a function
Hi SSBGH If you want a closed form solution, use the symbolic toolbox to define f(x) syms x f(x) = .... Once you have that y...

| accepted

Fast Fourier Transform function
The code in the question is plotting R vs index, but it sounds like the goal is to plot R vs frequency (Hz). So we need to comp...

| accepted

How do I keep f as a variable in a state space model, as f is supposed to be the input variable??
Hello Arya, Assuming that f(t) is the input to the system, it is not included in the B-matrix. Rathe, the B-matrix should be: ...

Undefined function 'hinfnorm' for input arguments of type 'tf'.
Hi Kashish, If you don't have the Robust Control Toolbox for hinfnorm(), you can probably use sigma() from the Control Systems ...

| accepted

Plotting Root Locus In matlab with K being a positive integar
Can use r = rlocus(sys,K) where sys is the transfer function with numerator of 1 and K is a vector of positive integers of int...

I am getting unrealistic values using vpasolve
I had to rearrange the second and third equations to muliply through by Vcl and sqrt(f) respectively and restrict the range of t...

Fourier series plotting and improving
Hi Faisal Running the code: syms n t T = 2; w = 2*pi/T; % Exponential Fourier series C(n) = (1/T)*(int(4*exp(-1i*w*n*t),t,...

| accepted

What determines the maximum frequency when taking the FFT?
Hi L'O.G. If X is the output of Matlab's fft() function, then the corresponding frequency vector is: N = numel(X); wn = (0:(N...

I am trying to plot an input function u(t) along a time interval [0,1] but I keep getting errors when i try to plot. Can anyone help me fix the code so it can plot?
Hello Royson, Looks like some sort of control problem involving the controllability Gramian? Maybe this is what you're looking ...

| accepted

Symbolically solving vectors simultaneously
q0 = sym('q0_',[1 3]) moment=[5*q0(1)+25 6*q0(2)+45 5*q0(3)+45] solve(moment==0,q0)

NEED HELP with fir1
Hi Justin, Executing the top part of the code ... f1 = 100; f2 = 300; f3 = 800; fs = 5000; m = (0.3*f1)/(f2/2); M = rand...

Why does h2syn double the number of states?
Hey Ali, Both of those functions produce a dynamic controller that has the same number of states as the plant. So the resulting...

| accepted

Trouble finding more than one solution with solve
Here's a brute force approach that at least seems to give a reasonable result for V0 = 0. This code runs very slowly because erf...

| accepted

Need help initializing variable x and plotting probability density function
Hi Stirling, To plot the curve mean = 1930000; standardDeviation= 64000; % note theat ^2 is changed to .^2 for elementwise o...

how can I plot a bode plot of transfer function that has one variable that changes the frequency response?
Tunable Models might be worth a look. They basically let a single model be defined with a one or more parameters that have defau...

symbolic integration needed in the present code
syms h q x f(x) f0 f1(x) f2(x) f3(x) f = f0 + f1*q + f2*q.^2 + f3*q.^3; Eq0 = diff(f,x) - diff(f0,x); Eq1 = diff(f,x) - f^2; ...

sum of dirac delta function
syms d(x) t d(x) = dirac(x); t0 = 1:3:10; % a couple of terms for example y(x) = sum(d(x-t0))

| accepted

Why is tan(a) not equal to sin(a)/cos(a)
It's kind of a pain, but .... syms b alpha eq1 = b == sin(alpha)/cos(alpha); eq1 = rewrite(eq1,'tan') eq1 = lhs(eq1) == simp...

Need help finding the stability of a closed loop system and storing the values.
Assuming continuous time: Gred = tf([7.385,476.7],[1,19.39,491.5,3835,0]); Kp = 0:1:200; Kd = 0:1:10; Preallocate the matrix...