QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2356|回复: 1
打印 上一主题 下一主题

求助,逼近积分

[复制链接]
字体大小: 正常 放大
splelvis 实名认证       

3

主题

2

听众

20

积分

升级  15.79%

该用户从未签到

自我介绍
追求完美,强迫症。
跳转到指定楼层
1#
发表于 2010-2-25 08:55 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
u''(x)=f(x),( h2 b# p% N! E& N0 j9 o
边界boundary conditions u(a)=0,u(b)=0.: Q& ~2 f$ [& j- U7 p5 [; h, i
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);, i0 D4 E& Q" q! w- P

3 Q& h: @" t, Y1 T' }8 G: ^% a% \maltab code:* G$ c& H- i6 y: [' ^; c2 c
& z4 n1 L) ~- C. y* r0 @
clear all3 F+ D) b/ [0 U1 |& r6 x
a=0;* A5 L) Z" M/ s  T
b=1;- @9 l) a4 i) s! D2 d0 N- M
n=10;+ w1 E, ^3 }6 y9 K: Q
h=(b-a)/(n+1);5 {/ O$ W6 V7 k, v* O: O) C) P
x_with_boundary=linspace(a,b,n+2)';. K) Z6 }8 H2 o. n. @
x=x_with_boundary(2:n+1);4 b- y8 Z6 [, A) T7 B+ i* b
A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
6 ?7 z/ t7 S2 C1 ]; urhs=f4(x);
1 |$ u; G6 T6 s+ E% Gsol=A\rhs;/ }5 Z$ i( Z7 ?" c: _6 m) Y( H
sol_with_boundary_conds=[0;sol;0];0 b/ X) o9 j- \. A! D/ v* Y
plot(x_with_boundary,sol_with_boundary_conds);: @/ P; q4 x8 P+ {4 c  t( P
& ?& \  B3 h2 m/ d; j$ T! b/ X. F

9 [& I; d  [) |- {3 ~0 B' j另开一文档,f4,
. J3 I3 g+ `8 g! R4 Efunction y=f4(x)
6 V0 u8 A8 M, n8 {$ W4 W. oy=ones(length(x),1);
0 A/ U4 l) f3 W* x% [( U( Y6 U# D% e
出来后是抛物线。
8 z! K0 q( s* ]% w5 a. m! q& _" K# V$ B# j. O" X! E4 k% c
4 o. l4 J$ j) c8 i: n# y) H
现在问题是,
8 n) _9 `% A8 E& @8 ~  b0 Vu''(x)=sin(2pix);' F) E8 q) x1 a$ K/ L; D+ F
u(-1)=0;1 L- _" C3 O3 S1 W& x3 j
u(1)=0;$ @& m( L" z' E. Z
范围[-1;1]
# @+ e- l' _. y% d+ B) G! ~5 _; F7 \: I# [& _2 }
请问各位要怎么把问题代入原题目中,哪些数值需要修改。
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
splelvis 实名认证       

3

主题

2

听众

20

积分

升级  15.79%

该用户从未签到

自我介绍
追求完美,强迫症。
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

手机版|Archiver| |繁體中文 手机客户端  

蒙公网安备 15010502000194号

Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

GMT+8, 2026-5-25 20:34 , Processed in 0.427401 second(s), 55 queries .

回顶部