数学建模社区-数学中国
标题:
求助，逼近积分
[打印本页]
作者:
splelvis
时间:
2010-2-25 08:55
标题:
求助，逼近积分
u''(x)=f(x),
! C5 F) B9 a6 z* n
边界boundary conditions u(a)=0,u(b)=0.
9 o/ L1 S! \8 `. |4 A2 l k {
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
9 x# o" B* S% i+ p. q
' a7 d% o- {3 A6 y) m
maltab code:
2 ?' r5 e, k1 ~. b
1 p" x* k# B5 d0 M9 x
clear all
, K4 g+ l. F$ q8 @( U/ o
a=0;
- U9 [ |; [# b
b=1;
5 W& g0 T) m0 D
n=10;
2 a" ^+ I5 \' [) u
h=(b-a)/(n+1);
- N! Q: G& Y% ?+ w& S% O8 d
x_with_boundary=linspace(a,b,n+2)';
! v" D3 Z0 q- V. R. K5 u
x=x_with_boundary(2:n+1);
# w7 A. @% R+ l8 w% |
A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
3 n. R4 |$ w5 l! }
rhs=f4(x);
% F; `' L k1 T; |4 L
sol=A\rhs;
. G. ^" z4 m1 Q
sol_with_boundary_conds=[0;sol;0];
0 F5 d. h/ c9 C6 T u
plot(x_with_boundary,sol_with_boundary_conds);
& z' o$ m6 H. L- R9 }9 O
# C6 K1 E% I* @9 C3 p3 s
+ i" L% y' W. @
另开一文档，f4,
7 F) i l7 o- A4 E ]
function y=f4(x)
$ p. d- R! c/ g, }! R
y=ones(length(x),1);
1 @. {3 l& N" P/ H. r
: c! |5 F! R( `
出来后是抛物线。
6 }/ l7 c' I$ ^ ?- I* Y
( n# M2 U# a4 v, t
t# Q" Z9 u# y' H7 i' L9 d I
现在问题是，
8 O+ d( C: R' z5 S# }0 c
u''(x)=sin(2pix);
! P8 y+ b4 ^5 m& d9 B V4 a
u(-1)=0;
& X. h9 e: ]. \5 }/ K9 H& X' S# p
u(1)=0;
" l2 A) O2 t) W9 H' t) j
范围[-1;1]
3 b S0 \4 W# F) `, G' n' k+ N, c) `
) _2 H0 u; Y D( w! M K
请问各位要怎么把问题代入原题目中，哪些数值需要修改。
作者:
splelvis
时间:
2010-2-26 04:07
没人能帮忙一下吗，
5 N7 P- x, v2 _1 }% d& K
还是有什么资料不明确？
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5