| 24小時(shí)熱門(mén)版塊排行榜 |
| 查看: 1378 | 回復(fù): 14 | ||
wkxj新蟲(chóng) (初入文壇)
|
[求助]
matlab中將含有變量“w”的復(fù)雜多項(xiàng)式存入矩陣元素,無(wú)法生成矩陣。哪里出問(wèn)題了?
|
|
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; %空氣縱波的拉梅常數(shù) lml2=290.24e9; %鋼縱波的拉梅常數(shù) lmn1=0; %空氣中沒(méi)有橫波 lmn2=71e9; %鋼橫波的拉梅常數(shù) 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 |
木蟲(chóng) (著名寫(xiě)手)
新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
|
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 |
木蟲(chóng) (著名寫(xiě)手)
新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
|
clc; tic; epssys=1.0e-6; syms w a1=w*0.003; a2=w*1.376e-004; kt2=w*3.315e-004; b2 = kt2; 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); tg_1=A11; tg(n2,n3)=str2double(vpa(tg_1)); end; end; tg |
木蟲(chóng) (著名寫(xiě)手)
新蟲(chóng) (初入文壇)
木蟲(chóng) (著名寫(xiě)手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +10 | Zhao anqi 2026-03-22 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +4 | 元夕元 2026-03-20 | 4/200 |
|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +4 | 幸運(yùn)的醬醬 2026-03-22 | 5/250 |
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
[考研]
|
簡(jiǎn)木ChuFront 2026-03-19 | 8/400 |
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變,給編輯又寫(xiě)郵件了,沒(méi)回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|