| 查看: 496 | 回复: 0 | |||
zwb565055403新虫 (小有名气)
|
[交流]
编程为什么结果出来不对啊
|
 我写的程序如下（虽然能运行出来，结果却不对 大家帮我看下程序是否有问题）： % Example 1 clear;clc %%%%%%%%%%%%%%%%%%%% 方程里的参数设置 %%%%%%%%%%%%%%%%%%%%%%%%%%%% alpha=0.5; beta=0.5; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% L=1; %区间x的长度[0,1] Mx=L/0.0025; %区间x的等分的份数 h=L/Mx; %区间x的划分小区间长度 x=[0:Mx]*h; % x的点值 %上面给出了x的点值，结合t点的点值和函数 u(x,t)的值就可以画出 曲面u(x,t),三维数组 N=1000; %时间t的等分份数 tau=0.0025; % 时间t的划分小区间长度 t=[0:N]*tau; %时间t的点值 0点和1点 101个 %%%%%%%%%%%%%%%%%%%%%%%%%%% 系数矩阵里面的常数定义 r R(i) %%%%%%%%%%%%%%%%%% R=(tau^alpha)*(gamma(2-alpha))/(h*h); % r常数 %计算 r(j) for j=1:Mx-1 r(j)=beta*R/j; % 变数 r_j end r; %%%%%%%%%%%%%%%%%%%%%%%%%%%%分数阶系数矩阵%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% v0=1+2*R-r; %主对角线上元素 A=diag(v0); v1=r-R ; v1=v1(2:end);%下次对角线上的元素 B=diag(v1,-1); v3=-R*ones(1,Mx-2); %上次对角线上的元素 C=diag(v3,1); D=A+B+C; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %非齐次项f(x,t)=2\sqrt(t/pi) sin(pi x)-pi t ((cos pi x)/(2x)-pi sin(pi x)) for k=1:size(t,2) %size(a,2)指行向量元素个数 F(k,  =2*(t(k)^(1/2)/pi^(1/2)).*sin(pi*x) + pi*t(k).*(pi*sin(pi.*x) - cos(pi*x)./(2*x));%第k行向量 时间不变而空间位置变化 end F=F(2:end-1,2:end-1); F=F'; %去掉第一行和最后一行 去掉第一列和最后一列 对应 t=0 x=0 %%%%%%%%%%%%%%%%%%下面用递推关系式得到关于 u 的数表 %%%%%%%%%%%%%%%%%%%%%%%%%%% u(:,1)=inv(D)*((tau^alpha)*gamma(2-alpha)*F(:,1)); % D*u1=u0+ tau^alpha)F1 例一中 u0=0 即初值函数为0函数 %定义 w(k)=(1+k)^(1-alpha)-k^(1-alpha);d(k)=w(k)-w(k+1) for k=1:N w(k)=(1+k)^(1-tau)-k^(1-tau); end for k=1:N-1 d(k)=w(k)-w(k+1); %系数 d(k) end d; %%%%%%%%%%%%%%%%%%%%%%%%%% 单独算出 u2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% u(:,2)=inv(D)*((1-w(1))*u(:,1)+(tau^alpha)*gamma(2-alpha)*F(:,1)); %%%%%求和 \sum d(k)*u^(n-k) 列和 sum(A(:,a:b);2) for n=2:N-2 for k=1:n-1 s(:,k)=d(n-k)*u(:,k); end u(:,n+1)=inv(D)*((tau^alpha)*gamma(2-alpha)*F(:,n+1)+(1-w(1))*u(:,n)+sum(s(:,1:n-1),2)); % u0=0 end h=x(2:end-1); g=t(2:end-1); [h,g]=meshgrid(h,g);%生成网格 mesh(h,g,u');%画出曲面图 xlabel('x') %添加坐标轴 ylabel('t') zlabel('u(x,t)') for k=1:N-1 for j=1:Mx-1 plot3(x(j),t(k),u(j,k)) e(j,k)=u(j,k)-t(k)*sin(pi*x(j)); %误差函数 end end e=max(max(abs(e))) %最大误差 |
回复此楼» 猜你喜欢
2027广东省杰青
已经有8人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有5人回复
-
投稿咨询
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有97人回复
-
时间戳变了，能看出什么问题？
已经有6人回复
-
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
找到一些相关的精华帖子，希望有用哦~
在下编程菜鸟，这是我编的一个c++的程序，运行不了总是有一个错误，求大神！
已经有6人回复
- 提示矩阵维数必须一致，不清楚在哪里编程出现了错误，求指教
已经有5人回复
- 研究结果表明不好，能否发出去文章？
已经有7人回复
- 毕业论文急需要ansys模型，算出来结果不对啊，
已经有3人回复
- 计算圆周率得到的结果错误
已经有10人回复
- 程序错误求教！！！！
已经有3人回复
- 请大侠指教matlab的这个错误提示到底是个神马意思？
已经有12人回复
- Fortran编程过程中遇到错误，求大侠帮忙看下
已经有5人回复
- Fortran提示如下错误：
已经有10人回复
- 怎么做MOFs连粉末都没出来，出来的晶体颜色都不对。。。是什么问题啊？
已经有15人回复
- fortran 程序有write和没有write两种情况计算的结果为什么会不一样?
已经有14人回复
- 大家爱帮我看看这段MATLAB程序为啥出错了
已经有5人回复
- 两次调用同一子程序 第二次调用得出的结果是错误的。。
已经有10人回复
- 急求fortran运行错误原因，在线等
已经有7人回复
- 求助一下程序调试过程出现的错误
已经有4人回复
- 请教 fortran 运行错误的原因
已经有13人回复
- 谭浩强C语言教材，P298页照着输入都不对！求原因
已经有12人回复
- 【求助】DLPOLY计算出来的结果，是否要编程才能分析？
已经有15人回复
- 【求助】有很多错误的程序，找不出来了
已经有18人回复
- 【求助】程序运行后出现的错误 access violation
已经有15人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
