Paul

Last seen: Today

198 total contributions since 2011

Paul's Badges

Answered
Simulink computational load of each subsystem (or block)
Have you looked at the doc pages for the the Simulink Profiler, Solver Profiler, and various pages on how to optimize performanc...

2 hours ago | 0

Answered
How do you use rlocus with multiple feedback loops?
Start with the correct characteristic equation: 1 + D(s)*G(s)*I(s) + D(s)*G(s)*Kt = 0 (not sure why the picture shows the last...

2 hours ago | 0

Answered
How to implement an Impulse Delay in Simulink
Not sure I've ever seen that called an "impulse delay." Doesn't matter what it's called. If Td is constant, the Transport Delay...

1 day ago | 0

Answered
Integrate a bivariate normal distribution over y when x is fixed
One simple approach is to evaluate f(a,y) for vector of y from y1 to y2 and then use cumtrapz() or trapz() on the result.

1 day ago | 1

Answered
How to detrend the signal?
Have you looked at doc detrend

1 day ago | 0

Answered
How to Solve a Second-Order Differential Equation?
>> syms x y(x) eqn2 = x*diff(y,x,2) + 2*diff(y,x) - x*y == 2*exp(x); sol = dsolve(eqn2); sol = rewrite(sol,'exp'); sol = sim...

2 days ago | 1

| accepted

Answered
Creating Symbolic state space model and transforming to canonical form?
Yes, at least in principle. For example: >> A=sym('a',3) A = [ a1_1, a1_2, a1_3] [ a2_1, a2_2, a2_3] [ a3_1, a3_2, a3...

2 days ago | 0

| accepted

Answered
Designing LQR controller error
The system might be controllable in principle, but in reality seems to be problematic for using state feedback. I had trouble ev...

4 days ago | 0

| accepted

Answered
How can I plot the transfer function? HELP Please.
Because you want to make plots, I'm going to assume that you have numerical values for each parameter in your model. For exampl...

6 days ago | 1

| accepted

Answered
How should I create a Transfer function with nth order polynomial term in denominator in simulink?
The power operator works for SISO tf and zpk objects. For example: >> h=exp(-0.2*tf('s'))*zpk([],-1,1)^5 h = ...

8 days ago | 0

| accepted

Answered
why my plot is not correct?
If you're looking for a plot with odd symmetry, try: plot(x,nthroot(x,3)) % check doc nthroot for details

8 days ago | 1

Answered
Display rotation matrix with s1, c1, etc. in place of sind(theta1), cosd(theta1), etc. respectively.
Use subs: >> doc subs >> syms theta1 s1 c1 >> A=[cos(theta1) sin(theta1);-sin(theta1) cos(theta1)] A = [ cos(theta1)...

9 days ago | 0

| accepted

Answered
Why is the convolution so different from the analytical answer?
When approximating the continuous convolution with a discrete convoluton, the discrete convolution needs to be multilplied by dt...

9 days ago | 1

| accepted

Answered
Plotting responses other than step responses to transfer function
Looks like an fplot() confusion issue. This code results in the expected plot: >> syms t; y(t) = ilaplace(1/s*poly2sym(T.Numera...

9 days ago | 0

Answered
Find the 5-point DFT of the sequence:
The code posted has a few issues. The second line is: >> u = sin(j) u = 0.0000e+00 + 1.1752e+00i >> whos u Name...

11 days ago | 0

| accepted

Answered
Cannot compute Riccati solution for LQI controller
I'm not sure it's relevant, but should the fifth row of SYS.c be all zeros? I suspect the issue is with stabilizability of the ...

12 days ago | 0

| accepted

Answered
How to solve a symbolic complex equation with real and imaginary parts?
I'm not sure what "change the paramters of the solve function" means. If a general apprach is desired, maybe something along the...

13 days ago | 1

| accepted

Answered
Problem with integrating Sin to Cos
If you just want ode in terms of cos() >> ode = int(sin(n*x),x,0,pi/2) ode = (2*sin((pi*n)/4)^2)/n >> rewrite(ode,'...

16 days ago | 0

Answered
Simbolic simplification of trigonometric functions
How about this: >> simplify(f,500) ans = -(sin(x)*(- sin(x)^2 - c)^(1/2)*(sin(x)^2 - 1)^2)/(2*d^2)

16 days ago | 1

| accepted

Answered
How to generate correlated samples from a single distribution? Each sample should be correlated with the previous.
I don't think the questionn is well posed. The question and the code make it sound like there is a single distribution and the s...

20 days ago | 0

Answered
how can i use symsum to make function with a variable apart from syms k?
This code works for the specified sequences. It would need to be modified if either x[n] or y[n] take on non-zero values for n <...

20 days ago | 0

Answered
Problems with the controllability matrix.
ctrb() in the Control System Toolbox only accepts numeric inputs. If you have numeric values for the variables, you can substit...

20 days ago | 0

| accepted

Answered
Making a transformation matrix in matrix format
You're not getting 4, 1 x 4 matrices. It just looks like that because of the formatting to the diplay. >> syms th d alph l A...

20 days ago | 0

| accepted

Answered
How to shift state space models by equilibrium points using ss and lsim?
The states, inputs, and outputs of Sys are all perturbations from the equilibrium solution. The intiial condition to lsim (x0) w...

22 days ago | 0

| accepted

Answered
Differntial equations and Initial conditions
The initial conditions are set to represent the initial conditions of the system to be simulated. What other considerations migh...

22 days ago | 0

| accepted

Answered
Matlab Code for DTFT to get frequency response
n, h, and w are all defined as row vectors and the * operator is trying to do matrix multiplication. But what you really want i...

24 days ago | 0

Answered
Help with simplifying a complex symbolic expression
This code eliminates U_j0: [num,den]=numden(term1); newnum = combine(num,'log','IgnoreAnalyticConstraints',true); newterm1 = ...

24 days ago | 0

| accepted

Answered
Laplace Transform for ODE - RC Circuit
I think that last subs command is getting confused because q is defined as the solution of the ode. I'm not sure why that's a pr...

26 days ago | 1

| accepted

Answered
Fourier Transformation of dirac comb
There is one too many points in the signal y. If you plot([y y]) and zoom in you'll see the pulse at sample 1001 immediately ...

28 days ago | 1

| accepted

Answered
Convert s in j*ω in a transfer function
No need to change anything in sys. It is a tf object and either freqresp() or bode() can evaluate sys at s = j*w doc freqresp ...

28 days ago | 1

| accepted