| 查看: 765 | 回复: 1 | ||||
[求助]
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 报错/////////////??? 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循环，将每一个P的表达式写出，不报错。但是因为N wp值的不确定，所以写成通用模式好修改。但是一直报错 |
回复此楼» 收录本帖的淘帖专辑推荐
 matlab编程绘图 |
» 猜你喜欢
河北省自然科学基金
已经有8人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有5人回复
-
论文撤稿了
已经有5人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
-
26/27申博自荐
已经有9人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
-
揭秘青基评审内幕：几个A才能顺利中标
已经有4人回复
-
青B发送上会通知了吗
已经有7人回复
-
博士申请
已经有3人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68144.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.6小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
小小梨落: 金币+5, ★★★★★最佳答案 2015-09-06 10:47:06
感谢参与，应助指数 +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)写错了，应该是w.*wp(i-1) 吧?（或者w*w(i-1)，不需要用点乘） w被定义为一维的符号变量，w(1)=w，w(2)及以上还未定义。 修改后： 结果： 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楼2015-09-06 10:21:44
10