|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
M=8; N=5;
kesai=0.1:1/10:10;
for s=1:M
for k=1:N
Bessel1(k,:)=besselj(1,(k/N*kesai));
Bessel2(k,:)= besselj(1,((k-1)/N*kesai));
Amedium_k(s,2:101)=2./kesai*exp(i*1/3*pi).*[(k/N)*Bessel1(k,:)-((k-1)/N)*Bessel2(k,:)];
syms x
A=2/x*exp(i*1/3*pi)*[(k/N)*besselj(1,((k/N)*x))-((k-1)/N)*besselj(1,((k-1)/N)*x)];
Amedium_k(s,1)=subs(limit(A,x,0,'right'))
A(s,:)=A(s,:)+Amedium_k(s,:);
end;
I(s,:)=abs(A(s,:)).^2;
end;
??? Assignment has more non-singleton rhs dimensions than non-singleton
subscripts.
Error in ==> zehnfufenbu at 11
A(s,:)=A(s,:)+Amedium_k(s,:); |
|