matlab 解常微分方程自用2

锦鲤投食机

已于 2022-06-08 22:11:33 修改

阅读量67

点赞数

于 2022-06-08 21:11:23 首次发布

本文链接：https://blog.csdn.net/skdsky/article/details/125192550

版权

微分方程 MATLAB 几何绘图符号计算偏微分

关键词由CSDN通过智能技术生成

clc;clear;
syms L d g a;
V =  (-d*d/(2*tan(a)*sin(a))+(L-d/sin(a))^2*(cos(a)/2))*g;
solve(diff(V,a,1)==0,a) %求关于a的一阶偏微分等于0时的解,多变量时指定变量
%%
% 画出 过（0，1）且导数为横坐标的函数
for j=1:10;
    Deltax =j/10;
    A=[0,1];
    for i=1:1:10*j
        Deltay=A(2)*Deltax;
        B=[A(1)+Deltax, A(2)+Deltay];
        line([A(1),B(1)],[A(2),B(2)],'color', [j/10, 0,0]);%color GB为0 R随j变大
        hold on;
        A=B;
    end 
end
line([3 3],[0 20],'color',[0 0 1])%画出（3，0）和（3，30）的直线
axis([0 5 0 20])

%%
%解微分方程
%dy/dx=1+y^2
dsolve('Dy=1+y^2');%求一阶导
%二阶导为本身，且过原点，一阶导为0的函数
dsolve('D2Y=Y', 'DY(0)=0','Y(0)=2');