QQ登录

只需要一步,快速开始

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

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

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

1

主题

2

听众

35

积分

升级  31.58%

该用户从未签到

自我介绍
20100103重要的日子
跳转到指定楼层
1#
发表于 2009-12-17 20:17 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
%运用Jacobi迭代拟合出来的u关于x和y的矩阵
. S0 z5 I9 A# D- o# M# }2 G5 Lfunction J=jac(A,b,u0,eps)
$ M8 r4 y6 h& A9 c7 G/ A( {if nargin==3# p: j; z9 H* X% ?; [7 c! Z# w
    eps=1.0e-8
0 O, b  _/ p, E( nelseif nargin<37 x5 F4 r# U3 _" _
    'error'# ^! @( T! g' `- ^& e5 E- P
    return
1 }5 j8 b8 _0 O( \$ R& Hend
8 T+ G' t2 f; _( M! I8 V- ~  [
; ^& ^: T6 M* c. ~0 n%定义内部节点矩阵u0
7 q) @# m6 m$ m" d+ z6 v+ C' f; |h=1;k=1;
% @1 l  ?: R) M/ ^6 ], E1 W0 [x=0:h:17;y=0:k:10;  T. c4 l+ b0 R  M
e=length(x)-2;f=length(y)-2;
$ q* Z9 j& C1 G$ }u=zeros(e,f);, s  N3 N  O3 P
u0=u;4 ]0 G! b9 O* J2 j2 d7 G! H
* X8 E6 O) v7 M3 ~7 q" z
%定义外部节点p
! b, Q( a2 o- J! Ap=zeros(e+2,e+2);4 _' ~7 ^4 Y" B1 c% Q: A5 [/ s" f
p(1,1:f+2)=0;p(e+2,1:f+2)=0;      
" Z5 `. }* y8 d! Lp(1:e+2,1)=100;p(1:e+2,f+2)=100;
9 B' V) ~! \+ m! k+ @' x4 b1 ~0 {6 V/ a8 A
%定义系数矩阵A
$ q. z6 |2 F7 l' t$ ~A=zeros(e*f,e*f);
9 M, S/ o4 t2 \% F, P0 c5 YB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);# V" M. \; H: y4 l- [, s  V1 C. Q
d1=ones(e,1);d2=ones(e-1,1);2 R+ g  s% B, J
M=4*diag(d1)-diag(d2,1)-diag(d2,-1);8 A. l8 J6 P: @8 a7 K8 C, H
N=-eye(e);
: l& f- ^9 P0 S# v7 tB{1,1}=M;B{e,e}=M# W# [' q2 }2 [. e  Q
for i=2:e-1
- k( R% H5 r6 a9 @$ A; T. ?( k    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
4 h" H# G; c, p6 ?0 f" Wend. \# T3 I1 C# m! i4 k) k
A=cell2mat(B);  
! X; H: ~4 Y7 P( G+ _6 n这里总是显示
4 \( t: C6 J5 ]5 c! y??? function J=jac(A,b,u0,eps)
# v) ^+ ^. T3 [& c    |
! U+ m8 U7 H6 h( oError: Function definitions are not permitted at the prompt or in scripts.

' i9 b9 r0 B9 X" J: P, ^) J
  t. E, g; Q; M* M%定义b! N3 D; R8 c: T: W/ j
b=ones(e*f,1);- k6 |  X# Q( u" R1 P$ _$ F  T
for i=1:e
' H+ i; i& o5 M# l. |% T    for j=1:f
" g4 f/ u( C  x        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
$ l9 U% X. O/ m: Q1 c: a5 a    end
* N2 ?) S  g% A& W- X  U- M3 Eend
+ g1 s$ A. u7 x' g- N; M" C%运用Jacobi迭代法计算
$ D# `; [5 T" \+ ~( U* xD=diag(diag(A));
* n- l, L* P" b( jD=inv(D);3 _4 _( Q: _5 S6 ^. G
L=tril(A,-1);7 \& F6 U1 e0 ?0 |/ u% V# c
U=triu(A,1);
& Y, E2 v9 F6 r6 U! e8 KB=-D*(L+U);, d! [  \. K" a$ j& q  H+ f) a
f=D*b;
$ W7 {  b# |: ~5 d* B+ DJ=B*u0+f;$ F+ T& ]/ }. J- l/ K, ?! f
while norm(J-u0)>=eps
+ X! @4 @8 C0 `5 @x0=J;
9 k% K, e  I) c1 iJ=B*u0+f;3 o* K5 S- F/ u. B" O9 U3 B
end
1 u) `  k. J, g0 Nreturn
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-5-26 18:25 , Processed in 0.389461 second(s), 73 queries .

    回顶部