楼主，我只计算时间方向的有限差分，为什么得出的图跟实际解还是_matlab辅导吧

06月28日漏签0天

matlab辅导吧关注：10贴子：97

3回复贴，共1页

<返回matlab辅导吧

楼主，我只计算时间方向的有限差分，为什么得出的图跟实际解还是

只看楼主收藏回复

送TA礼物

IP属地:辽宁

1楼2021-07-01 10:11回复

方程发一下看看

IP属地:湖北

来自Android客户端2楼2021-07-01 10:16

clc
clear all
tic;N=16;
T=1;
tau=T/N;
t=0:tau:T;
[T]=meshgrid(t);
F=0.25.*(1-T.^2)./(1+T.^2).^2;
u=zeros(N+1,N+1);
I=eye(N+1);
for kk=2:N+1
u(:,kk)=I\(u(:,kk-1)+tau*F(:,kk));
end
plot(T',u')
figure
uex=0.25.*T'./(1+T'.^2);
plot(T',uex)
err=max(max(abs(u-uex)));