标题: 求助，逼近积分 [打印本页] 作者: 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还是有什么资料不明确？