QQ登录

只需要一步,快速开始

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

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

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

1

主题

2

听众

35

积分

升级  31.58%

该用户从未签到

自我介绍
20100103重要的日子
跳转到指定楼层
1#
发表于 2009-12-17 20:17 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
%运用Jacobi迭代拟合出来的u关于x和y的矩阵
4 f7 Q2 Y% k# h3 U5 k, g+ @function J=jac(A,b,u0,eps)
' ]6 @6 Q4 s2 s  I! bif nargin==3
5 Q: ]2 `* @) u$ B. l8 J    eps=1.0e-8) d% V4 x6 P0 H7 [- i0 I/ ~$ V
elseif nargin<3! I  x" x2 l/ B# S$ g) U
    'error'
. N  a* a+ A' I1 s% J! I    return
* P3 m5 Y5 m' ]end
0 l( R! ~- H9 J( w- k
6 r+ U+ L* @. b  f' i%定义内部节点矩阵u0. r6 V! m& K  i  @
h=1;k=1;
- I; q8 }+ ]8 ix=0:h:17;y=0:k:10;
: J% i$ s' N7 v; y( o3 _- W, w2 Ge=length(x)-2;f=length(y)-2;% S: o! n, t& t8 f2 i3 L
u=zeros(e,f);
! y) f, {% b' O' `+ z9 K. \u0=u;
4 _+ Y; ^* M! `& f  ~& p% B7 ~2 C( ^* l
%定义外部节点p' V3 D" a  `# ^
p=zeros(e+2,e+2);
2 B# O4 _" C+ K! P& _, o* Fp(1,1:f+2)=0;p(e+2,1:f+2)=0;      
( ~* Q8 e' E2 k8 F6 G5 Q7 V: T8 Hp(1:e+2,1)=100;p(1:e+2,f+2)=100;
+ l& ?# Q3 Y3 r* K% C
) j' {+ m& U( j%定义系数矩阵A
* ~5 P$ f  h3 B0 dA=zeros(e*f,e*f);8 ~6 T) D0 \3 l7 Y, Q* r
B=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
& H" G5 ]8 S8 J# fd1=ones(e,1);d2=ones(e-1,1);
- S9 U- |9 \3 r0 s9 sM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
) B* x, l5 q% t( N9 e% _& X% |N=-eye(e);. E- F7 K9 R+ q) M  f  O/ [
B{1,1}=M;B{e,e}=M
* D' {; b/ f' X9 Y9 C' H  J+ `& y7 @for i=2:e-17 T8 S. c# w5 [+ L* R
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N1 K4 T4 u0 j% F+ O) \
end
; s: Y; j( {% m( Q" \A=cell2mat(B);  - O+ L8 U7 n& m; Z) t9 g- O
这里总是显示
) V7 Q" t" V1 g1 P2 y2 c??? function J=jac(A,b,u0,eps)- H" D+ x+ Q7 k8 }# t/ x  c6 l* e
    |: d+ b0 w1 k% U* s+ h
Error: Function definitions are not permitted at the prompt or in scripts.
4 v! r1 U" q' s: \/ J
7 V; w5 G; A) E6 c" [
%定义b
1 ?7 O0 }, v+ v2 a3 k6 N! [) }! lb=ones(e*f,1);
: z' P0 K8 Z' f# q$ Jfor i=1:e
5 K$ q4 y  w* [3 w    for j=1:f
. u# {: g; X. A! v7 x        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
/ ]6 h& s0 S3 v( J    end
7 {0 F2 D& e% R. Q" V/ aend
* G+ _, h0 c' T0 X%运用Jacobi迭代法计算
/ K2 D: c  m% oD=diag(diag(A));
! A8 k$ h( l4 T) n8 P0 ]& @/ Q# SD=inv(D);& `& J. ^) {7 G) o. s8 ?
L=tril(A,-1);4 b1 d4 I1 w3 R5 X9 C9 {
U=triu(A,1);$ W, ^3 W7 s* \4 J+ }' Q
B=-D*(L+U);$ a' V# K; f1 L
f=D*b;4 B, q5 u0 z7 b5 w' t. G& `
J=B*u0+f;
- b7 s* d0 j1 D) v# cwhile norm(J-u0)>=eps& y8 O2 W  G9 z% A0 p( U
x0=J;; T) h5 a3 h) {
J=B*u0+f;3 Y) d2 b+ L4 [
end) ^6 O: F0 X# G( J  |" z
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-8-25 05:23 , Processed in 0.911821 second(s), 73 queries .

    回顶部