QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3077|回复: 3
打印 上一主题 下一主题

求助:关于分块矩阵的还原,为什么实行不了

[复制链接]
字体大小: 正常 放大

1

主题

2

听众

35

积分

升级  31.58%

该用户从未签到

自我介绍
20100103重要的日子
跳转到指定楼层
1#
发表于 2009-12-17 20:17 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
%运用Jacobi迭代拟合出来的u关于x和y的矩阵( w2 B" F7 e; m% w* q* `$ D
function J=jac(A,b,u0,eps)
$ Z2 h' z  Y% Z* Q6 {if nargin==3
5 r3 o6 b) Q6 c8 O; ]$ a    eps=1.0e-8. N+ w9 @% X3 P9 }, x+ q- C
elseif nargin<3  I! y5 |* X- W& B
    'error'8 y5 I/ N! X7 n( R/ T7 n& t
    return
, R, X5 n5 W1 M& pend7 D) w! G0 S* m( [9 ]1 p9 S5 R7 B

! O! e" x; j. F! y) b/ K# n5 O%定义内部节点矩阵u0
  d! Y2 c9 W- g; t2 F* x5 j3 z! _h=1;k=1;
& {3 Z1 }5 X6 {2 v: k- jx=0:h:17;y=0:k:10;
% N: ]$ _4 J2 De=length(x)-2;f=length(y)-2;4 o) H5 A% m- T2 M1 l4 ?5 |1 N. }
u=zeros(e,f);
0 [& p6 l* ]0 X) b" [' Ru0=u;
# ?1 w; {- l. T' Y# S: n$ a
* v; [8 C; v6 ~* a5 Z%定义外部节点p% Q! l$ N. ?' I7 g9 l' k; f# ^% ?
p=zeros(e+2,e+2);7 s# y7 S" t' U; ]7 `
p(1,1:f+2)=0;p(e+2,1:f+2)=0;      
  u2 ]: b- i, y2 S8 pp(1:e+2,1)=100;p(1:e+2,f+2)=100;3 u5 O6 m/ {) M  r. p+ a! u& s

" T* c! J8 w1 c3 q%定义系数矩阵A* P" N( p0 H3 Z7 T9 k
A=zeros(e*f,e*f);
# j* L6 [1 y+ h. `/ I+ QB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
$ E7 O1 |6 y/ ~4 c# q* Sd1=ones(e,1);d2=ones(e-1,1);
7 Y. b! X- s5 A, o1 s, nM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
3 ~6 x: Z+ d' O4 w! aN=-eye(e);
' t# u4 l9 }$ @: U3 EB{1,1}=M;B{e,e}=M
: x# q* ^3 e" kfor i=2:e-1
! ~$ P& P, D! q) I; b- _* t    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
' {3 R; |7 L6 m' y8 {3 Vend
4 y  \; {& I# f. LA=cell2mat(B);  0 h1 K: }5 ]: m* n4 V0 C
这里总是显示
+ `7 o8 ^1 H& M1 Y0 ]: L6 S/ ~??? function J=jac(A,b,u0,eps)
; e' j. X: g4 l; G" ]    |* a, M1 o. K; @, p9 D/ H& A* H
Error: Function definitions are not permitted at the prompt or in scripts.
" n3 X& s' y& W+ }( E/ ]5 ]( N8 Z

% R/ ]9 F: [; q4 y8 V%定义b$ t) M" I4 J; }# ~5 n9 y, s
b=ones(e*f,1);0 H$ G6 b/ Q  N8 k. _
for i=1:e
  ~' p2 I- z4 p1 [) H4 d- \    for j=1:f" v0 t+ \+ z% N
        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
3 m' t4 C/ F/ s$ m7 H& U" t4 j    end3 u! {+ ^: s: X% u9 d  a; ]/ u9 N
end( m* d) ^& P$ {$ o
%运用Jacobi迭代法计算+ K% Z0 w; O, w4 x; a" y# _0 K
D=diag(diag(A));
, t% A5 A+ c9 g1 w8 C) vD=inv(D);
" v& M1 c3 [9 N* B' ~# M. n" XL=tril(A,-1);+ p  _7 y6 Y5 J" w- d0 E( |
U=triu(A,1);
3 g% g9 ?/ F' r* @B=-D*(L+U);! }( ?# ~7 w* \
f=D*b;1 g; }3 y+ `! T. a6 @/ `
J=B*u0+f;
/ X7 u+ g" q# B( t8 I' p  }while norm(J-u0)>=eps1 ?' h1 d0 R" V2 u& r- q* q7 @
x0=J;( L/ O4 g9 e4 n/ U* P" _& W
J=B*u0+f;7 k- S$ G5 o" j4 o$ L
end
3 v& J. `% K2 j) y- creturn
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

1

主题

2

听众

35

积分

升级  31.58%

该用户从未签到

自我介绍
20100103重要的日子
自己顶一下,拜托哪位高手指点一下,纠结这个矩阵的还原,想了好多方法还是不行~~实在想不出哪里出错了
回复

使用道具 举报

BenCam 实名认证       

9

主题

6

听众

89

积分

该用户从未签到

自我介绍
200 字节以内
不支持自定义 Discuz! 代码
回复

使用道具 举报

madio        

3万

主题

1312

听众

5万

积分

  • TA的每日心情
    奋斗
    2024-7-1 22:21
  • 签到天数: 2014 天

    [LV.Master]伴坛终老

    自我介绍
    数学中国站长

    社区QQ达人 邮箱绑定达人 优秀斑竹奖 发帖功臣 风雨历程奖 新人进步奖 最具活力勋章

    群组数学建模培训课堂1

    群组数学中国美赛辅助报名

    群组Matlab讨论组

    群组2013认证赛A题讨论群组

    群组2013认证赛C题讨论群组

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-8-25 04:43 , Processed in 0.378044 second(s), 73 queries .

    回顶部