- 在线时间
- 0 小时
- 最后登录
- 2010-3-3
- 注册时间
- 2010-2-19
- 听众数
- 2
- 收听数
- 0
- 能力
- 0 分
- 体力
- 52 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 20
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 8
- 主题
- 3
- 精华
- 0
- 分享
- 0
- 好友
- 0
升级   15.79% 该用户从未签到
- 自我介绍
- 追求完美,强迫症。
 |
u''(x)=f(x), T- u( k9 Z$ b# E& H
边界boundary conditions u(a)=0,u(b)=0.( T( x' f+ H9 e) I! \3 A& u
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
! l2 O" ?* L' b1 t/ l T7 x' ~
: e1 e, K6 L# q$ N, u% B" bmaltab code:
3 h2 L% @2 t. q6 B5 ^& [4 R$ l( I/ \
7 z$ I7 K! g8 bclear all1 f% q# |2 {/ O4 O' m
a=0;
# z& v7 w0 j& b5 hb=1;
6 ^5 ^5 e6 t- E. s" Nn=10;. f E! |3 f3 i- e1 F; P
h=(b-a)/(n+1);/ u1 ]5 w3 \. P
x_with_boundary=linspace(a,b,n+2)';& H: I0 @1 F' I; V
x=x_with_boundary(2:n+1);
! S( {3 E: N1 E1 x7 lA=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
4 N2 Y# V0 Z. _rhs=f4(x);
4 G2 ]. ~3 v# G8 }4 Tsol=A\rhs;
* G1 Y+ }, P esol_with_boundary_conds=[0;sol;0];
' v$ h' d( ?6 r( q- g5 O$ @! aplot(x_with_boundary,sol_with_boundary_conds);
$ [& F4 p6 z4 c p& K D) E) \( l
' p) P& \4 m- x4 E
另开一文档,f4,
3 p; I+ Q; J) cfunction y=f4(x)
$ G5 V6 X/ P9 D) ^! Q0 S5 }0 Zy=ones(length(x),1);
* t+ V N1 A$ [! q$ E; w& |" J. |
! c0 B8 z! q o! M( q出来后是抛物线。5 } ]; `4 t* T& q. Y/ A
# m6 I4 ]" Z& E3 k
/ A* R F- @( q- |: B9 X- G9 S现在问题是,
: q! W/ f% @7 d2 h4 Ru''(x)=sin(2pix);
6 _" n% s; b+ bu(-1)=0;
9 i1 p; z0 M; L8 pu(1)=0;6 {3 \! M# a w, x6 b' |2 _6 w
范围[-1;1]+ q/ I! n8 r8 U5 a6 T, P: A( {* F+ @
+ f( ^- r' g& u2 N
请问各位要怎么把问题代入原题目中,哪些数值需要修改。 |
zan
|