我想用ode45函數(shù)做一個(gè)定加速度的物體三維運(yùn)動(dòng)的模型，自己寫了框架，在最后畫圖時(shí)，圖形不出來，大

做簡諧運(yùn)動(dòng)的物體速度加速度都符合正弦定理。

創(chuàng)新互聯(lián)公司2013年成立，公司以成都做網(wǎng)站、成都網(wǎng)站建設(shè)、系統(tǒng)開發(fā)、網(wǎng)絡(luò)推廣、文化傳媒、企業(yè)宣傳、平面廣告設(shè)計(jì)等為主要業(yè)務(wù)，適用行業(yè)近百種。服務(wù)企業(yè)客戶上千多家，涉及國內(nèi)多個(gè)省份客戶。擁有多年網(wǎng)站建設(shè)開發(fā)經(jīng)驗(yàn)。為企業(yè)提供專業(yè)的網(wǎng)站建設(shè)、創(chuàng)意設(shè)計(jì)、宣傳推廣等服務(wù)。 通過專業(yè)的設(shè)計(jì)、獨(dú)特的風(fēng)格，為不同客戶提供各種風(fēng)格的特色服務(wù)。

1、如果質(zhì)點(diǎn)的位移與時(shí)間的關(guān)系遵從正弦函數(shù)的規(guī)律，即它的振動(dòng)圖像是一條正弦曲線，這樣的振動(dòng)叫做簡諧運(yùn)動(dòng)。簡諧運(yùn)動(dòng)是最基本也最簡單的機(jī)械振動(dòng)。當(dāng)某物體進(jìn)行簡諧運(yùn)動(dòng)時(shí)物體所受的力跟位移成正比并且總是指向平衡位置。它是一種由自身系統(tǒng)性質(zhì)決定的周期性運(yùn)動(dòng)。實(shí)際上簡諧振動(dòng)就是正弦振動(dòng)。故此在無線電學(xué)中簡諧信號實(shí)際上就是正弦信號。

2、簡諧運(yùn)動(dòng)加速度就是回復(fù)力產(chǎn)生的加速度，回復(fù)力就是沿著物體運(yùn)動(dòng)方向的力，彈簧振子的回復(fù)力就是彈簧形變的力，單擺要做受力分析，回復(fù)力沿與擺線垂直方向。F=ma ，a就是加速度。

3、正弦定理是三角學(xué)中的一個(gè)基本定理，它指出在任意一個(gè)平面三角形中，各邊和它所對角的正弦值的比相等且等于外接圓半徑的2倍，即a/sinA = b/sinB =c/sinC = 2R（R為外接圓半徑）。

ode45的語法

[T,Y] = ode45(odefun,tspan,y0)

[T,Y] = ode45(odefun,tspan,y0,options)

[T,Y,TE,YE,IE] = ode45(odefun,tspan,y0,options)

sol = ode45(odefun,[t0tf],y0...)

[T,Y] = ode45(odefun,tspan,y0)

odefun 是函數(shù)句柄,可以是函數(shù)文件名,匿名函數(shù)句柄或內(nèi)聯(lián)函數(shù)名

tspan 是區(qū)間 [t0 tf] 或者一系列散點(diǎn)[t0,t1,...,tf]

y0 是初始值向量

T 返回列向量的時(shí)間點(diǎn)

Y 返回對應(yīng)T的求解列向量

[T,Y] = ode45(odefun,tspan,y0,options)

options 是求解參數(shù)設(shè)置,可以用odeset在計(jì)算前設(shè)定誤差,輸出參數(shù),事件等

[T,Y,TE,YE,IE] =ode45(odefun,tspan,y0,options)

在設(shè)置了事件參數(shù)后的對應(yīng)輸出

TE 事件發(fā)生時(shí)間

YE 事件發(fā)生時(shí)之答案

IE 事件函數(shù)消失時(shí)之指針i

sol =ode45(odefun,[t0 tf],y0...)

sol 結(jié)構(gòu)體輸出結(jié)果

用ode45求解微分方程

clear all

clc

k=5e4*0.5;

r0=0.0125;

Density=2650;

m=3.1416*r0^3*Density*4/3;

g=9.81;

cc=2*sqrt(k*m);

b=0.5;

c=cc*b;

h0=1;

v0=sqrt(2*(h0-r0)*g);

dydt=@(t,y)([y(2);m*g+k*(r0-y(1))-c*y(2)]);

t_end=.9;

x0=[g;v0];

[t,y]=ode45(dydt,[0,t_end],x0);

d2y=m*g+k*(r0-y(:,1))-c*y(:,2);

subplot(1,3,1),plot(t,y(:,1)),xlabel('t'),legend('y');

subplot(1,3,2),plot(t,y(:,2)),xlabel('t'),legend('dy/dt');

subplot(1,3,3),plot(t,d2y),xlabel('t'),legend('d2y/dt');