| 5 | 1/1 | 返回列表 |
| 查看: 2161 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
一生快樂新蟲 (初入文壇)
|
[求助]
matlab程序運(yùn)行總是出現(xiàn)下面的錯誤,請各位高手指點(diǎn)一下,謝謝了!急 急急
|
|
|
NIND=40; %個體數(shù)目 MAXGEN=500; %最大遺傳代數(shù) NVAR=7; %變量個數(shù) GGAP=0.8; %代溝 %創(chuàng)建區(qū)域描述器 FielDR=[0.01,0.01,0.01,0.01,0.01,0.01,0.01;0.9,0.9,0.9,0.9,0.9,0.9,0.9] Chrom=crtrp(NIND,FielDR)%創(chuàng)建初始種群 A=Chrom(i,j); for j=1:40; sum(j)=0; for m=1:7; sum(j)=A(j,m)+sum(j); end; end for j=1:40 for i=1:7 Chrom(i,j)=A(i,j)/sun(j); end end; gen=0; %代計(jì)數(shù)器 ObjV=objfun_1(Chrom) %計(jì)算初始種群個體的目標(biāo)函數(shù)值 while gen FitnV=ranking(ObjV) %分配適應(yīng)度值 SelCh=select('sus',Chrom,FitnV,GGAP) %選擇(隨即遍歷抽樣) SelCh=recombin('xovsp',SelCh,0.7) %重組(單點(diǎn)交叉) SelCh=mutbga(SelCh,FielDR) %變異 ObjVSel=objfun_1(SelCh) %子代函數(shù)值 [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel) %重插入 gen=gen+1; [Y,I]=min(ObjV); %最優(yōu)解,種群序號 Y w=(Chrom(I, ) %輸出自變量值end w 該段程序 運(yùn)行錯誤顯示 Undefined function or method 'crtrp' for input arguments of type 'double'. function f=objfun_1(w) R51=[7.5 0.15 3.0 0.002 1.0 0.010 0.1 12.80 0.13 0.65 2.73 3.60 0.07 2.92 11.20 0.27 0.5 0.68 5.70 0.09 0.98 13.20 0.44 0.3 2.18 5.30 0.08 2.78 9.90 0.20 0.5 1.09 5.40 0.06 3.04 11.70 0.31 0.35 0.55 5.70 0.15 3.57 9.00 0.27 0.35 1.37 7.30 0.09 1.41 10.70 0.43 0.35 3.82 7.40 0.14 2.70 8.70 0.27 0.3 2.05 6.20 0.10 2.69 11.60 0.47 0.2 1.91 9.50 0.22 2.21 9.70 0.28 0.35 4.50 9.80 0.17 1.81 9.90 0.39 0.3 3.69 7.60 0.18 2.23 6.80 0.60 0.3 4.91 9.60 0.27 1.36 6.00 0.35 0.65 7.10 11.80 0.16 1.40 6.70 0.63 0.5 4.37 9.80 0.30 1.62 14.00 1.55 0.5 9.56 8.30 0.22 2.99 9.00 0.73 0.6 1.09 9.60 0.21 2.67 10.90 0.50 0.65 12.01 9.20 0.28 2.02 10.90 0.55 0.35 5.19 6.00 0.16 1.72 8.10 0.78 0.5 8.19 8.90 0.22 2.12 14.90 0.55 0.5 8.46 10.80 0.19 2.27 11.30 0.24 0.35 5.46 6.60 0.15 1.01 12.50 0.90 0.5 1.75 5.70 0.10 3.19 9.60 0.49 0.35 3.28 8.00 0.14 1.85 11.20 0.68 0.65 7.50 8.20 0.18 2.27]; A1=R51 sizes=size(A1)%統(tǒng)計(jì)行列數(shù) row=sizes(1); col=sizes(2); [aa,b]=size(w); %決策矩陣規(guī)范化 sum=0; for j=1:7; for m=1:25; a=A1(m,j)*A1(m,j); sum=a+sum; end; b=sqrt(sum); for i=1:25 R1(i,j)=A1(i,j)/b; end; end; %運(yùn)用目標(biāo)函數(shù)法確定權(quán)重系數(shù) sumx=0; x1(1:col)=0; for mm=1:aa; sum=0; for j=1:col for i=2:row xa(j)=(R1(i,j)-R1(1,j))^2; %每列中的每個元素與第一個元素的差取平方,即目標(biāo)函數(shù) x1(j)=x1(j)+xa(j);% end sumx(j)=x1(j)*w(m,j)*w(m,j);%目標(biāo)函數(shù)之和 sum=sum+sumx(j); end f(mm,1)=sum end 這段程序運(yùn)行錯誤顯示??? Input argument "w" is undefined. |

至尊木蟲 (著名寫手)
驃騎將軍

新蟲 (初入文壇)
|
我的這個程序就是調(diào)用的objfun_1(w)那個程序 現(xiàn)在的問題就是想讓w之和為1,也就是從A=Chrom(i,j);到 gen=0; 上一行的end這一段循環(huán)使得chrom里每行的7個數(shù)之和為1,共40行 NIND=40; %個體數(shù)目 MAXGEN=500; %最大遺傳代數(shù) NVAR=7; %變量個數(shù) GGAP=0.8; %代溝 %創(chuàng)建區(qū)域描器 FielDR=[0.01,0.01,0.01,0.01,0.01,0.01,0.01;0.9,0.9,0.9,0.9,0.9,0.9,0.9] Chrom=crtrp(NIND,FielDR);%創(chuàng)建初始種群 A=Chrom(i,j); for j=1:40; sum(j)=0; for m=1:7; sum(j)=A(j,m)+sum(j); end; end for j=1:40 for i=1:7 Chrom(i,j)=A(i,j)/sum(j); end end; gen=0; %代計(jì)數(shù)器 ObjV=objfun_1(Chrom); %計(jì)算初始種群個體的目標(biāo)函數(shù)值 while gen FitnV=ranking(ObjV); %分配適應(yīng)度值 SelCh=select('sus',Chrom,FitnV,GGAP); %選擇(隨即遍歷抽樣) SelCh=recombin('xovsp',SelCh,0.7); %重組(單點(diǎn)交叉) SelCh=mutbga(SelCh,FielDR); %變異 ObjVSel=objfun_1(SelCh); %子代函數(shù)值 [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel); %重插入 gen=gen+1; [Y,I]=min(ObjV); %最優(yōu)解,種群序號 Y w=(Chrom(I, ) %輸出自變量值end w ??? Subscript indices must either be real positive integers or logicals. Error in ==> yichuansuanfa at 8 A=Chrom(i,j); |
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑,一志愿北林食品與營養(yǎng)095500,301分,已過六級,有科研經(jīng)歷 +3 | 快樂儲蓄罐 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
[論文投稿]
chinese chemical letters英文版投稿求助
130+4
|
Yishengeryi 2026-03-30 | 4/200 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +6 | hanamiko 2026-03-29 | 6/300 |
|
|
[考研] 材料與化工272求調(diào)劑 +25 | 阿斯蒂芬2004 2026-03-28 | 25/1250 |
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +7 | 朵朵- 2026-03-26 | 9/450 |
|
|
[考研] 生物考研337分求調(diào)劑 +4 | cgxin 2026-03-30 | 6/300 |
|
|
[考研] 347求調(diào)劑 +11 | 山頂見α 2026-03-25 | 11/550 |
|
|
[考研] 一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑 +8 | 披星河 2026-03-30 | 8/400 |
|
|
[考研] 福建理工大學(xué)材料學(xué)院先進(jìn)合金團(tuán)隊(duì)招收考研調(diào)劑學(xué)生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
[考研] 285求調(diào)劑 +6 | AZMK 2026-03-29 | 9/450 |
|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +10 | Koxui 2026-03-28 | 10/500 |
|
|
[考研] 328求調(diào)劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚 2026-03-30 | 5/250 |
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|