检查: 640|回复: 0|重视: 0
打印 上一主题 下一主题

[未答复] 求助大神!!过错运用 symengine !!

[仿制链接]

新手

9 麦片

财富积分


050


7

主题

23

帖子

0

最佳答案
跳转到指定楼层
1#
宣布于 2019-4-25 16:27:11 | 只看该作者 回帖奖赏 |倒序阅览 |阅览形式
function [X1,X2] =pl2(N,wx2,q,w,e,r,miga,A)
syms  Z X Y  k ;
Z=floor(t.*N);
t=0:0.1:10;
A1=[1-0.5*wx2*(t-Z/N).^2 t-Z/N 0.5*wx2*(t-Z/N).^2 0];
A2=[0.5*wx2*(t-Z/N).^2 0 1-wx2*(t-Z/N).^2 t-Z/N];
p=[1-wx2/(2*N^2) 1/N wx2/(2*N^2) 0;-wx2/N 1 wx2/N 0;wx2/(2*N^2) 0 1-wx2/N^2 1/N;wx2/N 0 -2*wx2/N 1];
alfa(k)=[A/(2*N^2)*cos(miga*(k-1)/N);A/N*cos(miga*(k-1)/N);0;0];
vo=[q;w;e;r];
Y=p^(Z-k)*alfa(k);
X=symsum(Y,k,1,'Z');
X1=A1*(p^Z*v0+X)+(A/2)*((t-Z/N)^2)*cos(miga*Z/N);
X2=A2*(p^Z*v0+X);
end

过错信息:
>> [X1,X2] =pl2(100,1000,0.5,0.2,0.2,0.1,100,100)
过错运用 symengine
The column index is missing.

犯错 sym/symsum (line 68)
   rSym = mupadmex('symobj::map',fsym.s,'symobj::symsum',x.s,a.s,b.s);

犯错 pl2 (line 11)
X=symsum(Y,k,1,'Z');


求大神辅导!不胜感激!!



您需求登录后才能够回帖 登录 | 注册

本版积分规矩

封闭

站长引荐上一条 /3 下一条

快速回复 回来顶部 回来列表