[nospecial]

我是一个MATLAB的初学者,昨天老师留了个实验,用matlab求一个系统传递函数的未知参数.下面是我想的框架.主要的不是求解的原理问题,就是程序的可行性,有些注释,希望给些指点.

syms s
y=ilaplace(0.35/((s+0.5)*(s+0.7)));
for t=0:3
g(t)=y;
end /*t=0到3 时,将Y的值赋予g(t)*/

A=[g(1),g(2);g(2),g(3)]; /*方程组-g(0)=a1*g(1)+a2*g(2) ;-g(1)=a1*g(2)+a2*g(3)求a1.a2 */
B=[-g(0),-g(1)];
C=[a1,a2];
C=B/A; /*求出a1,a2 */
x=solve('1+ai*x+a2*x^2') /*1+ai*x+a2*x^2=0 的根*/
s1=log(x1);
s2=log(x2);
D=[1,1;exp(s1),exp(s2)];
E=[g(0),g(1)];
F=[c1,c2];
F=E/F; /*求出c1,c2*/