- 在线时间
- 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),
1 V& r# d& u1 J2 h( O R边界boundary conditions u(a)=0,u(b)=0.
5 _( C) X% i" \( ]. f( w% a(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
. [; V! [* k9 b+ G" E' T
* X# T7 L8 I/ a0 V, Y# jmaltab code:. Q0 C3 ^3 x- @8 w+ s [
' z( v% {! i a( e( K9 x
clear all, f5 U* t, }7 x4 T# g$ W7 ]( I8 f
a=0;
' {6 Z$ [" E. Q' l i8 e# Pb=1;
, O2 R) J2 j, R% n4 U) on=10;
: p! Y, j/ B$ g7 I, k1 dh=(b-a)/(n+1);" ^$ d" {1 ~1 J7 [( e& r8 e: f
x_with_boundary=linspace(a,b,n+2)';* W9 Y' \+ p5 A) U! T5 `
x=x_with_boundary(2:n+1);. Q9 f, y6 o5 A) u5 e4 o+ k0 S
A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));- s f" a* ?- h
rhs=f4(x);
, V+ v' K. i7 t# z7 lsol=A\rhs;
' W! F3 T" T" I. vsol_with_boundary_conds=[0;sol;0];
. i( D9 h F( Y! S+ `plot(x_with_boundary,sol_with_boundary_conds);
# v: O7 b3 l, S7 `; i2 r. ~! j! O+ h( o) v& a
( H) A0 h& i* {* I/ B8 U* j. N另开一文档,f4,9 o) R$ Q+ w0 e! p4 `1 I. L
function y=f4(x)& s6 h9 y- ]) M
y=ones(length(x),1);
% A: ^: U6 b+ _/ E$ ?/ y: U$ m4 U
出来后是抛物线。
. U) [. _$ C4 m9 f4 y6 L( b+ |
, ]* y2 r& s, k; v( j5 F; ^1 W' g, N6 c
! H( f9 H2 @& Q- V! O" W现在问题是, j; z6 G d3 [% G" m
u''(x)=sin(2pix);
) G$ ], {7 r, k8 v- `5 R1 ru(-1)=0;. O0 M$ L5 k" L" z
u(1)=0;
) A8 ]& p5 `7 O/ `4 e3 K& {7 ]范围[-1;1]0 B; N8 u! K Y C' E4 a" K
/ p+ W- k5 f3 s0 ^5 m$ u x请问各位要怎么把问题代入原题目中,哪些数值需要修改。 |
zan
|