MATLAB的曲线拟合具体实例要具体的来点新的呗！急！各位帮帮忙最好附带相关历史和历史人物

MATLAB的曲线拟合具体实例要具体的来点新的呗！急！各位帮帮忙最好附带相关历史和历史人物
MATLAB的曲线拟合具体实例
要具体的
来点新的呗！
急！各位帮帮忙
最好附带相关历史和历史人物

MATLAB的曲线拟合具体实例要具体的来点新的呗！急！各位帮帮忙最好附带相关历史和历史人物
哈哈哈,冬瓜把网址也给了,
来个实际例子
混凝土的抗压强度随养护时间的延长而增加,现将一批混凝土作成12个试块,记录了养护日期x（日）及抗压强度y（kg/cm2）的数据：
养护时间x 2 3 4 5 7 9 12 14 17 21 28 56
抗压强度y 35+r 42+r 47+r 53+r 59+r 65+r 68+r 73+r 76+r 82+r 86+r 99+r
建立非线性回归模型,对得到的模型和系数进行检验.
注明：此题中的+r代表加上一个[-0.5,0.5]之间的随机数
%模型为：y=a+k1*exp(m*x)+k2*exp(-m*x);
clc;clear;
x=[2 3 4 5 7 9 12 14 17 21 28 56];
r=rand(1,12)-0.5;
y1=[35 42 47 53 59 65 68 73 76 82 86 99];
y=y1+r
myfunc=inline('beta(1)+beta(2)*exp(beta(4)*x)+beta(3)*exp(-beta(4)*x)','beta','x');
beta=nlinfit(x,y,myfunc,[0.5 0.5 0.5 0.5]);
a=beta(1),k1=beta(2),k2=beta(3),m=beta(4)
%test the model
xx=min(x):max(x);
yy=a+k1*exp(m*xx)+k2*exp(-m*xx);
plot(x,y,'o',xx,yy,'r')
结果：
a = 87.5244
k1 = 0.0269
k2 = -63.4591
m = 0.1083
图形见：
http://hi.baidu.com/chemical%5Fengineering/album/item/e3f4b05126c1611b377abea8.html

x=-10:10;%x的取值范围
N=length(x)%x的向量长度
a=2+rand(1,N);%产生一个以2～3之间的随机数组
b=5+rand(1,N);%产生一个以5～6之间的随机数组
y=a+b.*x.^2 ;%产生y
plot(x,y,'o')%画图，蓝色圈圈
pp=polyfit(x,y,2)%2次拟合，pp为回归系数 <...

全部展开

x=-10:10;%x的取值范围
N=length(x)%x的向量长度
a=2+rand(1,N);%产生一个以2～3之间的随机数组
b=5+rand(1,N);%产生一个以5～6之间的随机数组
y=a+b.*x.^2 ;%产生y
plot(x,y,'o')%画图，蓝色圈圈
pp=polyfit(x,y,2)%2次拟合，pp为回归系数
xx=-10:0.1:10; %再产生新的xx
yy=polyval(pp,xx);%通过回归系数和xx得到yy
hold on% 在同一个地方画图
plot(xx,yy,'r')%用回归系数和xx作图，红色

收起

路过。。。。。。

这是电脑问题吗?

%编了个例子
x=-10:10;%x的取值范围
N=length(x)%x的向量长度
a=2+rand(1,N);%产生一个以2～3之间的随机数组
b=5+rand(1,N);%产生一个以5～6之间的随机数组
y=a+b.*x.^2 ;%产生y
plot(x,y,'o')%画图，蓝色圈圈
pp=polyfit(x,y,2)%2次...

全部展开

%编了个例子
x=-10:10;%x的取值范围
N=length(x)%x的向量长度
a=2+rand(1,N);%产生一个以2～3之间的随机数组
b=5+rand(1,N);%产生一个以5～6之间的随机数组
y=a+b.*x.^2 ;%产生y
plot(x,y,'o')%画图，蓝色圈圈
pp=polyfit(x,y,2)%2次拟合，pp为回归系数
xx=-10:0.1:10; %再产生新的xx
yy=polyval(pp,xx);%通过回归系数和xx得到yy
hold on% 在同一个地方画图
plot(xx,yy,'r')%用回归系数和xx作图，红色
图形见：
http://hi.baidu.com/chemical%5Fengineering/album/item/f00308fa1e54ba9658ee903e.html

收起

x