matlab ignoreanalyticconstraints,MATLAB函数随笔之计算篇

前言：这些都是我在学习MATLAB与数学建模时遇到的一些常用函数，就在这里写了下。当然MATLAB十分强大，函数众多，这些只是凤毛麟角，这篇博客我也会一直更新有新的函数就接着写上，希望对大家有所帮助。

前言：这些都是我在学习MATLAB与数学建模时遇到的一些常用函数，就在这里写了下。当然MATLAB十分强大，函数众多，这些只是凤毛麟角，这篇博客我也会一直更新有新的函数就接着写上，希望对大家有所帮助。有不足之处还望指出。

这里是MATLAB画图篇。

数值计算篇

1. 微分方程

解上式x(t)=?

x=dsolve("Dx=r*x*(1-x/m)","x(0)=x0","t")

x =

-m/(exp(log(1 - m/x0) - r*t) - 1)

上面化简一下就是下面这个

解这个微分方程

求出sym来不要慌，s.x就能给结果了，注意这里是s=，要是ans就是ans.x

本题来源B站大仙一品堂

无约束问题的MATLAB解法

2. 无约束极值问题的符号解

求多元函数f(x,y)=x3?y3+3?x2+3?y3?9?xf(x,y)=x^3-y^3+3*x^2+3*y^3-9*xf(x,y)=x3?y3+3?x2+3?y3?9?x的极值

具体为啥这样算，都是高数上的知识

clear

syms x y

% x=-100:100;

% y=-100:100;

f=x.^3-y.^3+3*x.^2+3*y.^2-9*x;

% plot3(x,y,z)

df = jacobian(f);%各求一阶偏导数

d2f = jacobian(df);%号求Hessian阵

[xx,yy] =solve(df) %求驻点

xx = double( xx) ;yy = double(yy);%名转化成双精度浮点型数据,下面判断特征值的正负须是数值型数据

for i =1: length( xx)

a=subs(d2f,{x,y},{xx(i),yy(i)});

b=eig