数学建模社区-数学中国
标题:
求助，逼近积分
[打印本页]
作者:
splelvis
时间:
2010-2-25 08:55
标题:
求助，逼近积分
u''(x)=f(x),
3 m- `/ y# a1 ]4 X; B3 n
边界boundary conditions u(a)=0,u(b)=0.
( i; A- n) W; N% Q
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
# ~' M& K- |2 ~6 b' a
3 [! B/ B0 ]8 L" e. Q" ]' b& T
maltab code:
; B' S- h' v* c; l+ v, v9 x/ K" m1 s
4 P# i. d( v; l1 Y
clear all
& o1 k' Y4 m1 J+ B
a=0;
3 R$ c) _7 S/ ?& t# b
b=1;
/ _' D, R( T' e
n=10;
; y) `1 K; @2 f# v' h! Y
h=(b-a)/(n+1);
& P1 t; z- }0 y X4 m' X4 W
x_with_boundary=linspace(a,b,n+2)';
' s& E( U$ T& [3 f
x=x_with_boundary(2:n+1);
4 Q/ N K* M6 L
A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
% a, f" l; ~' C3 s) r' O
rhs=f4(x);
6 W) g/ t" s- |
sol=A\rhs;
( k5 ?* S' X5 Y: s& D
sol_with_boundary_conds=[0;sol;0];
5 a2 \+ I8 l$ v& j- ^
plot(x_with_boundary,sol_with_boundary_conds);
^0 K8 c( E0 p G3 Z7 b
0 |- c7 X- j" J" j( t; z: n- T, }7 V
/ [! l$ k) P0 `
另开一文档，f4,
+ \, E" e) o/ z5 L# d4 }/ }
function y=f4(x)
7 r/ V) u: ~$ _+ c4 r3 S* {+ p
y=ones(length(x),1);
1 r) E, O8 `6 g0 c4 b8 P
' _* ?: J2 [$ w% R( y7 Q6 n2 B) g& r$ i
出来后是抛物线。
/ W! I7 i+ A$ z v6 `
. K# N+ S H S1 ^' J
+ U7 R7 }8 ^/ ?" u- R5 \9 N
现在问题是，
- A: w: h- `" v0 L
u''(x)=sin(2pix);
3 n' U+ r3 D1 U: L. M
u(-1)=0;
) t3 V* `4 I' T0 o$ o+ [% \
u(1)=0;
; {, V! [, Q8 W* O; s
范围[-1;1]
) c+ g9 |+ g7 }, D, i3 T7 L/ L$ G9 R
! S2 ^5 h& J% N% ?. N
请问各位要怎么把问题代入原题目中，哪些数值需要修改。
作者:
splelvis
时间:
2010-2-26 04:07
没人能帮忙一下吗，
}& e; {( Q0 Z0 p' n
还是有什么资料不明确？
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5