| 查看: 1691 | 回复: 14 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
matlab中将含有变量“w”的复杂多项式存入矩阵元素，无法生成矩阵。哪里出问题了？
|
||
|
clc; tic; epssys=1.0e-6; syms w %w=1000; a =40e-3; aa1 = a* [1 0 0]; aa2 = a* [0 1 0]; aa3 = a* [0 0 1]; ra11 = (2*pi)*cross(aa2,aa3)/dot(aa1,cross(aa2,aa3)); ra22 = (2*pi)*cross(aa3,aa1)/dot(aa1,cross(aa2,aa3)); ra1 = ra11(1:2); ra2 = ra22(1:2); lml1=0.0001419e9; %空气纵波的拉梅常数 lml2=290.24e9; %钢纵波的拉梅常数 lmn1=0; %空气中没有横波 lmn2=71e9; %钢横波的拉梅常数 rou1=1.29; %空气密度 rou2=7800; %钢密度 a1=w*0.003; a2=w*1.376e-004; b2=w*3.315e-004; r=0.01; bh2 = @(nu,z)besselj(nu,z)-1i*bessely(nu,z); % tg=[]; %tg=num2str(ones(5, 5)); tg=ones(5, 5); for n2=1:5 for n3=1:5 n=n2-3; n1=n3-3; A11=n*besselj(n,a1*r)/r-a1*besselj(n+1,a1*r); B11=n*bh2(n,a1*r)/r-a1*bh2(n+1,a1*r); C11=n*besselj(n,a2*r)/r-a2*besselj(n+1,a2*r); C21=1i*n*besselj(n,b2*r); C31=1i*kz*(n*besselj(n,b2*r)/r-b2*besselj(n+1,b2*r))/kt2; %kt2^2=ω^2*ρ2/μ2 A121=-kz^2*lml1*besselj(n,a1*r); A122=(2*n*lml1*a1+2*lml1*a1)*besselj(n+1,a1*r)/r; A123=lml1*a1^2*besselj(n+2,a1*r); A12=A121-A122+A123; B121=-kz^2*lml1*bh2(n,a1*r); B122=(2*n*lml1*a1+2*lml1*a1)*bh2(n+1,a1*r)/r; B123=lml1*a1^2*bh2(n+2,a1*r); B12=B121-B122+B123; C121=(2*lmn2*n^2-2*lmn2*n-lml2*r^2*kz^2)*besselj(n,a2*r)/r^2; %α2=a2 C122=(4*lmn2*n+2*lmn2+2*n*lml2+2*lml2)*a2*besselj(n+1,a2*r)/r; C123=(2*lmn2*a2^2+lml2*a2^2)*besselj(n+2,a2*r); C12=C121-C122+C123; C221=(2*lmn2*1i*n^2-2*lmn2*1i*n)*besselj(n,b2*r)/(kt2*r^2); %β2=b2 C222=2*lmn2*1i*n*b2*besselj(n+1,b2*r)/(kt2*r); C22=C221-C222; C321=(2*lmn2*1i*kz*(n^2-n))*besselj(n,b2*r)/(kt2*r^2); C322=(2*lmn2*1i*kz*(2*n+1))*b2*besselj(n+1,b2*r)/(kt2*r); C323=2*lmn2*1i*kz*b2^2*besselj(n+2,b2*r)/kt2; C32=C321-C322+C323; C131=2*lmn2*(1i*n^2+n^2)*besselj(n,a2*r)/r^2; C132=2*lmn2*1i*n*a2*besselj(n+1,a2*r)/r; C13=C131-C132; C231=lmn2*(2*n-2*n^2)*besselj(n,b2*r)/r^2; C232=lmn2*2*n*b2*besselj(n+1,b2*r)/r; C233=lmn2*b2^2*besselj(n+2,b2*r); C23=C231+C232+C233; C331=lmn2*(2*n*kz-2*n^2*kz)*besselj(n,b2*r)/(kt2*r^2); C332=lmn2*2*n*kz*b2*besselj(n+1,b2*r)/(kt2*r); C33=C331+C332; C141=lmn2*2*1i*kz*n*besselj(n,a2*r)/r; C142=lmn2*2*1i*kz*a2*besselj(n+1,a2*r); C14=C141+C142; C24=-n*kz*lmn2*besselj(n,b2*r)/r; C341=b2^2*lmn2*n*besselj(n,b2*r)/(kt2*r); C342=b2^3*lmn2*besselj(n+1,b2*r)/kt2; C34=C341-C342; F1=C12*C23*C34; F2=C22*C33*C14; F3=C32*C13*C24; F4=C32*C23*C14; F5=C22*C13*C34; F6=C12*C24*C33; F7=C11*C23*C34; F8=C11*C33*C24; F9=C21*C33*C14; F10=C21*C13*C34; F11=C31*C13*C24; F12=C31*C23*C14; Tn1=A11*(F1+F2+F3-F4-F5-F6); Tn2=B12*(F7-F8+F9-F10+F11-F12); Tn=(Tn1/Tn2); tg_1=Tn; tg(n2,n3)=str2double(vpa(tg_1)); end; end; tg |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有20人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
13楼2015-05-01 22:11:08
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
2楼2015-04-29 22:48:31
3楼2015-04-30 10:15:55
|
clc; tic; epssys=1.0e-6; syms w kl1=w*0.003; kl2=w*1.376e-004; kt2=w*3.315e-004; kz=0; a1 = sqrt(kl1^2-kz^2); a2 = sqrt(kl2^2-kz^2); b2 = sqrt(kt2^2-kz^2); r=0.01; bh2 = @(nu,z)besselj(nu,z)-1i*bessely(nu,z); tg=zeros(5,5); for n2=1:5 for n3=1:5 n=n2-3; n1=n3-3; A11=n*besselj(n,a1*r)/r-a1*besselj(n+1,a1*r); B11=n*bh2(n,a1*r)/r-a1*bh2(n+1,a1*r); C11=n*besselj(n,a2*r)/r-a2*besselj(n+1,a2*r); C21=1i*n*besselj(n,b2*r); A121=0; A122=(2*n*lml1*a1+2*lml1*a1)*besselj(n+1,a1*r)/r; A123=lml1*a1^2*besselj(n+2,a1*r); A12=A121-A122+A123; B121=0; B122=(2*n*lml1*a1+2*lml1*a1)*bh2(n+1,a1*r)/r; B123=lml1*a1^2*bh2(n+2,a1*r); B12=B121-B122+B123; C121=(2*lmn2*n^2-2*lmn2*n)*besselj(n,a2*r)/r^2; C122=(4*lmn2*n+2*lmn2+2*n*lml2+2*lml2)*a2*besselj(n+1,a2*r)/r; C123=(2*lmn2*a2^2+lml2*a2^2)*besselj(n+2,a2*r); C12=C121-C122+C123; C221=(2*lmn2*1i*n^2-2*lmn2*1i*n)*besselj(n,b2*r)/(kt2*r^2); C222=2*lmn2*1i*n*b2*besselj(n+1,b2*r)/(kt2*r); C22=C221-C222; C131=2*lmn2*(1i*n^2+n^2)*besselj(n,a2*r)/r^2; C132=2*lmn2*1i*n*a2*besselj(n+1,a2*r)/r; C13=C131-C132; C231=lmn2*(2*n-2*n^2)*besselj(n,b2*r)/r^2; C232=lmn2*2*n*b2*besselj(n+1,b2*r)/r; C233=lmn2*b2^2*besselj(n+2,b2*r); C23=C231+C232+C233; C341=b2^2*lmn2*n*besselj(n,b2*r)/(kt2*r); C342=b2^3*lmn2*besselj(n+1,b2*r)/kt2; C34=C341-C342; F1=C12*C23*C34; F5=C22*C13*C34; F7=C11*C23*C34; F10=C21*C13*C34; Tn1=A11*(F1-F5)-A12*(F7-F10); Tn2=B12*(F7-F10)-B11*(F1-F5); Tn=vpa(Tn1/Tn2); tg_1=Tn; tg(n2,n3)=str2double(vpa(tg_1)); end; end; tg |
4楼2015-04-30 11:05:03
