QQ登录

只需要一步,快速开始

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

求助,逼近积分

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

3

主题

2

听众

20

积分

升级  15.79%

该用户从未签到

自我介绍
追求完美,强迫症。
跳转到指定楼层
1#
发表于 2010-2-25 08:55 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
u''(x)=f(x),
( i, {2 c& P! r) }2 Y- T$ k边界boundary conditions u(a)=0,u(b)=0." p/ t. S8 w2 D* ]; V" C% f) {
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
; Y; D6 a) q+ ?1 b1 e" h" h+ D# t2 @/ ?' k
maltab code:
4 b9 H) S4 s( {) ?# H
/ Z9 y( q) H5 ^, k; q& I* Yclear all
* n! v, v9 u  Y9 Da=0;
) M3 G  _$ @2 ?b=1;
, i0 J; O. V& D1 A2 j3 In=10;
; h$ F: Q! f% S* O; q7 v5 C. }. R) Qh=(b-a)/(n+1);
6 z( j2 e  Z' G$ ?2 k+ gx_with_boundary=linspace(a,b,n+2)';3 k% l5 l1 J9 R% @
x=x_with_boundary(2:n+1);
5 H% r7 r' E  U8 k1 o# }: D* wA=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
. a$ ]3 d/ f+ q6 X% |: [' u- hrhs=f4(x);, d1 Y) ~9 ?8 H
sol=A\rhs;
! k2 A% r) h; A( J; k7 T3 M8 Gsol_with_boundary_conds=[0;sol;0];- x. a, G' w( i- C! r
plot(x_with_boundary,sol_with_boundary_conds);
1 d' Z! i6 s4 c* a$ J# d$ g/ s, t' r6 b2 ]' k; |# f, |
5 Z0 n9 x  [/ e6 s
另开一文档,f4,
/ @6 V1 X) o) H9 \function y=f4(x)7 D- `2 s" r& t. G- l! ?
y=ones(length(x),1);
4 U1 b# T" {5 z0 _; U. m0 W  H* m
出来后是抛物线。
- o2 q8 o6 V/ i1 X  Q7 u2 e% h" C2 |6 m: J; H' N, O

- O4 S$ F0 L" o/ U" _# E! L; k现在问题是,
) @( ~6 ~5 P# x5 G' U) fu''(x)=sin(2pix);. s! L" w& B+ g6 @% ?( f6 Y
u(-1)=0;! t* i4 J$ I" C
u(1)=0;' K3 J8 J+ O- m$ m1 G9 W+ w% k
范围[-1;1]* _+ ~9 @( X& R1 X

' W/ R) g4 l2 I0 x9 i请问各位要怎么把问题代入原题目中,哪些数值需要修改。
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 21:02 , Processed in 0.370912 second(s), 55 queries .

回顶部