QQ登录

只需要一步,快速开始

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

求助,逼近积分

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

3

主题

2

听众

20

积分

升级  15.79%

该用户从未签到

自我介绍
追求完美,强迫症。
跳转到指定楼层
1#
发表于 2010-2-25 08:55 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
u''(x)=f(x),! ]/ B+ l' l# v
边界boundary conditions u(a)=0,u(b)=0.
) k/ G. x+ q8 {2 s$ l7 C2 }(u(x+h)-2u(x)+u(x-h))/h^2=f(x);* Z# ?# v2 x9 ]. i5 G7 u& b! _" b

$ m$ S! v9 d! V& t5 ?& }maltab code:
* t, G0 ~- h2 ?- Z; r* C* ~
8 [) j. D$ O5 d0 G/ y9 z7 ]7 @4 Hclear all' ]3 }# F. r$ L' \) N. [
a=0;
4 U$ z2 `  B0 o, H# T3 ]  x: Ub=1;+ w# c( g( }5 e1 }$ ]7 k
n=10;" `0 [$ i/ X3 d, G: d, S
h=(b-a)/(n+1);0 E2 K- I3 ?: e) a3 V" J5 o
x_with_boundary=linspace(a,b,n+2)';
, c- k. {1 t0 }% Mx=x_with_boundary(2:n+1);! X1 \8 O" }3 {. a
A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
! T! M% g2 h7 B' H8 v1 brhs=f4(x);
% \, v( M: J/ d. Qsol=A\rhs;
+ j9 N2 e! o* a  Z' A% psol_with_boundary_conds=[0;sol;0];
- M& k$ |3 m8 ]3 d% M' \plot(x_with_boundary,sol_with_boundary_conds);
* J' P, k5 e% f  {
' ?  L* L2 X" o9 e, H, m. B2 ~- c2 t( m  S
另开一文档,f4,
8 D# E$ l! m% {0 rfunction y=f4(x)1 x' `$ r( W5 R# }! x7 e3 x
y=ones(length(x),1);
, n/ M: |# E7 c9 Z
, W& u: z2 c- _, t- K+ r# G出来后是抛物线。
* g. H" \( A! `: N4 h/ _$ o8 Z( A7 p) |2 n8 g! l
& b1 `; M& b+ O3 k8 z8 i
现在问题是,
5 _: r& @  |) l& k( }$ z" n3 Uu''(x)=sin(2pix);
/ h4 z1 h6 M; q& ?u(-1)=0;
4 A" N8 \8 F& ]3 _u(1)=0;: @2 P! [, s3 S7 l  h
范围[-1;1]
& l! I* h9 m6 S9 Q! W( j3 S+ Q
, _% E0 `5 M, _! A' i2 e请问各位要怎么把问题代入原题目中,哪些数值需要修改。
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-26 06:13 , Processed in 0.397970 second(s), 56 queries .

回顶部