最小二乘法算法程序

chengtianmmz

呵呵，前几天在论坛上写了一个帖子给eight，感谢他一语惊醒梦中人，把我从一个提问者变成了动手者，而后又继续自己的编程之旅。经过一个多星期的学习，我又有了一些体会。现在，我的作业已经完成了大部分，回头再看当初贴在论坛上的程序觉得真丢人:@L ，居然写得那么“笨”，可惜帖子已经不能修改了，只能留在那里继续丢人了。不过，那倒也记录了我成长的历程:@)
其实最小二乘法也可以写得不像我原来的那么麻烦：

1. fid=fopen('input.txt','r+');
   
2. U=fscanf(fid,'%f',a);
   
3. fclose(fid);
   
4. fid1=fopen('output1.txt','r+');
   
5. Y=fscanf(fid1,'%f',a);
   
6. fclose(fid1);
   
7. U1=U(1:a-1);
   
8. U2=U(1:a-2);
   
9. Y1=Y(1:a-1);
   
10. Y2=Y(1:a-2);
    
11. D1=-1*[0;Y1];
    
12. D2=-1*[0;0;Y2];
    
13. D3=[0;U1];
    
14. D4=[0;0;U2];
    
15. D=[D1 D2 D3 D4];
    
16. Q=inv(D'*D)*D'*Y

复制代码

下面是递推最小二乘的算法：

1. m=input('m=')
   
2. fid=fopen('input.txt','r+');
   
3. U=fscanf(fid,'%f',a);
   
4. fclose(fid);
   
5. fid1=fopen('output1.txt','r+');
   
6. Y=fscanf(fid1,'%f',a);
   
7. fclose(fid1);
   
8. U1=U(1:m-1);
   
9. U2=U(1:m-2);
   
10. Y1=Y(1:m-1);
    
11. Y2=Y(1:m-2);
    
12. Y3=Y(1:m);
    
13. D1=-1*[0;Y1];
    
14. D2=-1*[0;0;Y2];
    
15. D3=[0;U1];
    
16. D4=[0;0;U2];
    
17. D=[D1 D2 D3 D4];
    
18. Q=inv(D'*D)*D'*Y3;
    
19. P=inv(D'*D);
    
20. %以上程序是用最小二乘法计算的初值，取前m个数据%
    
21. for i=m:a-1;
    
22.     x=[-1*Y(i);-1*Y(i-1);U(i);U(i-1)];
    
23.     y=Y(i+1);
    
24.     p=1/(1+x'*P*x);
    
25.     Z=Q+P*x*p*(y-x'*Q);
    
26.     P=P-P*x*p*x'*P;
    
27.     q=norm(Z-Q)/norm(Q);
    
28.     Q=Z;
    
29.     if q<10e-6;
    
30.         k=i-m
    
31.         q
    
32.         Q
    
33.         break
    
34.     end
    
35. end
    
36. ppp='help me!'
    

复制代码

[ 本帖最后由 eight 于 2007-6-20 15:23 编辑

w89986581

一猜就知道是好学心细的女孩子，欢迎。

chengtianmmz

谢谢:@) 还得努力

ginger

原帖由 chengtianmmz 于 2007-6-20 14:41 发表

                               
登录/注册后可看大图

呵呵，前几天在论坛上写了一个帖子给eight，感谢他一语惊醒梦中人，把我从一个提问者变成了动手者，而后又继续自己的编程之旅。经过一个多星期的学习，我又有了一些体会。现在，我的作业已经完成了大部分，回头 ...

出来乍到！以后要多向你学习-----多动手，少提问！

wyyf

不错，就是注释太少

aliu

我也学习了，，谢谢

opqwer999

谢谢分享

xiaoyi123

收藏以后留着参考 学习

xyf21thc

很好 很不错

didi-baba

这个大学时学过~~~~~~~~

aiying

回复 1 # chengtianmmz 的帖子

不错， 我初次来这里，很高兴以后一起讨论

女夭彦页

学习学习~

细雨梧桐

没事看看，多多学习

caoqihhu

这论坛牛人真多！

ChaChing

ginger 发表于 2007-11-28 20:40

                               
登录/注册后可看大图

出来乍到！以后要多向你学习-----多动手，少提问！

个人以为"多动手，多讨论", 或许更好些!
多动手, 增进自己实力
多讨论, 学习别人精髓, 也比较不会偏颇