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

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

24小時(shí)熱門版塊排行榜

返回列表

gyq274706322

銀蟲 (正式寫手)

應(yīng)助: 6 (幼兒園)
金幣: 2727.2
散金: 200
帖子: 747
在線: 255.8小時(shí)
蟲號(hào): 1084211
注冊(cè): 2010-08-29
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)

[求助] 大家看看我的程序怎么了？

我的M函數(shù)文件是
function xdot=dianci1(t,x)
global w
m=19.6;g=9.8;la=0.116;lb=0.144;l=la+lb;
Jox=0.393;Joy=0.393;Joz=0.023;
ca=0.0004;cb=0.0004;
Ioxa=4;Ioya=4;Ioxb=4;Ioyb=4;
ra=0.0442;Da=0.08;Ara=pi*Da*ra;
rb=0.0442;Db=0.08;Arb=pi*Db*rb;
axy=0.16;
Nra=57;Nrb=57;miu0=4e-7*pi;
K=1;%K=7800;
KP=3.8;KI=200;KD=0.01;
m1=lb*m/l;m2=la*m/l;
e1=0;e2=0;

ixa=K*(KP*x(1)+KI*x(2)+KD*x(1));
iya=K*(KP*x(3)+KI*x(4)+KD*x(3));
ixb=K*(KP*x(5)+KI*x(6)+KD*x(5));
iyb=K*(KP*x(7)+KI*x(8)+KD*x(7));
%u(k)=kp*x(1)+kd*x(2)+ki*x(3); %PID Controller;
%x(1)=error(k);                      %Calculating P
%x(2)=(error(k)-error_1)/ts;          %Calculating D
%x(3)=x(3)+error(k)*ts;             %Calculating I
%xi(k)=x(3);
fra=miu0*Nra^2*Ara*(Ioxa+ixa)^2/(4*(ca+x(1))^2);
fla=miu0*Nra^2*Ara*(Ioxa-ixa)^2/(4*(ca-x(1))^2);
fta=miu0*Nra^2*Ara*(Ioya+iya)^2/(4*(ca+x(3))^2);
fba=miu0*Nra^2*Ara*(Ioya-iya)^2/(4*(ca-x(3))^2);
frb=miu0*Nrb^2*Arb*(Ioxb+ixb)^2/(4*(cb+x(5))^2);
flb=miu0*Nrb^2*Arb*(Ioxb-ixb)^2/(4*(cb-x(5))^2);
ftb=miu0*Nrb^2*Arb*(Ioyb+iyb)^2/(4*(cb+x(7))^2);
fbb=miu0*Nrb^2*Arb*(Ioyb-iyb)^2/(4*(cb-x(7))^2);

Fxa=fra-fla+axy*(x(1)/ca)*(fta+fba);
Fya=fta-fba+axy*(x(3)/ca)*(fra+fla);
Fxb=frb-flb+axy*(x(5)/cb)*(ftb+fbb);
Fyb=ftb-fbb+axy*(x(7)/cb)*(frb+flb);
fx=m1*e1*w^2*cos(w*t)+m2*e2*w^2*cos(w*t);
fy=-m1*e1*w^2*sin(w*t)+m2*e2*w^2*sin(w*t);
A=-Fxa-Fxb+fx;
B=-Fya-Fyb+m*g+fy;
C=Fxa*la-Fxb*lb-Joz/l*w*x(2)+Joz*w/l*x(4);
D=Fya*la-Fyb*lb+Joz/l*w*x(6)+Joz*w/l*x(8);

xdot=[x(2);
   (Joy*A*l-C*l*la*m)/(Joy*la*m+Joy*lb*m);
   x(4);
   (Joy*A*l+ C*l*lb*m)/(Joy*la*m+Joy*lb*m);
   x(6);
   (Jox*B*l-D*l*la*m)/(Jox*la*m+Jox*lb*m);
   x(8);
   (Jox*B*l+D*l*lb*m)/(Jox*la*m+Jox*lb*m)];

在matlab中運(yùn)行程序
clc;
clear;
%global w
w=100;
x0=[0,0,0,0,0,0,0,0];
for n=1:60
   [t,x]=ode45('dianci1',[0 10],x0);
x0=x(end,

;
hold on
plot(x(:,1),x(:,2));
end
但是結(jié)果怎么會(huì)出錯(cuò)
錯(cuò)誤是
??? Error using ==> odearguments at 116
DIANCI1 returns a vector of length 4, but the length of initial conditions vector is 8.
The vector returned by DIANCI1 and the initial conditions vector must have the same number
of elements.

Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

Error in ==> jiedianci at 7
[t,x]=ode45('dianci1',[0 10],x0);
我很糾結(jié)。我的方程是8個(gè)，為什么說(shuō)是返回4個(gè)？
這個(gè)我該怎么調(diào)式，才能出結(jié)果？
希望大家給點(diǎn)建議。
謝謝了。

回復(fù)此樓

» 猜你喜歡

一志愿中南化學(xué)337求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有23人回復(fù)
265求調(diào)劑已經(jīng)有9人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有10人回復(fù)
招收調(diào)劑碩士已經(jīng)有12人回復(fù)
296求調(diào)劑已經(jīng)有8人回復(fù)
招收博士1-2人已經(jīng)有3人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有11人回復(fù)
274求調(diào)劑已經(jīng)有8人回復(fù)
290求調(diào)劑已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

大家?guī)臀铱匆粋€(gè)fortran的程序，我總是計(jì)算不出正確的結(jié)果已經(jīng)有20人回復(fù)
【活動(dòng)】大家找工作更看重啥已經(jīng)有114人回復(fù)
大家看看這個(gè)投稿過(guò)程怪不怪？已經(jīng)有16人回復(fù)
請(qǐng)大家?guī)臀铱纯次业碾娯堅(jiān)趺椿厥?/a> 已經(jīng)有4人回復(fù)

【求助】急~~~~請(qǐng)大家?guī)兔纯催@臺(tái)液相是怎么了！�。� 已經(jīng)有9人回復(fù)
【求助】【求助】大家?guī)臀铱纯�，這是怎么回事？推薦信~~~~ 已經(jīng)有8人回復(fù)
【求助】大家看看，matlab安裝過(guò)程的截圖，找出錯(cuò) 已經(jīng)有14人回復(fù)
【求助/交流】大家看看我的蛋白膠是怎么了？已經(jīng)有29人回復(fù)
投了篇CPL，請(qǐng)大家看看是怎么回事？已經(jīng)有5人回復(fù)

1樓 2011-09-08 10:15:15

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

gyq274706322

銀蟲 (正式寫手)

應(yīng)助: 6 (幼兒園)
金幣: 2727.2
散金: 200
帖子: 747
在線: 255.8小時(shí)
蟲號(hào): 1084211
注冊(cè): 2010-08-29
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)

運(yùn)行
clc;
clear;
%global w
w=100;
x0=[0,0,0,0,0,0,0,0];
for n=1:60
[t,x]=ode45('dianci1',[0 10],x0);
x0=x(end,

;
hold on
plot(x(:,1),x(:,2));
end

贊一下

回復(fù)此樓

2樓2011-09-08 10:17:02

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

gyq274706322

銀蟲 (正式寫手)

應(yīng)助: 6 (幼兒園)
金幣: 2727.2
散金: 200
帖子: 747
在線: 255.8小時(shí)
蟲號(hào): 1084211
注冊(cè): 2010-08-29
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)

怎么會(huì)出現(xiàn)那個(gè)小圖��？

贊一下

回復(fù)此樓

3樓2011-09-08 10:17:39

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

信彼南山

木蟲 (著名寫手)

應(yīng)助: 33 (小學(xué)生)
金幣: 4142.9
散金: 1221
紅花: 16
帖子: 1178
在線: 233.5小時(shí)
蟲號(hào): 1133529
注冊(cè): 2010-10-27
專業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)

★
xiegangmai(金幣+1): 謝謝提示！ 2011-09-11 15:53:50

可以用
[ code] [/ code] （去掉空格）
框住你的代碼

贊一下(1人)

回復(fù)此樓

4樓2011-09-11 11:51:08

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

信彼南山

木蟲 (著名寫手)

應(yīng)助: 33 (小學(xué)生)
金幣: 4142.9
散金: 1221
紅花: 16
帖子: 1178
在線: 233.5小時(shí)
蟲號(hào): 1133529
注冊(cè): 2010-10-27
專業(yè): 導(dǎo)航、制導(dǎo)與傳感技術(shù)

肯定是方程有錯(cuò)了。

贊一下

回復(fù)此樓

5樓2011-09-11 17:39:14

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

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

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 085700資源與環(huán)境308求調(diào)劑 +10	墨墨漠 2026-03-18	11/550	2026-03-20 11:57 by 學(xué)員8dgXkO
[考研] 317求調(diào)劑 +4	申子申申 2026-03-19	9/450	2026-03-20 11:08 by 申子申申
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +20	rare12345 2026-03-18	20/1000	2026-03-20 08:42 by 無(wú)際的草原
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 324分 085600材料化工求調(diào)劑 +3	llllkkkhh 2026-03-18	3/150	2026-03-19 14:22 by houyaoxu
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-19	3/150	2026-03-19 13:36 by houyaoxu
[考研] 0817調(diào)劑 +3	沒有答案_ 2026-03-14	3/150	2026-03-19 09:51 by Xu de nuo
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

gyq274706322

[求助] 大家看看我的程序怎么了？

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

gyq274706322

gyq274706322

信彼南山

信彼南山

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 大家看看我的程序怎么了？

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: