日本人穿衣风格:关于下列matlab程序的算法,多谢指教

来源:百度文库 编辑:杭州交通信息网 时间:2024/05/09 12:51:33
function p=jisuan(x,y,m)
m=input('times')
x=[0 0.9 1.9 3.0 3.9 5.0];y=[0 10 30 50 80 110];
psi=fliplr(eye(m+1,m+1));p=zeros(1,m+1);
psi(2,m+1)=-sum(x)/length(x);
for k=2:m
t=polyval(psi(k,:),x); t1=polyval(psi(k-1,:),x);
a=(x.*t)*t'/(t*t');
b=(t*t')/(t1*t1');
psi(k+1,:)=conv([1-a],psi(k,2:m+1))-b*psi(k-1,:);
end
for k=0:m
t=polyval(psi(k+1,:),x);
p(k+1)=y*t'/(t*t');
end
p=p*psi
能帮画算法框图更好,急用,谢谢!