- 在线时间
- 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),
$ P( |0 n2 L$ A+ [边界boundary conditions u(a)=0,u(b)=0.
7 B% j; I$ F) M' E1 _(u(x+h)-2u(x)+u(x-h))/h^2=f(x);4 P. ~, ^- L3 U5 [
|: W0 P' {0 _9 ?* Rmaltab code:
: m. q, H* I+ X+ W1 Y4 o. r. O% w( s/ q- t( |# `
clear all6 c% ^7 k; A' ]3 a/ g- `) ~* H
a=0;
& D# [/ Z/ e% c& q" w' w0 ^" D+ wb=1;. J) z4 ]: J8 B. h C; Q
n=10;
% N, c2 e8 {1 D1 ]8 h- L+ Dh=(b-a)/(n+1);
: i% s0 W" {6 q6 {6 C$ S1 m5 Gx_with_boundary=linspace(a,b,n+2)';
1 d# T+ D1 T) R9 H% v: Nx=x_with_boundary(2:n+1);
1 y2 a1 Z1 A* w- b8 BA=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
$ ?2 l. C% {) ]$ Prhs=f4(x);
; o: k }7 @+ usol=A\rhs;: M+ t& ^9 P- a7 \! |/ u
sol_with_boundary_conds=[0;sol;0];/ e' R1 F5 U/ [ B( s. y
plot(x_with_boundary,sol_with_boundary_conds);
7 V: j. ?- L( e, F$ @. @- |$ G
+ c6 B. s* H+ @' F6 t" u) Q5 ~0 L& n* ]3 e4 w
另开一文档,f4,
% R: t% @* i0 u$ x( ?- }5 H0 qfunction y=f4(x), Z. p9 B" S5 d$ ]% T' s. T
y=ones(length(x),1);
1 M$ {9 Q5 \# Q0 ~1 i
' O. v" ^& |! j出来后是抛物线。
7 N1 P8 {: r" m+ Y
P/ g3 X; K7 ^9 k: X
' {, ^) o5 Q8 _$ B8 Z. h现在问题是,1 X3 m$ l: F5 k' S
u''(x)=sin(2pix);; {: _: f# I+ S& O* I) S- R4 O* J
u(-1)=0;
/ Z1 C2 q# s4 eu(1)=0;# D+ W2 ~+ N ?
范围[-1;1]
# z1 i/ ?, q3 A4 n. O" f) y' [* s( s7 u( u" Y& A- R
请问各位要怎么把问题代入原题目中,哪些数值需要修改。 |
zan
|