QQ登录

只需要一步,快速开始

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

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

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

1

主题

2

听众

35

积分

升级  31.58%

该用户从未签到

自我介绍
20100103重要的日子
跳转到指定楼层
1#
发表于 2009-12-17 20:17 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
%运用Jacobi迭代拟合出来的u关于x和y的矩阵
# j9 w7 _) t4 k( D+ q, afunction J=jac(A,b,u0,eps)
6 p/ V5 b: Z4 _* v. o4 r% J7 Dif nargin==3
' B. n) O+ d( ?9 ]6 M    eps=1.0e-8( D% k& J8 j4 ~6 A( ^, T! @  u- j
elseif nargin<39 |6 g1 \! D0 L! O) L% F; G! H
    'error'
; a' q3 s. b: \    return$ R  D2 x& @) I
end% D6 B) e% w9 ^- ]2 [
! Y& R6 f. R! }8 o9 ?( s' K9 c) j
%定义内部节点矩阵u01 h1 o; v  G# U9 b, j! t
h=1;k=1;8 T4 p9 o3 E' D+ c: i
x=0:h:17;y=0:k:10;
) u9 [1 L8 G8 }# t( B9 z1 Y; we=length(x)-2;f=length(y)-2;) m/ Q% G8 D; o4 U# \' n8 e& o  \
u=zeros(e,f);
, h  C4 a; H; P* {! \: F9 m" s: Hu0=u;
7 Q# q6 l' m2 T0 f0 [; t! k
/ L2 u7 g9 ~) V, U%定义外部节点p
$ m$ p7 Z. [7 P6 K: cp=zeros(e+2,e+2);, G! \5 s! c! R- s7 |6 n
p(1,1:f+2)=0;p(e+2,1:f+2)=0;       6 X; r6 o3 v3 p5 D" L
p(1:e+2,1)=100;p(1:e+2,f+2)=100;
8 T$ {& V( B# Z8 a5 x: F% n$ D" n8 J: f6 y
%定义系数矩阵A5 c1 r, W4 z- R/ K4 Y$ i4 _- y
A=zeros(e*f,e*f);
9 r' u/ o+ B9 k3 l( p# Z" Y- IB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
. T3 e' X8 j3 R0 H( r8 \: Wd1=ones(e,1);d2=ones(e-1,1);  R" @) @7 V/ K7 U/ B+ i
M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
$ b2 d, W4 x8 yN=-eye(e);8 p: j+ s8 |: p3 b2 y4 b; t5 l
B{1,1}=M;B{e,e}=M6 D% e% Y0 G! t# w1 A
for i=2:e-1
: D8 Z& `* N# J4 ]    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
  W( b( f# R% _0 nend: F' x' B& w" A9 }
A=cell2mat(B);  . h3 X% ]# C* v: S. K6 C/ |- @
这里总是显示4 j  F3 t, J& D8 M
??? function J=jac(A,b,u0,eps)9 D4 j$ D- D8 @1 u' y+ P
    |( Y4 F9 D$ H4 p* N5 w. N0 G' C
Error: Function definitions are not permitted at the prompt or in scripts.

/ b  [% T2 y8 E- i- W, S* h' T
; J: V# K/ L  E9 D; i& K8 S& ]%定义b
/ |# a1 {! Y& M7 |b=ones(e*f,1);
9 Z' z3 B  A3 S- m) hfor i=1:e8 c' f6 A: z7 P+ w$ i. s
    for j=1:f% H6 Y, j2 d: R1 e
        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
( G$ N& c. ~3 t$ p    end
8 n& ~7 d$ N# V- x# D/ j; o4 aend! Q% M" ~( ]  k- _
%运用Jacobi迭代法计算
( T# R4 Y/ i. Z0 wD=diag(diag(A));
/ \3 {' e  E8 d: N, J4 JD=inv(D);( g7 u  Z% E3 z% r8 u+ ]# N
L=tril(A,-1);7 ?  u9 R: d$ X9 r* ?# v
U=triu(A,1);
0 x$ o8 q* F% m2 D4 T$ h5 IB=-D*(L+U);& g; `8 G6 `0 p  d# a0 I4 Y
f=D*b;2 x' I/ X9 A# ]2 M9 O& g
J=B*u0+f;
; h* M* i# R  c+ r7 Ywhile norm(J-u0)>=eps+ [( {/ B" H$ X5 \. G6 a" w+ A. ?
x0=J;
9 j$ n/ Z9 W  uJ=B*u0+f;& W; L" r4 p! k: P; e' L7 o
end
( ~3 Q4 n0 C& F/ [. ?& u0 r* vreturn
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-7-10 14:15 , Processed in 0.499428 second(s), 73 queries .

    回顶部