抛物线显示差分格式matlab,与u=e^(-t)sinx比较,我只给的起那么多分了,
来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/24 17:02:45
抛物线显示差分格式matlab,与u=e^(-t)sinx比较,我只给的起那么多分了,
显示格式的效果比较差,而且对时间步长有限制,不然解就比较差了.下面是matlab的代码,其实主要部分是离散格式,掌握好网格划分与差分格式,其他的都是编程问题了.
%不知道你想比较什么,直接运行出图吧
N = 10;%建议取值10与30进行对比,取10时网格比较粗糙但效果可以,取30时结果发散.
%希望不要追问隐式格式的,那个又是另一种了.
T = 1;
X = pi;
dt = T/(N-1);
dx = X/(N-1);
t = 0:dt:T;
x = 0:dx:X;
co = dt/dx/dx;
u = zeros(N);
u(1,:) = sin(x);
u(:,1) = 0;
u(:,N) = 0;
for i = 2:N
for j = 2:N-1
u(i,j) = co*(u(i-1,j-1)+u(i-1,j+1)-2*u(i-1,j))+u(i-1,j);
end
end
v = zeros(N);
for i = 1:N
for j = 1:N
v(i,j) = exp(-t(i))*sin(x(j));
end
end
figure();
surf(x,t,u);
figure();
surf(x,t,v);
再问: 我想要的是数据的比较,发我邮箱吧,942069078@qq.com
再答: 你要的是分析结果吗?还是什么?上面程序中的u就是计算的近似值,v就是准确解,想比较就用 abs(u)-abs(v) 看看这个矩阵的最大值和最小值就是你想要的了 如果按误差公式就是求不同的范数就行了 norm(u-v,?); 主要问题是你想要用哪种方法算误差。
再问: 我也不知道老师要我们做什么的,能不能一堆u值,一堆V值,再加上他们的绝对值
再答: 估计就是比较数值解与精确解的误差,如果是一个未知数当然就是差的绝对值比较,现在是矩阵当然是按矩阵范数为准。范数是等价的,所以你选哪个是无所谓的。就直接norm(u-v)/norm(v) 表示个相对还是绝对误差?我也记不得了。就这样了。交上去应该没问题。主要就是想让你知道离散的过程和显示格式的缺点。
%不知道你想比较什么,直接运行出图吧
N = 10;%建议取值10与30进行对比,取10时网格比较粗糙但效果可以,取30时结果发散.
%希望不要追问隐式格式的,那个又是另一种了.
T = 1;
X = pi;
dt = T/(N-1);
dx = X/(N-1);
t = 0:dt:T;
x = 0:dx:X;
co = dt/dx/dx;
u = zeros(N);
u(1,:) = sin(x);
u(:,1) = 0;
u(:,N) = 0;
for i = 2:N
for j = 2:N-1
u(i,j) = co*(u(i-1,j-1)+u(i-1,j+1)-2*u(i-1,j))+u(i-1,j);
end
end
v = zeros(N);
for i = 1:N
for j = 1:N
v(i,j) = exp(-t(i))*sin(x(j));
end
end
figure();
surf(x,t,u);
figure();
surf(x,t,v);
再问: 我想要的是数据的比较,发我邮箱吧,942069078@qq.com
再答: 你要的是分析结果吗?还是什么?上面程序中的u就是计算的近似值,v就是准确解,想比较就用 abs(u)-abs(v) 看看这个矩阵的最大值和最小值就是你想要的了 如果按误差公式就是求不同的范数就行了 norm(u-v,?); 主要问题是你想要用哪种方法算误差。
再问: 我也不知道老师要我们做什么的,能不能一堆u值,一堆V值,再加上他们的绝对值
再答: 估计就是比较数值解与精确解的误差,如果是一个未知数当然就是差的绝对值比较,现在是矩阵当然是按矩阵范数为准。范数是等价的,所以你选哪个是无所谓的。就直接norm(u-v)/norm(v) 表示个相对还是绝对误差?我也记不得了。就这样了。交上去应该没问题。主要就是想让你知道离散的过程和显示格式的缺点。
抛物线显示差分格式matlab,与u=e^(-t)sinx比较,我只给的起那么多分了,
求一篇语文试卷分析急,我只考了76分,就连平时比我差的都考了80多分,:'(我等会就要,求大神帮我,写的好了再给分.急.
翻译一句很简单的..我英语太差了 ..所有分都给了..
我的名字叫周正、谁给我起个英文名.没分了.
E.T.的全称是?我是给分的,
急 帮我看一下这个分段差分方程作图哪里错了,matlab
120的题只考了90分;语文差3分达到老师给我的合格目标(80分),要做很多作业TAT英语89= =我对物理有一点兴趣,
Matlab里如何只显示矩阵中的一个圆形区域?【200分求解答】
我的一篇作文,老师只给了46分(满分六十),
可以手打和给我其他比较相符的,我看过的出现了不给分
求raining on sunday 原版伴奏 Keith urban的!我就这么点分了 差不过都给了
关于社会活动调查社会活动调查应该怎么写,要做些什么活动?格式要怎么样?最好给我一篇比较经典的范文,回答得好的有加分