Is there a function like "y = filter(b,a,x,zi) uses initial conditions zi for the filter delays" in fftfilt? If not, what's the most efficient way to implement this?
Filtering is a linear operation so the iniital condition response and the input response can be computed separately and then add...

2 months ago | 1

| accepted

How to put this transfer function into simulink
The transfer function in question syms s H(s) = -(1 - 0.5*exp(-10*s))/(40*s + 1)/(15*s + 1) Starting from a single input sign...

2 months ago | 0

How to get simplified, symbolic solutions for definite integrals
Declare x w to be real syms x w real fn1 = exp(-abs(x))*cos(w*x) int(fn1, x, [-inf inf])/sym(pi)

2 months ago | 0

| accepted

I am trying to use the rectangular pulse function, but I am not sure how to use it correctly to solve a Fourier problem.
Hi Connor, This question brings up some interesting aspects of Fourier analysis, and how transforms of the discrete-time sample...

2 months ago | 0

Matlab and Routh criterion for evaluation of K for stability
Hi @omar khasawneh, Without a value for K I think you'll need to use the Symbolic Math Toolbox. Your code adapts easily. I did ...

2 months ago | 1

| accepted

Lsim inputs inversion provide the same result
Interesting. A quick look at the code, at least for the case where an output argument is specified, shows that the function tha...

2 months ago | 1

| accepted

FFT of the average vs average of the FFT
For a matrix input, fft() works down the columns. Maybe this: rng(100); S = rand(5); fft(mean(S,1)) % mean for each column,...

2 months ago | 0

Continuous Convolution in Time Domain
Hello @Masoud Nateghi The code needs to use the full convolution sum to approximate the convolution integral. Also, because the...

2 months ago | 0

| accepted

Impulse and step response differences
Hello @Arcadius, The differential equation in the Question is ambiguously written. Your solution for the impulse response assum...

3 months ago | 1

| accepted

FFT of Gaussian filter in 1D
Hello @Donya Khaledyan The code in the question doesn't define t nor FGHE2, so the plots can't be recreated. Code below might b...

3 months ago | 1

Find indices of multiple strings within another string
It looks like using string variables with an inner loop is much faster than a cell array with cellfun, at least here on Answers ...

3 months ago | 0

| accepted

How to do Pade approximation with different order of numerator and denominator in Matlab
The Sybolic Math Toolbox has a function pade() where the order of the numerator and denominator can be specified separately. The...

3 months ago | 0

| accepted

Root Locus system description at given point not showing up
Is there any chance that you have some local version of rlocus() that is shadowing the Control System Toolbox? What is the outpu...

3 months ago | 1

Analog Bessel filter output problem (nan)
The num and den inputs to filtfilt() are for a discrete-time filter. But besself() returns the numerator/denominator for a cont...

3 months ago | 0

Fourier expansion of dirac delta function
Increasing n to better approximate the infinte sum yields L=10; F0=10000; v0=1; t=4; x=-10:0.1:10; n=400; f=0; for i=1:n...

3 months ago | 0

| accepted

squeeze function returning "NaN"
squeeze() can't ignore anything. If the input to squeez() is 1 x 1 x 33 matrix of NaN, the ouptut will be 33 x 1 of NaN. What ...

3 months ago | 0

Sum in two variables for response of NON-LTI discrete time system
Hi @ALESSIO MANIÃ€ 1. The loop should be over n, beause each time through the loop it should be computing y(ni), i.e., a single ...

3 months ago | 0

| accepted

How can I normalize my filter?
Hi @Amy Lg Why is fs defined and then multiplied by 1e9 in the call to bilinear? Which I guess might be o.k. but perhaps isn't...

3 months ago | 0

| accepted

Question

Is There a Way to disp a sym Double Subscript?
dips() works fine with a single subscript disp(str2sym("E_omega")) Double subscript looks odd disp(str2sym("E_omega_n")) Any...

3 months ago | 1 answer | 0

### 1

What does my Root Locus response mean? Why is it different?
Hello @Joel Okanta The root locus is much more general than as suggested in the Question. Given a system described by an open l...

3 months ago | 0

Solving a nonlinear system differential equation using symbols
When using ReturnConditions, the outputs need to include Parameters and Conditions, or use only a single output syms a b c d e ...

3 months ago | 0

Transfer Function: Continuous approximation of discrete bode
Hello @John The second paragraph of the question states that P is a discrete-time transfer function. But later on the question ...

3 months ago | 0

Most relevant variables to a Simulink model
The Simulink Design Optimization toolbox can do sensitivity analysis. I've never used it, so can't say anything more than that.

3 months ago | 0

matlab function fft, Error using fft Invalid data type. First argument must be double, single, int8, uint8, int16, uint16, int32, uint32, or logical.
At this line metroX=abs(fft(x,N1)); the variable x is still a sym object. The input to fft() must be a sequence of actual num...

3 months ago | 0

Step Response of Transfer Function is Different Than Response of State Space
Is linsys1 a discrete-time or continuous-time, i.e., what is the output of linsys.Ts If the result of that command is not zero...

3 months ago | 0

| accepted

generating random numbers from poisson distribution
Check out doc poissrnd assuming that lamda*t in the question is the same as lamda in that function.

3 months ago | 0

Single extension not everywhere?
The linked blog post does not show assignment (=) included uder "Supported Operators and Functions." It still isn't: doc bsxfun...

3 months ago | 0

| accepted

How could I do a summation of (y) so I can create a signal that goes from 0 to 30 sec repeating (y) every 3 seconds in a graph.
Maybe this is the goal? Note that I modified the definition of y to make it in line with what I assumed the question means. A ...

3 months ago | 0

| accepted