24小时热门版块排行榜

返回列表

小小梨落

新虫 (小有名气)

应助: 2 (幼儿园)
金币: 27.8
帖子: 86
在线: 35.9小时
虫号: 2125417
注册: 2012-11-14
专业: 常微分方程与动力系统

[求助] 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编程绘图

» 猜你喜欢

假如你的研究生提出不合理要求已经有7人回复
实验室接单子已经有4人回复
全日制（定向）博士已经有4人回复
对氯苯硼酸纯化已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有12人回复
不自信的我已经有12人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复
26申博（荧光探针方向，有机合成）已经有4人回复

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68123.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +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)及以上还未定义。

修改后：

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'

结果：

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)

>>

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

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