如何用matlab画图求解
==replied by huxuec ，matlabsky转载仅为帮助初学者==%解微分方程！！ %归一化： %令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))
MATLAB绘图问题求解
这是三次样条差值函数，系数前面已经算出来了，现在要做出函数图像，x1就是函数表达式里面的未知数，symsx1是参考别人的，不知道有没有必要
matlab绘图求解。
假设动能和重力势能相互转换，则有
1/2 m1glsinO+ m2g(l+r)sinO = 1/2 (m1+m2)v^2
从而，w=v/(l+r), a=v^2/(l+r) (公式太长不写了，具体见m脚本)
m脚本如下，
g=9.8;
l=1;
r=0.2;
m1=2; m2=8;
phi = input('Please input angle phi:');
v = sqrt( (m1*g*l*sin(phi)+2*m2*g*(l+r)*sin(phi))/(m1+m2) );
w = v/(l+r);
a = (v^2)/(l+r);
subplot(2,1,1);
plot(phi,w,'r-.*');
text(phi+.1, w+.1, sprintf('(%.3f, %.3f)', phi, w), 'FontWeight', 'Bold');
title('角度-角速度');
hold on;
subplot(2,1,2);
plot(phi,a,'r-.*');
text(phi+.1, a+.1, sprintf('(%.3f, %.3f)', phi, a), 'FontWeight', 'Bold');
title('角度-加速度');
输入phi= 1.234 可得
MATLAB画图程序求解
就我所知，polar命令存在以下问题：
1、一次调用只能绘制一条曲线；
2、如果使用hold on命令多次绘图，polar不能自动调整坐标范围，所以，绘图时应考虑首先绘制坐标范围最大的那条曲线。
参考代码：
n=10;
x=0:0.01*pi:2*pi;
for i=n:-1:1
y=50*(1+sqrt(0.64*i^2+0.36))-32./(1-0.6*cos(x));
polar(x,y)
hold on
end