声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 955|回复: 2

[编程技巧] 求助MATLAB高手指点下程序

[复制链接]
发表于 2006-9-27 09:30 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
程序在附件上,谢谢指点啊~!!

orth.m

1.01 KB, 下载次数: 3

回复
分享到:

使用道具 举报

发表于 2006-9-27 10:52 | 显示全部楼层
function Q=orth(n,k,s0)
n=2;
k=3;
s0=1;
       R=[1,1;2,2];
       L=[1,1;1,0];
       G=[1,0;0,1];
       C=[1,2;3,1];
       T=[zeros(n),eye(n);eye(n),zeros(n)];
       M=[L,zeros(n);zeros(n),C];
       N=[R,zeros(n);zeros(n),G];
       A=[zeros(n),G+s0*C;R+s0*L,zeros(n)];
      
       W0=[zeros([n,1]);ones([n,1])];
       syms z;
       q0=expm(-A*z)*W0;
      
      
       fs=int(q0'*q0,z,0,1);
       Q0=q0/fs;
       Q=Q0;
        for j=1:k-1
          syms x
          q(:,j)=int(expm(-A*(z-x))*(-inv(T)*M*q0),x,0,z);
            fs=zeros([2*n,1]);
           for i=0:j-1
               q0=expm(-A*z)*W0;
             fs=fs+int(q(:,j)'*q(:,i),z,0,1)*q(:,i);
           end
               
              qj=qj-fs;
               
               fs=int(qj'*qj,0,1);  
               
              
               Qj=qj/fs;
               
               q=[q0,qj];
               
               Q=[Q,Qj];
               
               q0=qj;
               
         end
  什么也没说,一点解释性的语句也看到,遇到什么问题也不说,换位思考一下,如果别人是这样的,你会如何做?
发表于 2006-9-27 14:15 | 显示全部楼层
呵呵,没错,问问题应尽量把问题细化,才可能得到满意的答案
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-20 01:44 , Processed in 0.065673 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表