神马作文网如何用matlab画图求解
来源：学生作业帮 编辑：神马作文网作业帮 分类：数学作业 时间：2024/11/18 17:55:32
如何用matlab画图求解
能用matlab对这个函数进行画图：
边长为1m的方桌四角上各有一条小虫,每只虫同时以同样的速度v=0.1（cm/s）按逆时爬向它相邻的虫子,求各个虫子的爬行轨迹、爬行路程长及汇合时间?
分析：将左下角虫轨迹设成f（x）=y,然后设虫的位置为（0,0）（1,0）（1,1）（0,1）,然后经时间t后,（0,0）变为（x,y）,（1,0）变为（1-y,x）,左下角虫的速度方向应指向右下角的虫,则有dy/dx=（x-y）/（1-y-x）
因此我想用matlab对函数dy/dx=（x-y）/（1-y-x）进行模拟,请问该怎样做,希望各位高手写出代码.
能用matlab对这个函数进行画图：
边长为1m的方桌四角上各有一条小虫,每只虫同时以同样的速度v=0.1（cm/s）按逆时爬向它相邻的虫子,求各个虫子的爬行轨迹、爬行路程长及汇合时间?
分析：将左下角虫轨迹设成f（x）=y,然后设虫的位置为（0,0）（1,0）（1,1）（0,1）,然后经时间t后,（0,0）变为（x,y）,（1,0）变为（1-y,x）,左下角虫的速度方向应指向右下角的虫,则有dy/dx=（x-y）/（1-y-x）
因此我想用matlab对函数dy/dx=（x-y）/（1-y-x）进行模拟,请问该怎样做,希望各位高手写出代码.
%我想用matlab对函数dy/dx=（x-y）/（1-y-x）进行模拟,请问该怎样做,
%解微分方程!
%归一化：
%令y=y（1）
%x=t=y(2)
%dy(1)/dt=(y(2)-y(1))/(1-y(1)-y(2))
%dy(2)/dt=1;
%函数文件
founction dy=fun1(t,y)
dy=zeros(2,1);
dy=[(y(2)-y(1))/(1-y(1)-y(2));1];
%以上保存为fun1.m文件
%以下是脚本程序
clear
ts=0:0.01:%时间范围
y0=[?];%函数边界条件y(0),x(0)
[t,y]=ode45('fun1',ts,y0);
plot(y(:,1),y(:,2))
%解微分方程!
%归一化：
%令y=y（1）
%x=t=y(2)
%dy(1)/dt=(y(2)-y(1))/(1-y(1)-y(2))
%dy(2)/dt=1;
%函数文件
founction dy=fun1(t,y)
dy=zeros(2,1);
dy=[(y(2)-y(1))/(1-y(1)-y(2));1];
%以上保存为fun1.m文件
%以下是脚本程序
clear
ts=0:0.01:%时间范围
y0=[?];%函数边界条件y(0),x(0)
[t,y]=ode45('fun1',ts,y0);
plot(y(:,1),y(:,2))