QQ登录

只需要一步,快速开始

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

求助,逼近积分

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

3

主题

2

听众

20

积分

升级  15.79%

该用户从未签到

自我介绍
追求完美,强迫症。
跳转到指定楼层
1#
发表于 2010-2-25 08:55 |只看该作者 |正序浏览
|招呼Ta 关注Ta
u''(x)=f(x),
* z  t4 q) O- ~9 i边界boundary conditions u(a)=0,u(b)=0.3 C# O  n# y" s
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);  ?- z3 q! b( j/ G

  Q' x7 E( G* z' ^$ K) nmaltab code:
; H" t2 ~7 n9 w7 B
. i% E' w! p$ |; U6 D" wclear all4 j! @+ l! B; A) s
a=0;& j# P  L5 C2 C4 s
b=1;
2 a- H9 v. N2 S9 tn=10;2 y) S0 G* H$ @
h=(b-a)/(n+1);) o# C( ^6 k7 \) Y+ O8 G
x_with_boundary=linspace(a,b,n+2)';
( j4 O+ T9 I# i# Ux=x_with_boundary(2:n+1);
& b" J8 [" U" I. Z* _A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));$ i0 |  \, j; ~0 H" Y
rhs=f4(x);
3 E5 r) c; l: z8 bsol=A\rhs;
' [0 ?5 O$ Q1 }* `: R3 T* Isol_with_boundary_conds=[0;sol;0];* f- m5 d! r8 I5 J! X: H/ @2 }
plot(x_with_boundary,sol_with_boundary_conds);+ I: G, E+ T( Y1 ~

) J8 s: s$ O+ n$ N3 u! x
: W) J6 {& o! n* K) L另开一文档,f4,4 ~4 P0 u$ }6 T, N
function y=f4(x)
; _2 C2 e  a/ K% I& o" p: S& ^y=ones(length(x),1);
# W# k+ y% F/ v( i0 h" e; |) z* Y8 y/ F1 g, \7 g0 [- a% O- T4 D
出来后是抛物线。' c) W$ R( N8 M: D

8 [' R" O9 k  K; Q" B8 z& {! A) N& S
现在问题是,. [7 E4 i( t9 ~6 O/ G8 B8 G
u''(x)=sin(2pix);( Y% R! `8 V, y7 K/ \' F+ O
u(-1)=0;5 i6 B( T2 @' p( k- t6 e, i
u(1)=0;1 S9 U6 Y7 g2 h
范围[-1;1]
; x4 V# t" @4 d* w# v; e8 m8 O
+ M% s, d2 e# K9 T4 H请问各位要怎么把问题代入原题目中,哪些数值需要修改。
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 11:00 , Processed in 0.441788 second(s), 56 queries .

回顶部