版塊導(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í)熱門版塊排行榜

返回列表

小小梨落

新蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 27.8
帖子: 86
在線: 35.9小時(shí)
蟲號(hào): 2125417
注冊(cè): 2012-11-14
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

[求助] matlab程序代碼修改求助已有1人參與

代碼如下：
clear,clc
N=5;
syms w P
wp(1)=-1/1.6886;
wp(2)=1/1.3199;
wp(3)=1/1.7433;
for i=4:N
   wp(i)=0;
end
P(1)=1;
P(2)=w-wp(1);
for i=3:N
   P(i)=-P(i-2).*(1-w.*w(i-1)).^2*((1-wp(i).^2)./(1-wp(i-1).^2)).^(1/2)+P(i-1).*(w-wp(i)+(w-wp(i-1)).*((1-wp(i).^2)./(1-wp(i-1).^2)).^(1/2));
end
報(bào)錯(cuò)/////////////??? Error using ==> mupadmex
Error in MuPAD command: Index exceeds matrix
dimensions.

Error in ==> sym.sym>sym.subsref at 1366
         B =
         mupadmex('mllib::subsref',A.s,inds{:});

Error in ==> Untitled3 at 13
   P(i)=-P(i-2).*(1-w.*w(i-1)).^2*((1-wp(i).^2)./(1-wp(i-1).^2)).^(1/2)+P(i-1).*(w-wp(i)+(w-wp(i-1)).*((1-wp(i).^2)./(1-wp(i-1).^2)).^(1/2));

/////////////不用for循環(huán)，將每一個(gè)P的表達(dá)式寫出，不報(bào)錯(cuò)。但是因?yàn)镹 wp值的不確定，所以寫成通用模式好修改。但是一直報(bào)錯(cuò)

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

matlab編程繪圖

» 猜你喜歡

求調(diào)劑已經(jīng)有3人回復(fù)
265求調(diào)劑已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有6人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有12人回復(fù)
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有3人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù) 已經(jīng)有4人回復(fù)
材料學(xué)碩318求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有6人回復(fù)
081700化工學(xué)碩調(diào)劑已經(jīng)有3人回復(fù)

1樓 2015-09-06 09:28:07

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

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專家經(jīng)驗(yàn): +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時(shí)
蟲號(hào): 1122189
注冊(cè): 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

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

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
小小梨落: 金幣+5, ★★★★★最佳答案 2015-09-06 10:47:06

這一句中：
P(i)=-P(i-2).*(1-w.*w(i-1)).^2*((1-wp(i).^2)./(1-wp(i-1).^2)).^(1/2)+P(i-1).*(w-wp(i)+(w-wp(i-1)).*((1-wp(i).^2)./(1-wp(i-1).^2)).^(1/2));
其中的w.*w(i-1)寫錯(cuò)了，應(yīng)該是w.*wp(i-1) 吧?（或者w*w(i-1)，不需要用點(diǎn)乘）
w被定義為一維的符號(hào)變量，w(1)=w，w(2)及以上還未定義。

修改后：

CODE:

clear all,clc

N=5;

syms w P

wp(1)=-1/1.6886;

wp(2)=1/1.3199;

wp(3)=1/1.7433;

for i=4:N

     wp(i)=0;

end

P(1)=1;

P(2)=w-wp(1);

for i=3:N

      P(i)=-P(i-2).*(1-w.*wp(i-1)).^2*((1-wp(i).^2)./(1-wp(i-1).^2)).^(1/2)+P(i-1).*(w-wp(i)+(w-wp(i-1)).*((1-wp(i).^2)./(1-wp(i-1).^2)).^(1/2));

end

P=P'

結(jié)果：

P =

                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 1
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                               conj(w) + 5000/8443
                                                                                                                                                                                                                                                                                                                                                                                          (conj(w) + 5000/8443)*((2538911247833447*conj(w))/1125899906842624 - 24683612236193007209375/16191719049372408741888) - (1413011340990823*((10000*conj(w))/13199 - 1)^2)/1125899906842624
                                                                                                                                                                                                      ((5000853307591915*conj(w))/2251799813685248 - 1718158433691666875/2453476634498433024)*((conj(w) + 5000/8443)*((2538911247833447*conj(w))/1125899906842624 - 24683612236193007209375/16191719049372408741888) - (1413011340990823*((10000*conj(w))/13199 - 1)^2)/1125899906842624) - (2749053493906667*(conj(w) + 5000/8443)*((10000*conj(w))/17433 - 1)^2)/2251799813685248
(1413011340990823*((10000*conj(w))/13199 - 1)^2)/1125899906842624 - 2*conj(w)*((2749053493906667*(conj(w) + 5000/8443)*((10000*conj(w))/17433 - 1)^2)/2251799813685248 - ((5000853307591915*conj(w))/2251799813685248 - 1718158433691666875/2453476634498433024)*((conj(w) + 5000/8443)*((2538911247833447*conj(w))/1125899906842624 - 24683612236193007209375/16191719049372408741888) - (1413011340990823*((10000*conj(w))/13199 - 1)^2)/1125899906842624)) - (conj(w) + 5000/8443)*((2538911247833447*conj(w))/1125899906842624 - 24683612236193007209375/16191719049372408741888)

>>

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無償私信求助。

2樓2015-09-06 10:21:44

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

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