模糊

来源：互联网发布：图片cms 网站系统编辑：程序博客网时间：2024/04/29 12:04

clear all
x1=[9.205 9.513 9.17126 13.127 9.513 9.924 8.8928 10.210 9.980 9.405];
n=10;
for i1=1:n
    x2(i1)=0;
end
for i3=1:n
    for i2=(i3+1-1):n
        x2(i3)=x2(i3)+x1(i2);
    end
    x2(i3)=x2(i3)/(n-i3+1);
end
x6=-x2(1,2:n);
x3(1)=x2(1);
for i4=2:n
    x3(i4)=x3(i4-1)+x2(i4);
end
x4(1)=0;
for i5=2:n
    x4(i5)=x2(i5)-x2(i5-1);
end
x5=x4(1,2:n);
for i6=1:n-1
    x7(i6)=-(x3(i6)+x3(i6+1))/2;
end
for i7=1:n-1
    x8(i7,1)=x6(i7);
    x8(i7,2)=x7(i7);
    x8(i7,3)=1;
end
x9=inv(x8'*x8)*x8'*x5';%求解白化方程参数
a=x9(1,1);
b=x9(2,1);
c=x9(3,1);
x10=solve('x10^2+a*x10+b=0');%特征二次方程不等实根
x10=subs(x10);
x11(1,1)=exp(x10(1,1));
x11(1,2)=exp(x10(2,1));
x11(1,3)=x2(1)-x9(3,1)/x9(2,1);
x11(2,1)=x10(1,1)*exp(x10(1,1)*(n-1));
x11(2,2)=x10(2,1)*exp(x10(2,1)*(n-1));
x11(2,3)=x2(n);%线形方程求解
x12=x11(1:2,1:2);
x13=x11(:,3);
x14=x12/x13;
for i8=1:n
    x15(i8)=x14(2,1)*exp(x10(2,1)*(i8))+x14(1,1)*exp(x10(1,1)*(i8))+x9(3,1)/x9(2,1);
end
x16(1)=x15(1);
for i9=2:n
    x16(i9)=x15(i9)-x15(i9-1);
end
for i10=1:n
    r(i10)=abs(x2(i10)-x16(i10))/x2(i10)*100;
end
for i11=1:n
    x17(i11)=x16(i11)*(n-i11+1);
end
x18(n)=x17(n);
for i12=1:n-1
    x18(i12)=x17(i12)-x17(i12+1);
end
for i10=1:n
    p(i10)=abs(x1(i10)-x18(i10))/x1(i10)*100;
end

x5
x11
x2
x16
r
x1
x18
p