版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (1444)
>考研 (153)
>導師招生 (132)
>論文投稿 (40)
>蟲友互識 (34)
>博后之家 (30)
>公派出國 (23)
>基金申請 (22)
>教師之家 (22)
>考博 (22)
>碩博家園 (11)
>找工作 (10)
>論文道賀祈福 (9)
>材料工程 (3)
>簽證指南 (3)
>材料綜合 (2)

返回列表

wkxj

新蟲 (初入文壇)

應助: 1 (幼兒園)
金幣: 17
帖子: 12
在線: 2.9小時
蟲號: 3321359
注冊: 2014-07-13
專業(yè): 聲學

[求助] 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;                                     %空氣縱波的拉梅常數(shù)
lml2=290.24e9;                                        %鋼縱波的拉梅常數(shù)
lmn1=0;                                              %空氣中沒有橫波
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

回復此樓

» 猜你喜歡

284求調劑已經(jīng)有10人回復
一志愿山東大學藥學學碩求調劑已經(jīng)有4人回復
07化學280分求調劑已經(jīng)有4人回復
298-一志愿中國農(nóng)業(yè)大學-求調劑已經(jīng)有12人回復
求材料，環(huán)境專業(yè)調劑已經(jīng)有3人回復
335求調劑已經(jīng)有5人回復
求調劑已經(jīng)有7人回復
一志愿吉大化學322求調劑已經(jīng)有4人回復
環(huán)境學碩288求調劑已經(jīng)有8人回復
341求調劑(一志愿湖南大學070300) 已經(jīng)有6人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

matlab循環(huán)生成的字符，如何隨著逐步放入字符矩陣中？已經(jīng)有6人回復
matlab中數(shù)組不能與矩陣相乘嗎？！已經(jīng)有5人回復
MATLAB子矩陣的問題已經(jīng)有11人回復
Matlab線性矩陣不等式分塊矩陣怎么描述已經(jīng)有4人回復
我用matlab求矩陣的特性向量矩陣，怎么出錯了？已經(jīng)有13人回復
matlab中fmincon的fun是矩陣運算怎么辦？已經(jīng)有5人回復
matlab求解正交矩陣。已經(jīng)有10人回復
用matlab組集剛度矩陣時，總剛度矩陣對應的行列式不為零已經(jīng)有5人回復
求助一個簡單的矩陣行列式問題~~~ 已經(jīng)有12人回復
matlab簡單矩陣運算問題請教已經(jīng)有7人回復
matlab三維數(shù)組矩陣延拓問題已經(jīng)有6人回復
請教matlab用數(shù)組給參數(shù)矩陣賦值的問題已經(jīng)有11人回復
matlab如何生成不同的隨機數(shù) 已經(jīng)有9人回復
matlab循環(huán)計算轉為矩陣 30金幣求助已經(jīng)有7人回復
matlab稀疏矩陣問題求教已經(jīng)有3人回復
如何用matlab求解矩陣系數(shù)的二階微分方程已經(jīng)有13人回復
如何用Matlab畫矩陣的散點圖？已經(jīng)有3人回復
菜鳥求助~matlab怎么樣讓一個1xn的向量中的每個元素是個矩陣？已經(jīng)有11人回復
初學者傷不起啊——關于matlab... 已經(jīng)有12人回復
matlab矩陣創(chuàng)建問題已經(jīng)有10人回復
關于計算機內(nèi)存和matlab產(chǎn)生矩陣能力的關系的困惑已經(jīng)有3人回復
【求助】關于matlab的高維數(shù)組已經(jīng)有6人回復

1樓 2015-04-29 16:39:37

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

信彼南山

木蟲 (著名寫手)

應助: 33 (小學生)
金幣: 4142.9
散金: 1221
紅花: 16
帖子: 1178
在線: 233.5小時
蟲號: 1133529
注冊: 2010-10-27
專業(yè): 導航、制導與傳感技術

bessel函數(shù)應該都是調用的besselmx（這個好像是一個c編譯成的mex）
可能是mex不支持符號運算吧。

程序太長了

贊一下

回復此樓

2樓2015-04-29 22:48:31

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wkxj

新蟲 (初入文壇)

應助: 1 (幼兒園)
金幣: 17
帖子: 12
在線: 2.9小時
蟲號: 3321359
注冊: 2014-07-13
專業(yè): 聲學

引用回帖:

2樓: Originally posted by 信彼南山 at 2015-04-29 22:48:31
bessel函數(shù)應該都是調用的besselmx（這個好像是一個c編譯成的mex）
可能是mex不支持符號運算吧。

程序太長了

其實就是中間代換的幾個比較繁瑣，但是沒有多大的意義，可以忽略。最后得到的矩陣行列式det=0，要解出變量w的�，F(xiàn)在矩陣一直出不來，不知道為什么。

贊一下

回復此樓

3樓2015-04-30 10:15:55

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wkxj

新蟲 (初入文壇)

應助: 1 (幼兒園)
金幣: 17
帖子: 12
在線: 2.9小時
蟲號: 3321359
注冊: 2014-07-13
專業(yè): 聲學

引用回帖:

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

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

信彼南山

木蟲 (著名寫手)

應助: 33 (小學生)
金幣: 4142.9
散金: 1221
紅花: 16
帖子: 1178
在線: 233.5小時
蟲號: 1133529
注冊: 2010-10-27
專業(yè): 導航、制導與傳感技術

問題就在于你的程序太長了，至少要找到問題出在哪么
可以考慮把無關的語句刪掉，找到出問題的那句

贊一下

回復此樓

5樓2015-04-30 20:11:31

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wkxj

新蟲 (初入文壇)

應助: 1 (幼兒園)
金幣: 17
帖子: 12
在線: 2.9小時
蟲號: 3321359
注冊: 2014-07-13
專業(yè): 聲學

引用回帖:

5樓: Originally posted by 信彼南山 at 2015-04-30 20:11:31
問題就在于你的程序太長了，至少要找到問題出在哪么
可以考慮把無關的語句刪掉，找到出問題的那句

我已經(jīng)整理了，該怎么聯(lián)系你

贊一下

回復此樓

6樓2015-04-30 21:13:01

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wkxj

新蟲 (初入文壇)

應助: 1 (幼兒園)
金幣: 17
帖子: 12
在線: 2.9小時
蟲號: 3321359
注冊: 2014-07-13
專業(yè): 聲學

引用回帖:

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

贊一下

回復此樓

7樓2015-04-30 21:16:23

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

信彼南山

木蟲 (著名寫手)

應助: 33 (小學生)
金幣: 4142.9
散金: 1221
紅花: 16
帖子: 1178
在線: 233.5小時
蟲號: 1133529
注冊: 2010-10-27
專業(yè): 導航、制導與傳感技術

剛看了一下
你那句str2double(vpa(tg_1))
str2double的操作要做什么沒看明白

符號變量轉化為數(shù)字？？那結果似乎一定是NAN吧？

贊一下

回復此樓

8樓2015-04-30 23:26:38

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

wkxj

新蟲 (初入文壇)

應助: 1 (幼兒園)
金幣: 17
帖子: 12
在線: 2.9小時
蟲號: 3321359
注冊: 2014-07-13
專業(yè): 聲學

引用回帖:

8樓: Originally posted by 信彼南山 at 2015-04-30 23:26:38
剛看了一下
你那句str2double(vpa(tg_1))
str2double的操作要做什么沒看明白

符號變量轉化為數(shù)字？？那結果似乎一定是NAN吧？

w是一個未知數(shù)，我要生成一個矩陣tg，然后det=0求解出w的值。

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復此樓

9樓2015-05-01 00:12:16

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

信彼南山

木蟲 (著名寫手)

應助: 33 (小學生)
金幣: 4142.9
散金: 1221
紅花: 16
帖子: 1178
在線: 233.5小時
蟲號: 1133529
注冊: 2010-10-27
專業(yè): 導航、制導與傳感技術

[/code]
tg_1=A11;
tg(n2,n3)=str2double(vpa(tg_1));
[/code]
這兩句改一下

CODE:

 tg(n2,n3) = A11;

這樣應該就能創(chuàng)建那個5x5的矩陣了。

不過對這個5x5的矩陣求det的話估計希望不大。

個人覺得應該考慮修改算法。

贊一下

回復此樓

10樓2015-05-01 19:45:19

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主 wkxj 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 材料專業(yè)求調劑 +11	hanamiko 2026-03-18	11/550	2026-03-23 23:12 by peike
[考研] 生物學學碩求調劑 +4	小羊睡著了? 2026-03-23	4/200	2026-03-23 21:42 by hyzs6688
[考研] 0703化學求調劑 +4	奶油草莓. 2026-03-22	5/250	2026-03-23 19:37 by pswait
[考研] 333求調劑 +3	ALULU4408 2026-03-23	3/150	2026-03-23 19:04 by macy2011
[考研] 291求調劑 +5	孅華 2026-03-22	5/250	2026-03-23 09:20 by haoshis
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調劑 +4	幸運的醬醬 2026-03-22	5/250	2026-03-22 20:15 by edmund7
[考研] 生物學071000 329分求調劑 +5	我愛生物生物愛?/a> 2026-03-17	5/250	2026-03-22 16:42 by tcx007
[考研] 318求調劑 +4	plum李子 2026-03-21	7/350	2026-03-22 14:17 by ColorlessPI
[考研] 286分人工智能專業(yè)請求調劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 求調劑 +3	13341 2026-03-20	3/150	2026-03-21 18:28 by 學員8dgXkO
[考研] 278求調劑 +9	煙火先于春 2026-03-17	9/450	2026-03-21 17:47 by 學員8dgXkO
[考研] 307求調劑 +3	余意卿 2026-03-18	3/150	2026-03-21 17:31 by ColorlessPI
[考研] 材料與化工（0856）304求 B區(qū) 調劑 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00
[考研] 一志愿華中科技大學，080502，354分求調劑 +5	守候夕陽CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 一志愿西北大學，070300化學學碩，總分287，雙非一本，求調劑。 +3	晨昏線與星海 2026-03-18	3/150	2026-03-21 00:46 by JourneyLucky
[考研] 330求調劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 304求調劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 317求調劑 +5	申子申申 2026-03-19	9/450	2026-03-20 22:26 by JourneyLucky
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 【同濟軟件】軟件（085405）考研求調劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518

24小時熱門版塊排行榜

[求助] matlab中將含有變量“w”的復雜多項式存入矩陣元素，無法生成矩陣。哪里出問題了？

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

[求助] matlab中將含有變量“w”的復雜多項式存入矩陣元素，無法生成矩陣。哪里出問題了？

» 本主題相關價值貼推薦，對您同樣有幫助: