اشتباها نمودار سرعت رو رسم میکردم
حل شد

m=2;
k=500;
x0=0.1;
v0=5;
wn=sqrt(k/m);
y0=[x0;v0];
tspan=[0 1];
% y = zeros(2,2);
springmass=@(t,y)[y(2);-k/m*y(1)];
[t,y] = ode45(springmass,tspan,y0);
x = x0*cos(wn*t)+v0/wn*sin(wn*t);
plot(t,x)
figure
plot(t,y(:,2))

آیا کد زیر صحیح میباشد؟
clear all
clc
M = 300;m=50;%sabethaye jerm
k1 = 3000;k2=30000;%sabethaye fanar
D=600;%sabete damper
func = @(t,x) [x(3);x(4);-k1/M*x(1)-k1/M*x(2)-D/M*x(3)+D/M*x(4);k1/m*x(1)-(k1+k2)/m*x(2)-D/m*x(3)-D/m*x(4)]; %systeme talighe faale khodro
x0 = [0.2;0;0;0];%sharayete avali
tspan = [0 5];%bazeye zamanie shabih sazi
[t, x] = ode23(func,tspan,x0);%hale adadie moadelat diffranciel
plot(t,x(:,1));%rasme mogheiat bar hasbe zaman
title(‘Position vs. Time’);%name mehvar
hold
plot(t,x(:,2));%rasme sorat bar hasbe zaman
title(‘Velocity vs. Time’);%name mehvar

مثلاً در سیستم زیر چجوری باید کدنویسی کنیم؟
من هر کاری میکنم جواب نمیگیرم
x1’=x3
x2’=x4
x3’= -(k1/M) x1 – (k1/M) x2 – (D/M) x3 + (D/M) x4
x4’= (k1/m) x1 – (k1+k2) /m x2 – (D/m) x3 – (D/m) x4

لطفا توضیح بیشتری بدهید.

خوب ، متغییرهای حالت 4 تا هستند.

مشکل چیست ؟