QQ登录

只需要一步,快速开始

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

求助,逼近积分

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

3

主题

2

听众

20

积分

升级  15.79%

该用户从未签到

自我介绍
追求完美,强迫症。
跳转到指定楼层
1#
发表于 2010-2-25 08:55 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
u''(x)=f(x),
1 V& r# d& u1 J2 h( O  R边界boundary conditions u(a)=0,u(b)=0.
5 _( C) X% i" \( ]. f( w% a(u(x+h)-2u(x)+u(x-h))/h^2=f(x);
. [; V! [* k9 b+ G" E' T
* X# T7 L8 I/ a0 V, Y# jmaltab code:. Q0 C3 ^3 x- @8 w+ s  [
' z( v% {! i  a( e( K9 x
clear all, f5 U* t, }7 x4 T# g$ W7 ]( I8 f
a=0;
' {6 Z$ [" E. Q' l  i8 e# Pb=1;
, O2 R) J2 j, R% n4 U) on=10;
: p! Y, j/ B$ g7 I, k1 dh=(b-a)/(n+1);" ^$ d" {1 ~1 J7 [( e& r8 e: f
x_with_boundary=linspace(a,b,n+2)';* W9 Y' \+ p5 A) U! T5 `
x=x_with_boundary(2:n+1);. Q9 f, y6 o5 A) u5 e4 o+ k0 S
A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));- s  f" a* ?- h
rhs=f4(x);
, V+ v' K. i7 t# z7 lsol=A\rhs;
' W! F3 T" T" I. vsol_with_boundary_conds=[0;sol;0];
. i( D9 h  F( Y! S+ `plot(x_with_boundary,sol_with_boundary_conds);
# v: O7 b3 l, S7 `; i2 r. ~! j! O+ h( o) v& a

( H) A0 h& i* {* I/ B8 U* j. N另开一文档,f4,9 o) R$ Q+ w0 e! p4 `1 I. L
function y=f4(x)& s6 h9 y- ]) M
y=ones(length(x),1);
% A: ^: U6 b+ _/ E$ ?/ y: U$ m4 U
出来后是抛物线。
. U) [. _$ C4 m9 f4 y6 L( b+ |
, ]* y2 r& s, k; v( j5 F; ^1 W' g, N6 c
! H( f9 H2 @& Q- V! O" W现在问题是,  j; z6 G  d3 [% G" m
u''(x)=sin(2pix);
) G$ ], {7 r, k8 v- `5 R1 ru(-1)=0;. O0 M$ L5 k" L" z
u(1)=0;
) A8 ]& p5 `7 O/ `4 e3 K& {7 ]范围[-1;1]0 B; N8 u! K  Y  C' E4 a" K

/ p+ W- k5 f3 s0 ^5 m$ u  x请问各位要怎么把问题代入原题目中,哪些数值需要修改。
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-8-25 06:58 , Processed in 0.442118 second(s), 56 queries .

回顶部