matlab求系在y''(t)+5y'(t)+6y(t)=2f'(t)+8f(t) 在 y(0_)=-3,y'(0_)=
来源:学生作业帮 编辑:神马作文网作业帮 分类:综合作业 时间:2024/11/11 05:00:53
matlab求系在y''(t)+5y'(t)+6y(t)=2f'(t)+8f(t) 在 y(0_)=-3,y'(0_)=0 输入为f(t)=e^-t时的响应
这是我写的
clear;
a=[1 5 6];
b=[2 8];
t=0:0.01:10;
[A,B,C,D]=tf2ss(b,a);
sys=ss(A,B,C,D);
f=zeros(1,length(t));
zi=[0,-3]
yx=lsim(sys,f,t,zi);%零输入
m=exp(-t);
yf=lsim(sys,m,t);%零状态
y=lsim(sys,m,t,zi);%全响应
subplot(311);plot(yx);
subplot(312);plot(yf);
subplot(313);plot(y);
得到图形如下:
可是如果用yx=dsolve('D2y+5*Dy+6*y=0','y(0)=-3','Dy(0)=0','t')
可以得到yx=6*exp(-3*t)-9*exp(-2*t) 自己计算也是,那么t=0时,应该有yx=-3 ,但是图上却是-20多.
实在想不通.
这是我写的
clear;
a=[1 5 6];
b=[2 8];
t=0:0.01:10;
[A,B,C,D]=tf2ss(b,a);
sys=ss(A,B,C,D);
f=zeros(1,length(t));
zi=[0,-3]
yx=lsim(sys,f,t,zi);%零输入
m=exp(-t);
yf=lsim(sys,m,t);%零状态
y=lsim(sys,m,t,zi);%全响应
subplot(311);plot(yx);
subplot(312);plot(yf);
subplot(313);plot(y);
得到图形如下:
可是如果用yx=dsolve('D2y+5*Dy+6*y=0','y(0)=-3','Dy(0)=0','t')
可以得到yx=6*exp(-3*t)-9*exp(-2*t) 自己计算也是,那么t=0时,应该有yx=-3 ,但是图上却是-20多.
实在想不通.
原因在于,线性系统的状态空间模型不唯一,而默认的tf2ss得到的状态空间模型状态变量并不一定是y及其导数,这样就涉及到一个很重要的问题:你的初始值是按照哪组状态变量设置的?
事实上,对于本系统而言,ss2tf得到的状态空间模型是一种能控标准型,而要想把初值直接应用于系统,应该使用能观标准型才对.
事实上,对于本系统而言,ss2tf得到的状态空间模型是一种能控标准型,而要想把初值直接应用于系统,应该使用能观标准型才对.
matlab求系在y''(t)+5y'(t)+6y(t)=2f'(t)+8f(t) 在 y(0_)=-3,y'(0_)=
已知y'''(t)+4y''(t)+6y'(t)+2y(t)=f''(t)-f(t)用MATLAB如何求频率响应H(jw
怎么用matlab画y(t)=f(t+3)
以知2x+y+t=8,试求F(x,y,t)=5(x-y)^2+4y^2+3t^2最小值
一个MATLAB题 画出Y(T)=-3T^2+5 T>=0;3T^2+5 T
信号系统中的特解求法描述某系统的微分方程为 y”(t) + 3y’(t) + 2y(t) = 2f’(t) + 6f(t
MATLAB中,t=0:.1:2*pi; y=sint(t); plot(t,y)
关于信号与系统的问题,已知y''(t)+5y'(t)+6y(t)=f(t),y(0-)=1,y'(0-)=1 f(t)=
设y=f(x,t),而t是方程F(x,y,t)=0所确定的x,y的函数(F't(x,y,t)≠0),求dy/dx..
求参数方程导数x=f'(t),y=tf'(t)-f(t)
matlab中t=linspace(0,2*pi,100);x=[t;t]';y=[sin(t);cos(t)]';pl
设f(t)是二次可微函数且f''(t)不等于0 x=f'(t),y=tf'(t)-f(t),求dy/dx,d^2y/dx