| 查看: 607 | 回复: 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 >> 问大神怎么解决 |
回复此楼» 猜你喜欢
290调剂生物0860
已经有40人回复
-
求调剂推荐
已经有3人回复
-
085400电子信息类（川大控制工程）求调剂可跨专业 求老师联系
已经有6人回复
-
材料工程281还有调剂机会吗
已经有42人回复
-
273求调剂
已经有6人回复
-
药学305求调剂
已经有6人回复
-
求调剂学校
已经有12人回复
-
山东省基金2026
已经有8人回复
-
求助调剂，跨调
已经有19人回复
-
求调剂
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助MATLAB做三维图问题
已经有4人回复
- 关于Matlab画三维曲面图的问题，求助高手，感激感激
已经有3人回复
- 求助一个matlab关于分类讨论的问题
已经有9人回复
- matlab数据拟合问题求助
已经有7人回复
- 求助 Matlab 读取图上数据点 Data Cursor 显示数据问题
已经有6人回复
- 求助关于Matlab中字符串的转化问题
已经有4人回复
- matlab画图求助
已经有11人回复
- 求助 matlab交点问题
已经有19人回复
- matlab 积分求助
已经有8人回复
- matlab问题求助
已经有16人回复
- 求助有关Matlab有约束非线性最优化问题
已经有9人回复
- 求助使用MATLAB计算绿色叶片的覆盖度的问题
已经有24人回复
- 重金求助Matlab图像加噪声问题
已经有8人回复
- 求助matlab-simulink运行时报错问题
已经有3人回复
- 【求助】MATLAB多个mat中 数据画图问题
已经有7人回复
- 【求助】matlab 二次规划的优化的问题
已经有4人回复
- 【求助】问一个关于matlab中legend的问题
已经有11人回复
- 【求助】matlab中保存计算数据的问题
已经有8人回复
3楼2013-04-26 22:18:50
2楼2013-04-26 10:38:06
4楼2013-04-27 10:17:27
