matlab计算阻尼运动并画图

你先学下matlab基本知识吧

代码供参考：n=1.5;%x1=φ,x2=Rdx=@(z,x)[1/(x(2)*(n+(n+1)*(2*cos(x(1))+cos(x(1))*cos(x(1)))));-1/tan(x(1))];[

x=(0.1:0.1:80);y=zeros(1,800); fori=1:1:800   ifi<=100    

有求微分方程的解,还要有初值才行再问：补充了抱歉再答：解析解解法s=dsolve('D2x-0.05*Dy+0.15*x=0','D2y-0.05*Dx+0.15*y=0

clc;clearlamta=0.5*10^-4;r=0.3;s=pi*r^2;j=1;fork=0:10:40total=0;fori=0:k-1total=total+((lamta*s)^i/f

是阻尼运动,最终会停下来在中间某一位置两侧的来回运动叫往复运动.如果当中那一点是平衡位置,则这一往复运动叫做振动.振动一定是往复运动.这种运动

由于运动物体收到阻力而使运动位移的峰值逐渐减小的现象!

你需要把上面的部分保存为Lagrangenew.m文件.文件头应加 functions=Lagrangenew(x,y,x0)后面是nx=length(x);ny=length(y);.然后

symstv=dsolve('Dv=(190.708-90.64*v^2)/47.27','v(0)=0','t');t=0:0.00001:0.002;v=eval(v);plot(t,v)使用这样

>> y=dsolve('Dy=1-y^2','y(0)=0') y = tanh(t)>> 

x=0:0.1:100;y1=sin(x)./x;figure,plot(x,y1);y2=(1+x).^(1./x);figure,plot(x,y2);再问：我决定采纳你的答案你能在帮我接到题吗？

试试~发站内信

N=1000;>>fori=1:Na=randint(1,i);sum1(i)=sum(a);P(i)=sum1(i)/i;end>>n=1:N;>>plot(n,P);xlabel('次数n');y

由于方程比较复杂,解析解不能用初等函数表示只是要获得图像的话,用数值计算的办法可能更方便fun=@(t,y)1.44*(10^9)*(1-y).*exp(-109170./(8.314*t));[T,

知道x,y坐标,画图只需要plot(x,y)即可,但是要预测就需要看你数据的规律了,比如说最简单的要是直线,你可以把直线的方程求出来,就可以预测了,其他的曲线一般都可以用多项式等拟合,在预测,

这是一个一阶微分方程组的简化表达方式.设i的取值范围是1—N,那么j的取值范围是1-N；编程思路如下：alpha=zeros(N,N);%应该alpha已知%x=zeros(1,N);%x是所求向量,

参考代码：a=0.5; b=0.6;dIdt=@(t,I)a*(1-I)*I-b*I;ode45(dIdt,[0 10],0)存在的问题是,如果初值为0,可以计算出导数dI/dt=

clc;clear;n1=-10:1:-1;y1=0*n1;n2=0:1:50;y2=(0.8).^n2;n=[n1,n2]; y=[y1,y2];stem(n,y);

无阻尼运动的体征之一是受外界周期驱动力用于补充不断消耗的能量,使振幅保持不变.而简谐运动是物体所受合外力和偏离平衡位置的位移成正比.明显你可以加一个驱动力使得加上它以后物体所受合外力不满足和位移成正比

sol1=dsolve('D2x=-Dx-1','x(0)=1,Dx(0)=1','a');>>sol2=dsolve('D2x=-2*Dx-1','x(0)=1,Dx(0)=1','a');>>so