版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

xuanshan

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 51.1
帖子: 11
在線: 3.3小時
蟲號: 3792579
注冊: 2015-04-07
專業(yè): 機(jī)械摩擦學(xué)與表面技術(shù)

[求助] 程序問題已有1人參與

function dy=rigid(t,y)
a=0.06;%輪胎的接地印跡長度
w=63;%輪胎的轉(zhuǎn)速
R=0.317;%輪胎半徑
u=1;%外界風(fēng)速
ls=1.293;%空氣密度
lg=2.5*10^3;%磨粒密度
d=10^(-5);%磨粒直徑
g=9.806;%重力加速度
m=1.79*10^(-6);%空氣動力粘度
k1=3*ls/(4*lg*d);
g1=(1-ls/lg)*g;
dy=zeros(4,1);%一個列向量
V=sqrt(y(2)^2+y(4)^2);
U=w*sqrt(R^2-(y(1)^2+y(3)^2));
Re1=ls*d*abs(u-V)/m;
Re2=ls*d*abs(U-V)/m;
Cd1=24*(1+0.0975*Re1-0.636*Re1^2)/Re1;
Cd2=24*(1+0.0975*Re2-0.636*Re2^2)/Re2;
k2=3*Cd2*ls/(4*lg*d);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=k1*(Cd1*(u-y(2))^2+Cd2*(w*(y(1)-a)/(sqrt(R^2-((y(1)-a)^2+(y(3)+sqrt(R^2-a^2))^2))+y(2))^2));
dy(4)=-g1+k2*(w*(y(3)+sqrt(R^2-a^2))/sqrt(R^2-((y(1)-a)^2+(y(3)+sqrt(R^2-a^2))^2)+y(4))^2);
dy=[dy(1) dy(2) dy(3) dy(4)]';
end
這是我的程序代碼，但是在運(yùn)算的時候出現(xiàn)下列情況
tspan=linspace(0,1.5,1000);
[T Y]=ode23s(@rigid,tspan,[0 19.6 0 3.87]);
Warning: Failure at t=7.361808e-06. Unable to meet integration tolerances without reducing the step size below the smallest value
allowed (2.615440e-20) at time t.
> In ode23s at 402
麻煩哪位大神給指導(dǎo)一下，謝謝！

回復(fù)此樓

» 猜你喜歡

304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)
085410人工智能專碩317求調(diào)劑（0854都可以）已經(jīng)有3人回復(fù)
330求調(diào)劑已經(jīng)有3人回復(fù)

1樓 2015-10-13 10:17:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

512095518

專家顧問 (正式寫手)

專家經(jīng)驗(yàn): +52
仿真EPI: 1
應(yīng)助: 118 (高中生)
金幣: 7665.4
紅花: 78
帖子: 712
在線: 186.1小時
蟲號: 3195576
注冊: 2014-05-10
性別: GG
專業(yè): 大氣環(huán)境與全球氣候變化
管轄: 計算模擬

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

應(yīng)該是你的常微分方程組四個解寫得有問題吧，估計是太過于復(fù)雜才導(dǎo)致ode23求解不了。直接用數(shù)值會好一些�；蚴悄阍囅耾de45看看。具體你還需要自己去查下類似的方程該怎么寫。

下面是我用過的供你參考，其中a-n都是參數(shù)，都是提前算好的。
function dx = hanshu(t, x)
global a b c d e f g h i j k m l n
      dx = zeros(4,1);
      dx(1) = a + b*x(3) + c*x(2) - d*x(1);
      dx(2) = e + f*x(1) + g*x(3) + h*x(4) - i*x(2);
      dx(3) = j + k*x(1) - l*x(3);
      dx(4) = m*x(2) - n*x(4);
end

贊一下

回復(fù)此樓

沒有極限，只有超越！

2樓2015-10-14 15:31:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 xuanshan 的主題更新