📄 r11nl.m
字号:
x = [ 2 3 4 5 7 8 10 11 14 15 16 18 19 ];
y = [6.42 8.20 9.58 9.50 10.00 9.93 10.49 10.59 10.60 10.90 10.76 11.00 11.20];
l1=polyfit(x,y,1)
z1=l1(1)*x+l1(2);
s1=y-z1;
x1=1./x;
y1=1./y;
y2=log(y);
l2=polyfit(x1,y1,1)
z2=1./(l2(1)./x+l2(2));
s2=y-z2;
l3=polyfit(x1,y2,1)
nn=size(x);
n=nn(2);
e1=0;
e2=0;
e3=0;
e4=0;
a1=0;
a2=0;
for i=1:n
a1=a1+y(i)*exp(l3(1)/x(i));
a2=a2+exp(2*l3(1)/x(i));
end
A=exp(l3(2))
z3=A*exp(l3(1)./x);
s3=y-z3;
B=a1/a2
z4=B*exp(l3(1)./x);
s4=y-z4;
for i=1:n
e1=e1+s1(i)^2;
e2=e2+s2(i)^2;
e3=e3+s3(i)^2;
e4=e4+s4(i)^2;
end
e1,e2,e3,e4
plot(x,y,'.',x,z1,x,z2,x,z3,x,z4)
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -