实验名称微分方程与差分方程。
计算机与信息安全学院计算机科学与技术专业。
班第实验小组。
作者学号。同作者。
实验日期 2017 年 1 月 3 日。
一、 实验目的。
1. 掌握用matlab求微分方程和微分方程组的数值解的方法。
2. 进一步巩固加强差分方程模型的建模、求解能力;
3. 掌握用matlab求解差分方程数值解的相关命令。
二、 实验内容。
1. 微分方程及方程组的解析求解法;
2. 微分方程及方程组的数值求解法;
3. 直接使用matlab命令对微分方程(组)求解(解析解、数值解);
4. 利用图形对解的特征作定性分析;
5. 差分方程模型的建立与求解。
三、 实验任务。
1. 已知微分方程dy/dx=-2y/x+4x,y(1)=2
1) 求其解析解。
解:dsolve('dy=-2*y/x+4*x','y(1)=2','x')
ans =1/x^2 + x^2
2) 求其在[1,3]区间内的数值解并作图。
解:.m文件:
function dy = fun4(x,y)
dy = x^2+1/(x^2+1);
命令:[x,y] =ode45('fun4',[1 3],1.5);
plot(x,y,'-
结果:2. 用matlab求微分方程组在区间[0,2]的数值图形解。
m文件:function dy = fun2(t,y)
dy = zeros(4,1);
dy(1) =y(3)-3*y(1);
dy(2) =y(3)+2*y(4)-7*y(1)*y(2)-5*y(2)^2;
dy(3) =y(3)+2*y(1)*y(2);
dy(4) =y(4)+6*y(2)^2;
命令:[t,y] =ode45('fun5',[0 2],[1 1 0 0]);
plot(t,y(:,1),'t,y(:,2),'t,y(:,3),'t,y(:,4),'x')
结果:3. 下表是某公司的一则房产广告,不难算出,你向银行总共借了25.2万元,30年内总共要还51.696万,约为当初借款的两倍。贷款的年利率是多少?
解:设r为月利率。
clear
s = 252000;
m = 516960;
p = 360;
eq = r) s*(1+r)^p-m;
options = optimset('maxfunevals',1000000,'maxiter',100000);
r = fsolve(eq,0,options);
r结果:r=0.002
年利率为:12*r=2.4%
数学建模 实验报告
解 根据题意我们可以建立目标函数和约束条件。目标函数 约束条件 均大于0且为整数。所以周一开始工作到从周日开始工作的雇员数分别为 2 人 2 人 4 人 3 人 3 人 0 人 8 人时,总聘用费用最少为 11000 元。直接将区间范围改为 5,5 时,结果如下 最小值点取在x 1.5715处,然而...
数学建模实验报告
桂林电子科技大学2017 2018学年第1学期。数学建模。一 实验目的。学会使用matlab软件解决线性规划问题的最优值求解问题。学会使用matlab软件解决非线性规划问题的最优值求解问题。学会将实际问题归结为线性规划问题或者非线性规划问题用matlab软件建立巧当的数学模型来求解。二 实验内容。用...
数学建模实验报告
湖南科技学院实验报告。实验目的。1 学会操作matlab。2 应用软件建立层次分析模型,选择最佳方案。实验内容 算法 程序 步骤和方法 输入 输出 实验结果 实验结果。分析 问题 某企业由于生产效益较好,年底取得一笔利润,领导决定拿出一部分资金分别用于 1 为企业员工发年终奖金 2 扩建集体福利设施...