数学建模社区-数学中国

标题: 求助，逼近积分 [打印本页]
作者: splelvis    时间: 2010-2-25 08:55
标题: 求助，逼近积分
u''(x)=f(x),
- T4 a; N) k; n1 r- \' k5 ~边界boundary conditions u(a)=0,u(b)=0.
* y3 i+ ~: G1 w1 X(u(x+h)-2u(x)+u(x-h))/h^2=f(x);4 _1 ~; B2 T# }
) y2 S$ y& s, @; @! q2 d+ A% ~% f
maltab code:+ |/ Y$ i# V  E' ?) K
* C' L* O4 H3 ^9 t0 D, ^
clear all) e" _" Y" M4 Z5 `1 h& z8 d2 l
a=0;; Q& {& m  T1 \
b=1;7 A8 ~0 z6 P" c6 c# ^* t
n=10;* N$ n/ P3 E3 u
h=(b-a)/(n+1);& J/ B1 K+ U- _% I% A; A) j
x_with_boundary=linspace(a,b,n+2)';
6 s! x* Z+ b1 ^9 o/ v/ H2 Fx=x_with_boundary(2:n+1);
2 |9 E3 _( z# Y* d3 q) K; FA=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));% d; {) h; H: Z: b& h4 e5 Q4 ^0 h2 p
rhs=f4(x);
7 w( v# u' [' e+ r* [% xsol=A\rhs;
8 _* f, H9 @7 z0 Rsol_with_boundary_conds=[0;sol;0];: ?/ G* H. [" _7 B# ^! @
plot(x_with_boundary,sol_with_boundary_conds);
2 S1 V3 e5 c6 M, {+ h- R% {5 P2 m0 g7 I7 Y
. u2 ^! f) i; f1 s
另开一文档，f4,, @* p! C( G8 {& z, p4 B
function y=f4(x). o/ ]% d5 q8 b- h6 C2 a6 H( I0 D3 @
y=ones(length(x),1);
5 w# W6 q. n# G$ L' {& {) T" D. ^
- b- l8 ~$ E$ r) x5 `& W出来后是抛物线。
9 |: W0 u: P  ^" N" e
3 d4 F3 q$ t, c. L" ?7 o# Z7 T* x  Y/ C
现在问题是，, j4 r8 n8 O& S0 D6 }2 |
u''(x)=sin(2pix);: F6 E" m) h& a4 H" Z; `* ~3 o2 W
u(-1)=0;
; Y! {% p$ a! l: c  D3 ~u(1)=0;
1 b) e; F& R9 N2 k# H$ n8 t; k# K: n% y范围[-1;1]% W' ^3 b5 \8 C9 S$ a. C4 `, Q4 X
+ q! [2 n$ V. z8 L. y
请问各位要怎么把问题代入原题目中，哪些数值需要修改。
作者: splelvis    时间: 2010-2-26 04:07
没人能帮忙一下吗，
+ ~+ S, S* a% r. f还是有什么资料不明确？



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