24小时热门版块排行榜    

查看: 666  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +29 医学老男孩 2026-08-13 61/3050 2026-08-15 20:41 by cqbnu
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +8 majunge000 2026-08-11 10/500 2026-08-15 18:55 by Equinoxhua
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见