新建 文本文档 (2).txt

回归分析几个例子。都是源码。没有采用M文件。

X和Y的长度不等啊，X是16个数，Y是8个数。且X的前8个都是1。如果用X的后8个数与Y回归，则： 

clear;clc 
X=[1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0; 1.5 2.0 3.0 4.5 7.5 9.1 10.5 12.0]' 
Y=[5.6 6.6 7.2 7.8 10.1 10.8 13.5 16.5]' 

p=polyfit(X(:,2),Y,1) 
xx=linspace(1.5,12,30); 
yy=polyval(p,xx); 
plot(X(:,2),Y,'o',xx,yy) 
y=poly2sym(p,'x') 

运行结果： 
p = 
0.8950 4.1575 


y = 0.8950*x+4.1575 

clc;clear; 
x=1:10 
y=3*sin(2*x)./x 
myfunc=inline('beta(2)*sin(beta(1)*x)./x','beta','x'); 
beta=nlinfit(x,y,myfunc,[0.5 0.5]); 
a=beta(1),b=beta(2) 


clc;clear 
x=1:5; 
y=[2 4 5 7 9]; 
x2=x(1):0.05:x(end); 
pp=polyfit(x,y,1)%pp为回归系数 
y2=polyval(pp,x2);%代入回归方程 
figure%再画一个图 
plot(x,y,'*',x2,y2) 
y=poly2sym(pp,'x')%转化为代数方程形式