QQ登录

只需要一步,快速开始

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

求助,逼近积分

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

3

主题

2

听众

20

积分

升级  15.79%

该用户从未签到

自我介绍
追求完美,强迫症。
跳转到指定楼层
1#
发表于 2010-2-25 08:55 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
u''(x)=f(x),
+ B4 |+ Z. j' g1 M5 v边界boundary conditions u(a)=0,u(b)=0.
6 X  _2 i& Z, |  X6 _1 H* P(u(x+h)-2u(x)+u(x-h))/h^2=f(x);7 l) @$ S. `6 `0 ?4 W9 O

% G5 R( j  ]  x- Y  Hmaltab code:' L8 ]1 N1 ~4 r3 k; ?1 H6 }5 E

& B1 S# a+ W/ y7 ]6 l$ b3 }# \clear all
( T0 Y; x* D! M3 a3 A1 Ya=0;0 z  f3 ?  r1 c6 q4 W6 E
b=1;
# |( k3 }  t. @1 ^# M$ H5 M8 P9 w: [n=10;, @  {1 I: v2 I3 g7 o
h=(b-a)/(n+1);! n8 R* A1 R( i
x_with_boundary=linspace(a,b,n+2)';
5 M4 D& X  y; Y- i3 |7 z! Ix=x_with_boundary(2:n+1);- Q4 q" {$ `* Y2 L
A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
, g! g( d$ m, J" E4 ^" R( d1 jrhs=f4(x);  B1 g9 d0 @; u; U  Y! ~
sol=A\rhs;3 I4 K2 T' K6 R$ ?) `* D) U$ t* ~
sol_with_boundary_conds=[0;sol;0];# E' G( u6 }- m3 G! S$ a2 m  \' o' w
plot(x_with_boundary,sol_with_boundary_conds);3 v+ V- p1 |! O* x

' @' _. W2 T! ]
* E6 j; P$ i+ r. }另开一文档,f4,
; C- c, O4 o9 T3 Ofunction y=f4(x)4 T1 B! {- w3 J4 \
y=ones(length(x),1);8 F" f9 a: {* z( q0 v- v0 {
& S4 _9 ^- F, x' R! E3 Y
出来后是抛物线。
, W9 K! v4 R$ D  t( K+ B. l' k3 s8 y0 G- b

) {7 f; D+ b. z0 z现在问题是,( Z5 P; b+ T9 N3 k7 w  i8 [0 K
u''(x)=sin(2pix);$ @- I8 v  G7 S! a: A0 G, U
u(-1)=0;
/ g( C) s. \+ D% N) ~2 {u(1)=0;+ E& t; ?9 i; u' ^
范围[-1;1]% ~8 F% K3 N" i

* a+ D! }* l/ T7 s( I* Y请问各位要怎么把问题代入原题目中,哪些数值需要修改。
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-7-10 14:14 , Processed in 0.405301 second(s), 56 queries .

回顶部