clear;clc;
y0 = [0 1];
h = 0.1;
tspan = [0 10];
N = 1/h +1 ;
alphat = [0.01 0.99];
ah = 0.01;
MM = (alphat(2) -alphat(1))/ah ;
T = [];
global a
for a=alphat(1):ah:alphat(2)
varphi = sqrt(3)/pi*log(a/(1-a));
[y,~] = RepMilneHammin(@example2,tspan,y0,h);
[yy,~] = MilneHammin(@example2,tspan,y0,h);
[yyy,~] = RungeKutta4(@example2,tspan,y0,h);
[yyyy,t] = AdamsSimpson(@example2,tspan,y0,h);
M = length(t);
z = (exp(t+varphi*t)-1)/(1+varphi);
T = [T;y(N,1),yy(N,1),yyy(N,1),yyyy(N,1),z(N,1)];
% e = sum((y(:,1)-z).^2)/M;
% ee = sum((yy(:,1)-z).^2)/M;
% eee = sum((yyy(:,1)-z).^2)/M;
% eeee = sum((yyyy(:,1)-z).^2)/M;
% [e,ee,eee,eeee]
%
% figure
% plot(t,y(:,1),t,yy(:,1),t,yyy(:,1),t,yyyy(:,1),t,z);
end
alpha = alphat(1):ah:alphat(2);
alpha = alpha';
plot(T(:,1),alpha,'-o',T(:,2),alpha,'-.',T(:,3),alpha,'-*',T(:,4),alpha,'-x',T(:,5),alpha);
legend('RepMilneHammin','MilneHammin','RungeKutta4','AdamsSimpson','Exact solution')
xlabel('x_{t}^{\alpha}','FontName','times new Roman','FontSize',15)
ylabel('\alpha','FontName','times new Roman','FontSize',15);
MSE = [sum((T(:,1)-T(:,5)).^2)/MM,sum((T(:,2)-T(:,5)).^2)/MM,sum((T(:,3)-T(:,5)).^2)/MM,sum((T(:,4)-T(:,5)).^2)/MM];
MAE = [sum(abs(T(:,1)-T(:,5)))/MM,sum(abs(T(:,2)-T(:,5)))/MM,sum(abs(T(:,3)-T(:,5)))/MM,sum(abs(T(:,4)-T(:,5)))/MM];
MSE
MAE
MATLAB实现求解常微分方程数值解(完整源码+数据).zip
版权申诉
97 浏览量
2023-10-16
21:14:12
上传
评论 2
收藏 24KB ZIP 举报
不安分的小女孩
- 粉丝: 9213
- 资源: 2001
最新资源
- Flume进阶-自定义拦截器jar包
- Dubins曲线算法讲解和在运动规划中的使用.pdf
- 上市公司-股票性质数据-工具变量(民企、国企、央企)2003-2022年.dta
- 上市公司-股票性质数据-工具变量(民企、国企、央企)2003-2022年.xlsx
- Reeds+Shepp曲线算法讲解和实现.pdf
- 毕业设计基于SpringBoot+MyBatisPlus+MySQL+Vue的外卖配送信息系统源代码+数据库
- 词向量(Word Embeddings)是自然语言处理(NLP)领域的一种重要技术.txt
- Surfer,线性函数
- MyBatis 的动态 SQL 是其核心特性之一.txt
- 时代的sdddsddsddsd
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈