数学建模社区-数学中国

标题: 求助，逼近积分 [打印本页]
作者: splelvis    时间: 2010-2-25 08:55
标题: 求助，逼近积分
u''(x)=f(x),
) N( G2 H( @  ?. J2 ]) l# h% H边界boundary conditions u(a)=0,u(b)=0.. l+ R  f  b" \0 i% _. H* n
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
( w8 _/ s' i4 d# Q
9 }/ N7 v* m: s( B# qmaltab code:
% I4 x9 N) |4 m$ s, I2 s* d0 v
4 ~4 g9 F( g. L4 ~7 ~$ w8 V- kclear all- \4 B, ^  i6 g; \# T5 ^
a=0;: _! |9 F0 G% X8 L
b=1;" |8 q: a2 j% d7 p0 \, m
n=10;* T, q- ]. U" s* {) C7 U
h=(b-a)/(n+1);2 h- l- ]6 N' p% i) f
x_with_boundary=linspace(a,b,n+2)';1 N( ^( I0 O/ ~0 m+ j8 m8 b$ X
x=x_with_boundary(2:n+1);
! z, v! [2 i: \. f# ?- NA=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));: D0 W1 t2 Q: @9 B% q2 t/ z  v
rhs=f4(x);, z& }/ g1 \9 Y
sol=A\rhs;
4 c$ l2 S8 G; b! f5 isol_with_boundary_conds=[0;sol;0];
7 Q$ |  x3 R+ b# A6 e, [1 lplot(x_with_boundary,sol_with_boundary_conds);( v- U, k0 I3 y$ E+ S
6 v( k9 O# ?8 J& d& C2 q
9 _, U* R5 I& Y, H
另开一文档，f4,- f( e7 W! W+ v/ f/ h' q
function y=f4(x)
: p/ C3 O) g8 H$ p! A8 wy=ones(length(x),1);6 @+ ]9 i6 d. W" z7 v

  v1 A" j% \8 U8 l0 ~出来后是抛物线。
; Q; u0 p* Z8 x8 Y) _! O( ^5 u' O  M

& I1 W( v! ^: {, `2 D  V现在问题是，3 F7 Y' c* q5 J
u''(x)=sin(2pix);" O: G& G; q9 |4 f# [# F
u(-1)=0;
/ r" J+ o1 P- uu(1)=0;
6 a0 U: p7 k" W, y范围[-1;1]9 f- K' J7 E8 M) F

6 _, A* ]8 T/ u( U. y  m请问各位要怎么把问题代入原题目中，哪些数值需要修改。
作者: splelvis    时间: 2010-2-26 04:07
没人能帮忙一下吗，1 ^, k( J( S% Z$ B. f2 d. L0 t
还是有什么资料不明确？



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5