大家好，又见面了，我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。
Jetbrains全系列IDE稳定放心使用
clear;clc;
x=[2 3 4 5 7 8 10 11 14 15 16 18 19];
y=[106。42 108。26 109。58 109。5 110 109。93 110。49 110。
59 110。6 110。9 110。76 111 111。2];
x3=min(x):0。1:max(x);%可以放在循环外
%v=zeros(3,6);没有必要，需要的话可以，v=[];
for n=1:3
v(n,1:n 3)=polyfit(x,y,n 2);
y3(:,n)=polyval(v(n,1:n 3),x3);
text=strcat(num2str(n 2),’次多项式模拟函数为’);
disp(text);
disp(poly2str(v(n,1:n 3),’x’));
end
%用多项式拟合一般效果都不好，用分式拟合效果好
fun=@(a,t) (a(1)*t。
^2 a(2)*t a(3))。/ (t a(4))
a0=[0。1 100 -100 -1];
a=nlinfit(x,y,fun,a0)
yf=fun(a,x3);
plot(x,y,’o’,x3,y3,x3,yf,’*’)
结果：
3次多项式模拟函数为
0。
0032566 x^3 – 0。12242 x^2 1。5113 x 104。4824
4次多项式模拟函数为
-0。00048007 x^4 0。023554 x^3 – 0。40963 x^2 3。
0473 x 102。0444
5次多项式模拟函数为
9。6416e-005 x^5 – 0。0054588 x^4 0。11758 x^3 – 1。2012 x^2 5。9223 x
98。
5719
fun =
@(a,t)(a(1)*t。^2 a(2)*t a(3))。/(t a(4))
a =
0。0541 110。
2197 -117。8171 -1。0335。
全部
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/188920.html原文链接：https://javaforall.net
