24小时热门版块排行榜    

查看: 573  |  回复: 3

1017199818

铜虫 (初入文坛)

[求助] MATLAB问题求助

function F=TDCL(N,L)
%function A=TDCL1(N,ka,b,s1,s2,u,t)
for NUM=1:5
   if NUM==1
        s1=double(1.0);
    elseif NUM==2
        s2=double(0.5);
    elseif NUM==3
        s3=double(0.2);
    elseif NUM==4
        s4=double(0.1);
    elseif NUM==5
        s5=double(0.0);
   end
   ka=0;
   b=double(1.10);
   %s1=double(1.0);
   u=double(0.2731);
   s2=0;
   t=2.20;
   alpha=(0.25*(ka^2-s1^2+s2^2)^2+(u^(-2)-s1*s2)^2+0.50*(ka^2-s1^2+s2^2))^0.50;
   beta=(0.25*(ka^2-s1^2+s2^2)^2+(u^(-2)-s1*s2)^2-0.50*(ka^2-s1^2+s2^2))^0.50;
   Ka=alpha+i*beta;
   Kb=Ka*b;
   M=N+1;
   n=2*N+1;
   for J=1:51
       x=-b;
       rp=4*(J-1)/50;
       q(J)=4*(J-1)/50;
       Kr=Ka*(x^2+rp^2)^0.50;
       ta=pi+atan(rp/x);
       ta1=atan(rp/b);
   for J=1:L
    if J==1
       for I=1:n
           A(I,J)=(i*Ka-s1)*cos(ta)*exp(i*Ka*b)*exp((i*Ka-s1)*cos(ta))/(exp(-s1*cos(ta))*(s1*cos(ta)*qiubesselh(I-M,Ka))-(Ka*((I-M)*qiubesselh(I-M-1,Ka)-(I-M+1)*qiubesselh(I-M+1,Ka))*legendre(n,cos(ta))/(2*(I-M)+1)));
       end
这个程序提示错误为??? Subscripted assignment dimension mismatch.

Error in ==> TDCL at 37
           A(I,J)=(i*Ka-s1)*cos(ta)*exp(i*Ka*b)*exp((i*Ka-s1)*cos(ta))/(exp(-s1*cos(ta))*(s1*cos(ta)*besselh(I-M,Ka))-(Ka*((I-M)*besselh(I-M-1,Ka)-(I-M+1)*besselh(I-M+1,Ka))*legendre(n,cos(ta))/(2*(I
>> 问大神怎么解决
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

同心同德同舟楫,济人济事济天下。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mvpyqz

至尊木虫 (知名作家)

【答案】应助回帖


感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2013-04-26 23:33:19
貌似执行一次循环之后,A维数在变化,所以导致错误,贴出来的程序不全,只能猜测了!

[ 发自手机版 http://muchong.com/3g ]
2楼2013-04-26 10:38:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mocito

铜虫 (初入文坛)

【答案】应助回帖


感谢参与,应助指数 +1
jjdg: 金币+1, 感谢支招 2013-04-26 23:33:29
把A里面调用的几个function   besselh(),qiubesselh(),legendre(),还有N,L,i 都换成1之后测试可以运算。。。
3楼2013-04-26 22:18:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hua_007

铜虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2013-04-27 16:25:32
提示的直面意思是,维数不匹配,
江湖
4楼2013-04-27 10:17:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1017199818 的主题更新
信息提示
请填处理意见