| 2 | 1/1 | 返回列表 |
| 查看: 728 | 回復(fù): 1 | |||||
小小梨落新蟲(chó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á)式寫(xiě)出,不報(bào)錯(cuò)。但是因?yàn)镹 wp值的不確定,所以寫(xiě)成通用模式好修改。但是一直報(bào)錯(cuò) |
matlab編程繪圖 |
主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專(zhuān)家經(jīng)驗(yàn): +1059 |
|
這一句中: 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)寫(xiě)錯(cuò)了,應(yīng)該是w.*wp(i-1) 吧?(或者w*w(i-1),不需要用點(diǎn)乘) w被定義為一維的符號(hào)變量,w(1)=w,w(2)及以上還未定義。 修改后: 結(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) >> |

| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武漢理工材料工程專(zhuān)碩調(diào)劑 +5 | Doleres 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 5/250 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚(yú) 2026-03-18 | 6/300 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專(zhuān)利,省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[基金申請(qǐng)] 被我言中:新模板不強(qiáng)調(diào)格式了,假專(zhuān)家開(kāi)始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|