t=[];
y=[];
y=log(1-y);
%线性拟合
n=8;Lx=0;Ly=0;Lxx=0;Lxy=0;Lyy=0;
for i=1:n
Lx=Lx + t(i);
Ly=Ly + y(i);
Lxx=Lxx + t(i)*t(i);
Lxy=Lxy + t(i)*y(i);
Lyy=Lyy + y(i)*y(i);
end
a=?;
b=?;
%参数和模型求取
w1=?;
w2=?;
T=1/sqrt(w1*w2)
zeta=w1*w2*T
z=[];
p=[-w1-w2];
k=w1*w2;
不知道你这个问题是否已经解决, 如果还没有解决的话: