数学建模社区-数学中国

标题: 求助,逼近积分 [打印本页]

作者: splelvis    时间: 2010-2-25 08:55
标题: 求助,逼近积分
u''(x)=f(x),# q1 f7 `' U9 |
边界boundary conditions u(a)=0,u(b)=0.) y5 X) _, n6 t2 F5 Z/ N
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
6 e3 T4 Z- f# M% o! k( U3 w* K% r) M+ w% Y+ d) ^8 P. \
maltab code:
$ L$ t7 b: |0 R4 l) U6 g+ C$ q; ~: T, x8 S: ^
clear all
& F; w6 z# E8 a5 N( Va=0;* e4 ~& H4 @, u; ~, D
b=1;% X& v3 m( B  `& j8 h( O
n=10;. P7 g; u" g% b# @1 E
h=(b-a)/(n+1);
# J! s" n' k( L; lx_with_boundary=linspace(a,b,n+2)';& I0 }  ?& l5 }% R5 U
x=x_with_boundary(2:n+1);! T  S/ M! A& a/ O' T7 I
A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));" T9 u4 z4 H5 A& ?7 e
rhs=f4(x);( p* V) u5 Y2 w3 I3 w4 [9 \  v
sol=A\rhs;) S$ `8 b$ \- M* C
sol_with_boundary_conds=[0;sol;0];: `( S" ]( t- H: b/ }. j
plot(x_with_boundary,sol_with_boundary_conds);
/ m: ~7 n- x- D' l. d& |  M# j+ I% m  Z

8 s' ?5 R: z' Y) \( |' E. h另开一文档,f4,+ T" z. U' N, N2 ~( Y
function y=f4(x)
$ C& d- L9 S! o$ D  x2 Zy=ones(length(x),1);: u  F& o" ?9 P: R
, Y8 k! _5 Z1 t0 A  ^1 A% t
出来后是抛物线。" J- Q( n; S/ h. A/ a
( i1 z1 H* _! E5 x/ ^+ B& R

  s% h+ L) \" T0 `- S1 s' Q现在问题是,9 U2 F7 g, i6 y1 o% h& K
u''(x)=sin(2pix);
, G/ C! `+ l9 Du(-1)=0;
! j1 e3 q' w2 S4 Y- {4 ou(1)=0;4 Y. e; b/ ~5 d' N7 Z/ D7 E7 \
范围[-1;1]
% b$ S/ y4 N3 Z: ?) |; O' b- W9 i* i/ r( K+ o: F
请问各位要怎么把问题代入原题目中,哪些数值需要修改。
作者: splelvis    时间: 2010-2-26 04:07
没人能帮忙一下吗,
' V/ r$ W8 c4 a( I# b7 t  b还是有什么资料不明确?




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