声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1358|回复: 10

[编程技巧] 菜鸟求助编程循环问题(帮忙看看程序)

[复制链接]
发表于 2009-4-22 12:32 | 显示全部楼层 |阅读模式

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

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

x
:@(过程约束方程Ax=c  (我在编程序时候c=0了)
未命名.jpg
clear
for k=1:7
x0=[0  0  0  0  0  0  0  ];
x=[10 12 2 10 10 6 4];%真实值
y=[10.5 11.5 2.5 15 9.5 6.5  4.5 ];%测量值
A=[1 -1 1 0 0 0 0
    0 1 -1 -1 0 0 0
    0 0 0  1 -1 0 0
    0 0 0 0  1 -1 -1];
for i=1:7
    Q(i,i)=0.01*y(i);% 方差协方差矩阵Q对角线元素取测量值的0.01
end
       e= zeros(7,1);
e(k,1) = 1;%e为单位列向量第i个为1,其余为0
r=A*inv(y(k)); %每一个节点的残差为r(i),c=0
V=A*Q*A';%r的协方差
f=A*e ;   % ??
T=((f'*inv(V)*r).^2)/(f'*inv(V)*f); %T(k)=r'*inv(V)*r-inv(r-b*f)*inv(V)*(r-b*f);
  Tc=0.248445; %计算查表求检验统计量的临界值Tc
    if T(k)>Tc
           b=inv(f'*inv(V)*f)*(f'*inv(V)*r);         
        y=y-b*e;     %最后一步循环回去不会
         Ax=-b*e;  
          else               
    end
  end   

我感觉第八步循环不回去
??? Error using ==> minus
Matrix dimensions must agree.
Error in ==> Untitled0421 at 29
        y=y-b(k)*e;

[ 本帖最后由 prince00840 于 2009-4-22 18:10 编辑 ]
回复
分享到:

使用道具 举报

发表于 2009-4-22 15:04 | 显示全部楼层
看下本版规则! 出错提示?
LZ的编程好像亦未完成!? 给个建议编程急不得

[ 本帖最后由 ChaChing 于 2009-4-22 15:19 编辑 ]
 楼主| 发表于 2009-4-22 15:16 | 显示全部楼层
主任您好
我感觉第二步  那个e是 单位 列向量,第i个是 1,其他为零。这一步取得好像不对

还有就是第八步,实现不了循环到第一步:@(
发表于 2009-4-22 17:18 | 显示全部楼层
i=1:7      k=1:5   看下是不是这里的问题   i和k是不是得相等呀?
 楼主| 发表于 2009-4-22 18:11 | 显示全部楼层
:@(重新写的新的了,还是不知道那个把新的y转入到第一步
发表于 2009-4-23 08:28 | 显示全部楼层
再加一个循环
 楼主| 发表于 2009-4-23 11:16 | 显示全部楼层
感觉少点东西,就不知道往哪里填写:@o
发表于 2009-4-23 11:44 | 显示全部楼层
感觉你的编程思路很混乱啊。。。
1.A矩阵赋值第三行是不是漏掉了一个-1?
2. e(k,1) = 1;%e为单位列向量第i个为1,其余为0
你的注释和你的命令对不上啊??到底是行还是列??
3.for k=1:7 应该放在A赋值之后吧?否则,每次循环开始,y都会别重新赋值为初始值,当然没办法对y进行更新了。
4. else应该放在Ax=...这句前面??

[ 本帖最后由 ydlcsu 于 2009-4-23 11:49 编辑 ]

评分

1

查看全部评分

 楼主| 发表于 2009-4-23 14:02 | 显示全部楼层

回复 8楼 ydlcsu 的帖子

1、A矩阵没有错误啊
2、我用这两个式子   e= zeros(7,1);
e(k,1) = 1;得到e为单位列向量第i个为1,其余为0的 矩阵,不知道对不对
3.你说的有道理,呵呵,谢谢
4、我也不清楚
 楼主| 发表于 2009-7-21 18:04 | 显示全部楼层
怎么写成子函数形式啊?
发表于 2009-7-21 23:05 | 显示全部楼层
抱歉! 不确定"子函数形式"指什麽? 对应的英文是什麽?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-3 05:34 , Processed in 0.076744 second(s), 25 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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