求助，逼近积分

[复制链接]

字体大小: 正常放大

splelvis

3 主题	2 听众	20 积分

升级 15.79%

该用户从未签到

自我介绍: 追求完美，强迫症。

电梯直达

1^#

发表于 2010-2-25 08:55 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

u''(x)=f(x),
边界boundary conditions u(a)=0,u(b)=0.
(u(x+h)-2u(x)+u(x-h))/h^2=f(x);

maltab code:

clear all
a=0;
b=1;
n=10;
h=(b-a)/(n+1);
x_with_boundary=linspace(a,b,n+2)';
x=x_with_boundary(2:n+1);
A=h^(-2).*(diag(ones(1,n-1),-1)+diag(-2.*ones(1,n),0)+diag(ones(1,n-1),1));
rhs=f4(x);
sol=A\rhs;
sol_with_boundary_conds=[0;sol;0];
plot(x_with_boundary,sol_with_boundary_conds);

另开一文档，f4,
function y=f4(x)
y=ones(length(x),1);

出来后是抛物线。

现在问题是，
u''(x)=sin(2pix);
u(-1)=0;
u(1)=0;
范围[-1;1]

请问各位要怎么把问题代入原题目中，哪些数值需要修改。

zan