s20090281

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 460.9
散金: 177
紅花: 1
帖子: 235
在線: 31.2小時(shí)
蟲號: 1185039
注冊: 2011-01-07
性別: MM
專業(yè): 人工智能與知識工程

[求助] matlab中M文件的調(diào)用

我最近參考別人的程序?qū)懥艘粋€(gè)自己的程序但是代入不了fitness函數(shù)，請大家指點(diǎn)一下：我要做的是y(k)=a1*y(k-1)+b1*u(k-1)+a2*y(k-2)+a3*y(k-3)函數(shù)求a1,a2,a3,b1.已知y和u，我把適應(yīng)度取為399組數(shù)據(jù)的[y(k)-a1*y(k-1)-b1*u(k-1)-a2*y(k-2)-a3*y(k-3)]^2最小

CODE:

function [out]=fitness(x)

CODE:

ai=x(i,：)           %第幾組數(shù)據(jù)就是取第幾個(gè)

CODE:

a1=ai(1)

CODE:

a2=ai(2)

CODE:

a3=ai(3)

CODE:

a4=ai(4)

CODE:

u=load('F:\測試記錄控制.txt');u1=u(2:399,2);

CODE:

y=load('F:\測試記錄反饋.txt');y1=y(2:400,2);

CODE:

 c=0;

CODE:

for k=3:399

CODE:

f=[y(k)-a1*y1(k-1)-a2*y1(k-2)-a3*y1(k-3)-b1*u1(k-1)]^2        %求出各級數(shù)據(jù)的平方和并逐漸相加得到一組對對應(yīng)變量的適應(yīng)度

CODE:

c=c+f

CODE:

Msum(k)=c

CODE:

end

CODE:

out=Msum(:,k);                %輸出適應(yīng)度

那么在粒子群算法中應(yīng)如何調(diào)用這個(gè)函數(shù)呢？當(dāng)位置和速度更新后，怎么計(jì)算這時(shí)的適應(yīng)度呢？

CODE:

ffitness(i)=fitness(pop(i,:))

這樣調(diào)用為什么會出錯(cuò)呢？請各位幫我這個(gè)新手解答一下，謝謝

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

lijing870728

» 猜你喜歡

306求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級已過已經(jīng)有9人回復(fù)
289材料與化工（085600）B區(qū)求調(diào)劑已經(jīng)有4人回復(fù)
一志愿武理材料工程348求調(diào)劑已經(jīng)有6人回復(fù)
08工科 320總分求調(diào)劑已經(jīng)有11人回復(fù)
一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑已經(jīng)有6人回復(fù)
323求調(diào)劑已經(jīng)有5人回復(fù)
284求調(diào)劑已經(jīng)有5人回復(fù)
317求調(diào)劑已經(jīng)有16人回復(fù)
287求調(diào)劑已經(jīng)有9人回復(fù)

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

m文件調(diào)用問題已經(jīng)有4人回復(fù)
Matlab中如何調(diào)用M文件作圖已經(jīng)有9人回復(fù)
求助大俠幫忙，有關(guān)matlab模擬數(shù)據(jù) 已經(jīng)有29人回復(fù)
怎樣將txt文件讀入到matlab中，謝謝啦已經(jīng)有4人回復(fù)
歸一化問題已經(jīng)有11人回復(fù)
請教Matlab高手已經(jīng)有5人回復(fù)
【討論】關(guān)于一個(gè)m文件中函數(shù)調(diào)用及畫圖問題已經(jīng)有7人回復(fù)
【求助】matlab 中，幾個(gè)m文件調(diào)用已經(jīng)有14人回復(fù)
【求助】matlab生成旋轉(zhuǎn)體？已經(jīng)有9人回復(fù)
【求助】如何獲得matlab中當(dāng)前fig圖像的信息已經(jīng)有6人回復(fù)
【求助】matlab中如何將程序的結(jié)果輸出到一個(gè)文件中已經(jīng)有11人回復(fù)
【求助】VC++.NET/VC++ 調(diào)用EXCEL文件顯示數(shù)據(jù)再控件上【已解決】已經(jīng)有13人回復(fù)
MATLAB中如何導(dǎo)入mat文件?如何查看此mat文件已經(jīng)有4人回復(fù)

1樓 2011-09-20 15:40:42

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

信彼南山

木蟲 (著名寫手)

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

★
xiegangmai(金幣+1): 謝謝參與！ 2011-09-20 22:33:25

你這代碼寫的，放在一個(gè)code里面就行了

贊一下(1人)

回復(fù)此樓

2樓2011-09-20 16:08:20

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

s20090281

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 460.9
散金: 177
紅花: 1
帖子: 235
在線: 31.2小時(shí)
蟲號: 1185039
注冊: 2011-01-07
性別: MM
專業(yè): 人工智能與知識工程

引用回帖:

2樓: Originally posted by 信彼南山 at 2011-09-20 16:08:20:
你這代碼寫的，放在一個(gè)code里面就行了

你可否幫我解決一下這個(gè)問題，我的只能運(yùn)行一組當(dāng)運(yùn)行第二組時(shí)就會出現(xiàn)。我在正文中這么調(diào)用函數(shù)
for i=1:size
pop(i,

=2*rand(size,1)-1
v(i,

=rand(size,1)
fitness(i)=fitness(pop(i,

)
end
??? Subscript indices must either be real positive integers or logicals.
Error in ==> myPSO at 14
fitness(i)=fitness(pop(i,

)
我查了一下是下標(biāo)不正確，可是那怎么改呢？

贊一下

回復(fù)此樓

3樓2011-09-20 16:20:21

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

信彼南山

木蟲 (著名寫手)

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

你這代碼也太亂了吧
你不會重新整理一下啊

贊一下

回復(fù)此樓

4樓2011-09-20 19:08:37

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

s20090281

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 460.9
散金: 177
紅花: 1
帖子: 235
在線: 31.2小時(shí)
蟲號: 1185039
注冊: 2011-01-07
性別: MM
專業(yè): 人工智能與知識工程

引用回帖:

4樓: Originally posted by 信彼南山 at 2011-09-20 19:08:37:
你這代碼也太亂了吧
你不會重新整理一下啊

CODE:

function [out]=fitness(x)

          ai=x(i,：)           %第幾組數(shù)據(jù)就是取第幾個(gè)

          a1=ai(1)

          a2=ai(2)

          a3=ai(3)

          a4=ai(4)

          u=load('F:\測試記錄控制.txt');u1=u(2:399,2);) ;

          y=load('F:\測試記錄反饋.txt');y1=y(2:400,2);

          c=0;

         for k=3:399

              f=[y(k)-a1*y1(k-1)-a2*y1(k-2)-a3*y1(k-3)-b1*u1(k-1)]^2 ;      

 %求出各級數(shù)據(jù)的平方和并逐漸相加得到一組對對應(yīng)變量的適應(yīng)度       

             c=c+f;

            Msum(k)=c;

       end

       out=Msum(:,k);                %輸出適應(yīng)度

贊一下

回復(fù)此樓

5樓2011-09-20 20:04:17

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

s20090281

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 460.9
散金: 177
紅花: 1
帖子: 235
在線: 31.2小時(shí)
蟲號: 1185039
注冊: 2011-01-07
性別: MM
專業(yè): 人工智能與知識工程

引用回帖:

4樓: Originally posted by 信彼南山 at 2011-09-20 19:08:37:
你這代碼也太亂了吧
你不會重新整理一下啊

這樣可以了嗎？可以幫幫我了吧。

謝謝

贊一下

回復(fù)此樓

6樓2011-09-20 20:05:19

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

信彼南山

木蟲 (著名寫手)

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

★
xiegangmai(金幣+1): 謝謝參與！ 2011-09-20 22:33:42

ai=x(i,：) %第幾組數(shù)據(jù)就是取第幾個(gè)
你這第一句就有問題，i沒有定義

贊一下(1人)

回復(fù)此樓

7樓2011-09-20 20:45:20

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

信彼南山

木蟲 (著名寫手)

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

★ ★
xiegangmai(金幣+2): 辛苦了！ 2011-09-20 22:33:57

你試試這樣行不行？

CODE:

function [out]=fitness(a1,a2,a3,a4)

%          ai=x(i,：)           %第幾組數(shù)據(jù)就是取第幾個(gè)

%          a1=ai(1)

%          a2=ai(2)

%          a3=ai(3)

%          a4=ai(4)

          u=load('F:\測試記錄控制.txt');u1=u(2:399,2);) ;

          y=load('F:\測試記錄反饋.txt');y1=y(2:400,2);

          c=0;

         for k=3:399

              f=[y(k)-a1*y1(k-1)-a2*y1(k-2)-a3*y1(k-3)-b1*u1(k-1)]^2 ;      

%求出各級數(shù)據(jù)的平方和并逐漸相加得到一組對對應(yīng)變量的適應(yīng)度       

             c=c+f;

            Msum(k)=c;

       end

       out=Msum(:,k);                %輸出適應(yīng)度

贊一下(1人)

回復(fù)此樓

8樓2011-09-20 20:49:40

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

信彼南山

木蟲 (著名寫手)

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

★ ★
xiegangmai(金幣+2): 謝謝參與！ 2011-09-20 22:34:07

你這個(gè)編程的習(xí)慣簡直太不好了
在for循環(huán)里面怎么能放讀寫文件的指令呢？而且每次讀的還都是同一個(gè)文件。
完全可以放在主函數(shù)里面，讀一次就行了么

贊一下(1人)

回復(fù)此樓

9樓2011-09-20 20:51:45

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

s20090281

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 460.9
散金: 177
紅花: 1
帖子: 235
在線: 31.2小時(shí)
蟲號: 1185039
注冊: 2011-01-07
性別: MM
專業(yè): 人工智能與知識工程

引用回帖:

8樓: Originally posted by 信彼南山 at 2011-09-20 20:49:40:
你試試這樣行不行？
[code]
function [out]=fitness(a1,a2,a3,a4)
%       ai=x(i,：)          %第幾組數(shù)據(jù)就是取第幾個(gè)
%       a1=ai(1)
%       a2=ai(2)
%       a3=ai(3)
%       a ...

不好意思這個(gè)不行，對于你說的問題，我還不太懂，以前從來沒寫過程序，都是新寫的，呵呵，不過真的很感謝你，抽空幫我看

贊一下

回復(fù)此樓

10樓2011-09-20 22:37:27

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

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

[考研] 269專碩求調(diào)劑 +6	金恩貝 2026-03-21	6/300	2026-03-22 14:31 by ColorlessPI
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +4	zxts12 2026-03-21	7/350	2026-03-22 09:57 by zxts12
[考研] 0856材料專碩353求調(diào)劑 +4	NIFFFfff 2026-03-20	4/200	2026-03-22 09:49 by 2026paper
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4	七月-七七 2026-03-21	4/200	2026-03-21 18:20 by 學(xué)員8dgXkO
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 307求調(diào)劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[考研] 266求調(diào)劑 +3	哇呼哼呼哼 2026-03-20	3/150	2026-03-21 16:46 by barlinike
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 材料專碩英一數(shù)二306 +7	z1z2z3879 2026-03-18	7/350	2026-03-20 23:48 by JourneyLucky
[考研] 304求調(diào)劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考博] 26博士申請 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷