数学建模社区-数学中国
标题:
求助,逼近积分
[打印本页]
作者:
splelvis
时间:
2010-2-25 08:55
标题:
求助,逼近积分
u''(x)=f(x),
! Z0 P3 e5 o+ }% Q* p
边界boundary conditions u(a)=0,u(b)=0.
$ S1 r6 z* ?' Z& K
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
! J8 P$ F5 M/ M6 ^0 N
+ q/ P/ E( f# a$ a! u
maltab code:
8 @2 c" t: @8 o$ _
1 h# l3 C5 a0 Y0 T- q7 Z+ V
clear all
" g* W! c3 g. g* S/ x# n
a=0;
6 ?; K$ U5 j8 V3 M, Q
b=1;
. N( D$ b4 k) y: _6 U% m
n=10;
Z1 |; y0 s0 ^
h=(b-a)/(n+1);
0 a8 Q) X4 {- [! j
x_with_boundary=linspace(a,b,n+2)';
: r% ]$ k8 S4 E* x/ p* `
x=x_with_boundary(2:n+1);
3 Y* p, r4 X; `
A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
/ I8 R7 H9 [" T- l# d/ Y1 o
rhs=f4(x);
# n: k' K0 m J/ s0 V
sol=A\rhs;
! `6 z& `7 @+ w; P; E* {
sol_with_boundary_conds=[0;sol;0];
1 f5 v, v" }7 t" l* T: n: U
plot(x_with_boundary,sol_with_boundary_conds);
% Y& n( `" l' _$ J8 H, m2 h4 c
# H# p8 h% ^5 U) \ M3 u
) ?$ j6 ^; e! `
另开一文档,f4,
* m; n1 _( f( @- N: G; K& u
function y=f4(x)
$ t8 J; c* T. o' Y1 T- J: ~) E
y=ones(length(x),1);
" K R: y9 y; B, J& R
( X# }! C( o, \+ P2 ^
出来后是抛物线。
; }: e9 n4 H0 J/ A
( U" F/ {3 z( z1 z. H7 a& q
3 A3 e8 d/ b1 r+ T- N( `& v
现在问题是,
* ^% N0 X# R+ a% }4 F
u''(x)=sin(2pix);
" t9 l3 [0 J W% @" H
u(-1)=0;
3 e) _7 q" }/ H5 d
u(1)=0;
2 D% P6 Q3 t, m' `, m( u: d& o
范围[-1;1]
& N( w5 D) |( g C( c) E
& m6 t5 N1 L, u; t6 V8 C* O( q: ~. Q
请问各位要怎么把问题代入原题目中,哪些数值需要修改。
作者:
splelvis
时间:
2010-2-26 04:07
没人能帮忙一下吗,
; p8 b7 E8 v: y0 }- R. b! |8 I8 h/ n
还是有什么资料不明确?
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5