| 查看: 485 | 回复: 0 | ||
[求助]
matlab程序求助
|
|
代码如下： 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值的不确定，所以写成通用模式好修改。但是一直报错 求大神指教 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有189人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
