用MATLAB编写预估-校正法程序求分数阶常微分方程组数值解

王福昌 2017-10-5 10:49

主函数比较 % 根据文献自己编写 q = 0.5; %分数阶阶数 fdefun = @(t,y) ; %一元微分方程 y0 =0;%初值y0为列向量 h = 2^(-6);%步长 tspan = ; = fdewfc(q,fdefun,tspan,y0,h); figure(1) plot(t,y(1,1:end)) ; xlabel('t'); ylabel('y(t)'); hold on plot(t,t.^2-t+0.5,'r-.') ...

10653 次阅读|没有评论

用MATLAB模拟简单小球跳动

王福昌 2017-9-30 16:06

这里使用了ode求解器的event事件技巧 小球满足的微分方程 \begin{array}{cl} y'_1 = y_2\\ y'_2 = -9.8 . \end{array} 写作 function dydt = f(t,y) dydt = ; 小球在接触地面时会反弹，事件函数可写为 function = bounceEvents(t,y) value = y(1); &n ...

5343 次阅读|没有评论

用MATLAB绘制空间曲面曲面及其切平面

王福昌 2017-9-30 11:20

先定义函数 f = @(x,y) x.^2 + y.^2; 求偏导数 = meshgrid(-5:0.25:5); = gradient(f(xx,yy),0.25); 切点 (x0,y0) = (1,2) . 函数值 f(1,2) = 5 . x0 = 1; y0 = 2; t = (xx == x0) (yy == y0); indt = find(t); fx0 = fx(indt); fy0 = fy(indt); 创建切平面方程 z = ...

个人分类: 教学辅导|10500 次阅读|没有评论

MATLAB优化工具箱遗传算法函数ga()用法示例

王福昌 2017-9-11 22:09

用法格式 X = ga(FITNESSFCN,NVARS,A,b,Aeq,beq,lb,ub,NONLCON,options) minimizes with the default optimization parameters replaced by values in OPTIONS. OPTIONS can be created with the OPTIMOPTIONS function. See OPTIMOPTIONS for details. For a list of opt ...

24141 次阅读|没有评论

MATLAB中的两个常用数值优化函数fminsearch()和fminunc()

热度 1 王福昌 2017-9-8 23:10

MATLAB中提供了两个命令：fminsearch(),fminunc(), x = fminsearch(Fun,x0) = fminsearch(Fun,x0,opt,p1,p1,...) fminunc 与 fminsearch 用法相同。 测试版本 MATLAB 2016b. 例1 已知 ,求其最小值并显示求解过程。 f = @(x)(x(1)^2-x(1))*exp(-x(1)^2-x(2)^2-x(1)*x(2));x0 = ; = fm ...

个人分类: 教学辅导|17195 次阅读|4 个评论 热度 1

无约束最优化问题求解

王福昌 2017-9-8 22:29

无约束最优化问题是最简单的一类最优化问题，其一般数学描述为 其中 称为优化变量， 函数称为目标函数。 1、解析法和图解法 令梯度 得到驻点，即为极值可疑点。 图解法只是适合低维的一元、二元和三元函数。 例1 求一元函数最值 解析解方 ...

个人分类: 教学辅导|5836 次阅读|没有评论

LaTex编辑器WinEDT10.2与TeXLive2016套装配合

热度 2 王福昌 2017-8-29 21:40

由于CTex套装停止更新，于是装了Texlive2016，但是有一些期刊投稿模板用的仍是老的Tex模板，在Texlive2016下采用的新编辑器TexStudio下显示乱码，当然可以通过改为 utf8 格式显示正常，但是毕竟不如Winedt好使。于是单独安装了Winedt10.2，没想到可以直接编译成 pdf 格式。

个人分类: 教学辅导|5849 次阅读|4 个评论 热度 2

基于黄金分割法的最速下降和共轭梯度法与MATLAB程序实现

王福昌 2017-8-23 22:28

教学笔记，留存备用。 黄金分割法搜索极小值 function = goldsearch(fun,a,b,eps,itmax) %黄金分割法搜索极小值 % ---input % fun 所求的目标函数 % a 区间的下界 % b 区间的上界 % eps 区间的最小阈值长度,精度 % ---output % xmin 函数取极 ...

个人分类: 教学辅导|7418 次阅读|没有评论

毕业二十周年聚会校园随拍

热度 2 王福昌 2017-8-12 16:59

图1 毕业二十周年聚会合影 图 2 学而不厌 诲人不倦 &nb ...

4459 次阅读|2 个评论 热度 2

本页有 1 篇博文因作者的隐私设置或未通过审核而隐藏