1樓：匿名使用者
擬合引數，採用多元線性迴歸，程式如下：
clear;
clc;
x = [0.064 0.068 0.
072 0.074 0.078 0.
082 0.09 0.0961 0.
1153 0.13 0.21 0.
4 0.6 0.8 1];
y = [7.5 5 2.5 0 -2.
5 -5 -7.5 -8.75 -8.
75 -7.5 -5 -2.5 -1.
75 -1.375 -1.042];
x1 = 1./x;
x2 = 1./(x.^9);
myx = [ones(1,length(y)); x1; x2]';
myy = y';
[b,bint,r,rint,stats] = regress(myy,myx)
結果如下：
b =-1.8715
-0.4460
0.0000
bint =
-4.7448 1.0018
-0.8086 -0.0834
0.0000 0.0000
r =-3.3298
2.0318
3.7514
2.5732
1.7737
0.1965
-1.5877
-2.7445
-3.1088
-2.2312
-1.0052
0.4865
0.8648
1.0540
1.2755
rint =
-5.0347 -1.6250
-2.5935 6.6572
-0.6703 8.1731
-2.1996 7.3460
-3.1138 6.6611
-4.7980 5.1909
-6.4799 3.3046
-7.4448 1.9558
-7.8155 1.5979
-7.1867 2.7242
-6.0962 4.0858
-4.4445 5.4174
-3.9320 5.6616
-3.6615 5.7694
-3.3742 5.9251
stats =
0.7936 23.0733 0.0001 5.7103
其中b即為引數。
所以方程為：
y=-1.8715/x-0.4460/x^9
繪製影象：
yy = b(1)./x+b(2)./x.^9;
plot(x,y,'ro',x,yy,'b')
圖我就不貼了，自己執行就可以看到，但可以明確告訴你用你這組資料擬合成9階效果肯定非常的差，因為：y=-a/x+b/x^9 中1/x^9將是當x偏離1時一個非常大的數，而當x接近1時確實個很小的數，差值非常的大。另外要擬合到9階對資料量的需求也是非常大。
你這15組根本達不到好的要求
有問題歡迎追問，覺得有用，請給採納
2樓：匿名使用者
x = [1:0.1:10];
plot(x, y)
怎麼用matlab畫函式y=x/(x^2+1)^1.5的影象
3樓：匿名使用者
^x=-10:0.1:10;
y=x./(x.^2+1).^1.5;
plot(x,y)
grid on
%區間自己去定，我給出了-10到10
或者快速回畫圖答
ezplot('y=x/(x^2+1)^1.5')grid on
4樓：匿名使用者
x=(-100:10:100);
y=x/(x^2+1)^1.5;
plot(x,y)
用matlab繪製方程f=y/(1+x^2+y^2),在x=[-2,2],y=[-1,1]區間的圖形
5樓：匿名使用者
用baiplot3（）函
數可以繪出其空du間曲線。
x=-2:0.1:2;y=-1:0.05:1;
z=y./(1+x.^zhi2+y.^2);
plot3(x,y,z,'ro')
grid on
xlabel('x'),ylabel('y'),zlabel('z')
title('用daoplot3繪製z=y/(1+x^2+y^2)的影象');
用mesh（）函式可以繪出其空間曲面版。權
x=[-0.2:0.1:0.2];
y=[-1:0.1:1];
[x,y] = meshgrid(x,y);
z=y./(1+x.^2+y.^2);
mesh(x,y,z) %,view([0,30])
xlabel('x'),ylabel('y'),zlabel('z')
title('用mesh繪製z=y/(1+x^2+y^2)的影象');