Trying to plot steady state, transient, total response but I can't seem to get it to work.

29 次查看（过去 30 天）
Jeffrey
Jeffrey 2014-5-9
回答： A Jenkins 2014-5-9
Case 1
clc
clear all
%%Initial Conditions
zeta = 0.05;
r = 0.3;
naturalfreq = 3.198;
dampfreq = 3.194;
phi = 0.033;
X = 0.048;
omega = 0.96;
Xinitial = 0.048;
phizero = 0.0599;
%%Inputting time
for t = 1:2:20
transient(t) = Xinitial*exp(-zeta*naturalfreq*t)*cos((dampfreq*t)-phizero);
steadystate(t) = X*cos((omega*t)-phi);
total(t) = Xinitial*exp(-zeta*naturalfreq*t)*cos((dampfreq*t)-phizero) + X*cos((omega*t)-phi);
end
time = [1:19];
%%Plot
plot(time,transient,'k');
hold on
plot(time,steadystate,'b');
hold on
plot(time,total,'r');
legend('Transient','Steady State','Total');
xlabel('Time');
ylabel('Response');
title('Case I');
This is the code I'm using and I'm trying to graph to transient, steady state, total response with 10 time periods. It's suppose to look like a sinusoidal graph and similar to what this code produces. However, my code produces it in a triangle like shape when it should be curves.
I tried to solve it by trying to make more points from the for loop but I can't do 1:0.5:20 or something similar because it won't let me index non whole numbers. I'm completely lost on how I should approach this problem any help or suggestion would be greatly appreciated!

请先登录，再回答此问题。

回答（1 个）

A Jenkins
A Jenkins 2014-5-9
%%Initial Conditions
zeta = 0.05;
r = 0.3;
naturalfreq = 3.198;
dampfreq = 3.194;
phi = 0.033;
X = 0.048;
omega = 0.96;
Xinitial = 0.048;
phizero = 0.0599;
%%Inputting time
idx=1;
time = [1:.01:20];
for t = time
transient(idx) = Xinitial*exp(-zeta*naturalfreq*t)*cos((dampfreq*t)-phizero);
steadystate(idx) = X*cos((omega*t)-phi);
total(idx) = Xinitial*exp(-zeta*naturalfreq*t)*cos((dampfreq*t)-phizero) + X*cos((omega*t)-phi);
idx=idx+1;
end
%%Plot
plot(time,transient,'k');
hold on
plot(time,steadystate,'b');
hold on
plot(time,total,'r');
legend('Transient','Steady State','Total');
xlabel('Time');
ylabel('Response');
title('Case I');

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by