QQ登录

只需要一步,快速开始

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

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

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

1

主题

2

听众

35

积分

升级  31.58%

该用户从未签到

自我介绍
20100103重要的日子
跳转到指定楼层
1#
发表于 2009-12-17 20:17 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
%运用Jacobi迭代拟合出来的u关于x和y的矩阵0 Q0 K7 v- L. G  F/ |5 W5 l
function J=jac(A,b,u0,eps)# r! K2 ?5 A& M1 g) O
if nargin==38 F& U+ Y0 _+ n( Z- v( _
    eps=1.0e-8
6 I; b# C$ T/ a. Relseif nargin<3
# h* K8 C7 J. m9 ^) k# r    'error'
& k8 V; |, n$ p+ i' o    return
' P1 e7 S% u# }/ w2 ^end0 j/ v( z; z4 @' g+ U/ w1 a
  _9 o7 b* V: {) v8 W6 N
%定义内部节点矩阵u0* O4 \$ M- l) Z
h=1;k=1;9 @  A7 m( j' s& A) X6 D' ^# {
x=0:h:17;y=0:k:10;0 f+ x' R& I$ ?8 q7 _  ~3 A6 C
e=length(x)-2;f=length(y)-2;! Y3 O/ e) j- q! D% m( ]" [) z9 @
u=zeros(e,f);
: f& K+ b4 y3 \# `; @u0=u;; k) l2 W% M/ s7 c

# a# W# d% e- S# |, |%定义外部节点p+ w: E% e" n* A
p=zeros(e+2,e+2);
" T# L$ h' M& E- O% ip(1,1:f+2)=0;p(e+2,1:f+2)=0;       , I2 v% J! n0 x0 u9 ~" @3 h
p(1:e+2,1)=100;p(1:e+2,f+2)=100;( x, T+ U* O, v) k" X$ ~% Z
0 b, U- r% X7 k& G) D: _! a+ M( Q
%定义系数矩阵A! G2 d  u2 n6 @1 E6 u" {
A=zeros(e*f,e*f);4 q: o- }% k$ f  U9 V8 V: H: a6 X% g
B=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
/ F' ?4 q: S* O% ad1=ones(e,1);d2=ones(e-1,1);; I2 r$ t. Y" q$ R9 ?
M=4*diag(d1)-diag(d2,1)-diag(d2,-1);3 h9 \' E) t& u+ Y
N=-eye(e);/ d) p% _: S0 r9 r% b
B{1,1}=M;B{e,e}=M
2 ], V0 ~* b% o+ S0 D) I$ Afor i=2:e-1
3 D. a3 t* g, g    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N9 M1 j; M5 n) y. C7 |5 M0 S
end
6 e1 G2 A$ D3 y& R3 qA=cell2mat(B);  0 c4 f, v& n- ?" E) g
这里总是显示
* \1 ~% G5 q- D3 ?1 L0 y  e0 `??? function J=jac(A,b,u0,eps)7 u! g9 Y/ Q' c: C
    |! Y8 b( z. z& T9 H
Error: Function definitions are not permitted at the prompt or in scripts.

3 N1 G. }+ D1 t6 d  i) X( U8 h/ ?1 B& d1 U4 C3 |
%定义b
( o: X" T. I. Q& Y3 q) h5 Sb=ones(e*f,1);' z/ R* b9 Q. S5 y" V% R5 d6 e/ @
for i=1:e1 s' L* T0 J, [3 h$ `
    for j=1:f
5 T$ f# w3 S5 x5 {& M        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)5 a% w  t2 c0 p- O* Q' ^
    end7 p* z& h* f, T& l
end
8 t% g3 H; l% r" E( e' G%运用Jacobi迭代法计算
1 L( C& D1 s( U+ }  ~% `D=diag(diag(A));, e" b* j5 h. n1 y% b/ d
D=inv(D);% Q9 }- Y( j( o' |& y! i* L! c. X
L=tril(A,-1);
. u8 }2 S0 S$ C+ O! LU=triu(A,1);, B: P' l& F: C$ o# S( Z* A2 e
B=-D*(L+U);
' t9 V" V' K: u5 h! K2 Of=D*b;$ X# Y* n4 J( |% x6 O/ t4 D
J=B*u0+f;0 ~( c2 C: K9 j9 H; F7 v' J" R) C7 T7 X
while norm(J-u0)>=eps
; z0 U, u! {2 G7 n3 cx0=J;
, k, |+ R' E7 a, F1 T- |! YJ=B*u0+f;
! P# n+ i; T3 }0 \( yend4 ~% F+ D2 N- ?) m$ d" g. D
return
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 10:01 , Processed in 0.424459 second(s), 73 queries .

    回顶部