標題: 量子勢阱組 勢壘特征提取matlab代碼 [打印本頁]
作者: cxye1989    時間: 2017-4-10 00:29
標題: 量子勢阱組 勢壘特征提取matlab代碼
量子特征提取matlab代碼

量子勢阱組
量子勢壘特征提取
完整代碼下載：
量子特征提取matlab代碼.rar (33.95 KB, 下載次數(shù): 8)





部分代碼預覽：
  1. %公用參數(shù)
  2. hbar=1;
  3. m=1;
  4. y2=1;
  5. dz=0.00034;
  6. %an=a0+0.00034*n
  7. %vn=1+0.09*v0

  9. %起始量子勢壘V0（55Hz以下）
  10. e0=0:0.0000001:0.0142;
  11. v0=1.2;
  12. a0=0.001;

  14. k0=2*a0/hbar*sqrt(2*m*(v0-e0));
  15. l0=v0^2./(4*e0.*(v0-e0));
  16. t0=1./(1+l0.*sinh(k0).*sinh(k0));

  18. %透射系數(shù)接近1的能量
  19. %E=zeros(1,length(e));
  20. T0=zeros(1,length(e0));
  21. j0=1;
  22. for i0=1:length(e0)
  23.     if t0(1,i0)>0.9
  24.         T0(1,i0)=t0(1,i0);
  25.         E0(1,j0)=i0*0.0000001;
  26.         j0=j0+1;
  27.     end
  28. end

  30. %figure(1);
  31. %plot(e0,t0,e0,y2);
  32. %axis([0 0.0142,0 1.2]);
  33. %title('作為能量函數(shù)的透射系數(shù)');
  34. %xlabel('能量E');
  35. %ylabel('透射系數(shù)T');

  37. %第1量子勢壘V1（55Hz）
  38. e1=0:0.0000001:0.0156;
  39. v1=1.308;
  40. a1=a0+dz;

  42. k1=2*a1/hbar*sqrt(2*m*(v1-e1));
  43. l1=v1^2./(4*e1.*(v1-e1));
  44. t1=1./(1+l1.*sinh(k1).*sinh(k1));

  46. %透射系數(shù)接近1的能量
  47. %E=zeros(1,length(e));
  48. T1=zeros(1,length(e1));
  49. j1=1;
  50. for i1=1:length(e1)
  51.     if t1(1,i1)>0.9
  52.         T1(1,i1)=t1(1,i1);
  53.         E1(1,j1)=i1*0.0000001;
  54.         j1=j1+1;
  55.     end
  56. end

  58. % figure(2);
  59. % plot(e1,t1,e1,y2);
  60. % axis([0 0.0156,0 1.2]);
  61. % title('作為能量函數(shù)的透射系數(shù)');
  62. % xlabel('能量E');
  63. % ylabel('透射系數(shù)T');

  65. %第2量子勢壘V2(58.2705Hz)
  66. e2=0:0.0000001:0.0166;
  67. v2=1.416;
  68. a2=a0+dz*2;

  70. k2=2*a2/hbar*sqrt(2*m*(v2-e2));
  71. l2=v2^2./(4*e2.*(v2-e2));
  72. t2=1./(1+l2.*sinh(k2).*sinh(k2));

  74. %透射系數(shù)接近1的能量
  75. %E=zeros(1,length(e));
  76. T2=zeros(1,length(e2));
  77. j2=1;
  78. for i2=1:length(e2)
  79.     if t2(1,i2)>0.9
  80.         T2(1,i2)=t2(1,i2);
  81.         E2(1,j2)=i2*0.0000001;
  82.         j2=j2+1;
  83.     end
  84. end

  86. % figure(3);
  87. % plot(e2,t2,e2,y2);
  88. % axis([0 0.0166,0 1.2]);
  89. % title('作為能量函數(shù)的透射系數(shù)');
  90. % xlabel('能量E');
  91. % ylabel('透射系數(shù)T');

  93. %第3量子勢壘V3(61.7354Hz)
  94. e3=0:0.0000001:0.0175;
  95. v3=1.524;
  96. a3=a0+dz*3;

  98. k3=2*a3/hbar*sqrt(2*m*(v3-e3));
  99. l3=v3^2./(4*e3.*(v3-e3));
  100. t3=1./(1+l3.*sinh(k3).*sinh(k3));

  102. %透射系數(shù)接近1的能量
  103. %E=zeros(1,length(e));
  104. T3=zeros(1,length(e3));
  105. j3=1;
  106. for i3=1:length(e3)
  107.     if t3(1,i3)>0.9
  108.         T3(1,i3)=t3(1,i3);
  109.         E3(1,j3)=i3*0.0000001;
  110.         j3=j3+1;
  111.     end
  112. end

  114. %第4量子勢壘V4(65.4064Hz)
  115. e4=0:0.0000001:0.0186;
  116. v4=1.632;
  117. a4=a0+dz*4;

  119. k4=2*a4/hbar*sqrt(2*m*(v4-e4));
  120. l4=v4^2./(4*e4.*(v4-e4));
  121. t4=1./(1+l4.*sinh(k4).*sinh(k4));

  123. %透射系數(shù)接近1的能量
  124. %E=zeros(1,length(e));
  125. T4=zeros(1,length(e4));
  126. j4=1;
  127. for i4=1:length(e4)
  128.     if t4(1,i4)>0.9
  129.         T4(1,i4)=t4(1,i4);
  130.         E4(1,j4)=i4*0.0000001;
  131.         j4=j4+1;
  132.     end
  133. end

  135. %第5量子勢壘V5(69.2957Hz)
  136. e5=0:0.0000001:0.0197;
  137. v5=1.740;
  138. a5=a0+dz*5;

  140. k5=2*a5/hbar*sqrt(2*m*(v5-e5));
  141. l5=v5^2./(4*e5.*(v5-e5));
  142. t5=1./(1+l5.*sinh(k5).*sinh(k5));

  144. %透射系數(shù)接近1的能量
  145. %E=zeros(1,length(e));
  146. T5=zeros(1,length(e5));
  147. j5=1;
  148. for i5=1:length(e5)
  149.     if t5(1,i5)>0.9
  150.         T5(1,i5)=t5(1,i5);
  151.         E5(1,j5)=i5*0.0000001;
  152.         j5=j5+1;
  153.     end
  154. end

  156. %第6量子勢壘V6(73.4162Hz)
  157. e6=0:0.0000001:0.0209;
  158. v6=1.848;
  159. a6=a0+dz*6;

  161. k6=2*a6/hbar*sqrt(2*m*(v6-e6));
  162. l6=v6^2./(4*e6.*(v6-e6));
  163. t6=1./(1+l6.*sinh(k6).*sinh(k6));

  165. %透射系數(shù)接近1的能量
  166. %E=zeros(1,length(e));
  167. T6=zeros(1,length(e6));
  168. j6=1;
  169. for i6=1:length(e6)
  170.     if t6(1,i6)>0.9
  171.         T6(1,i6)=t6(1,i6);
  172.         E6(1,j6)=i6*0.0000001;
  173.         j6=j6+1;
  174.     end
  175. end

  177. %第7量子勢壘V7(77.7817Hz)
  178. e7=0:0.0000001:0.0221;
  179. v7=1.956;
  180. a7=a0+dz*7;

  182. k7=2*a7/hbar*sqrt(2*m*(v7-e7));
  183. l7=v7^2./(4*e7.*(v7-e7));
  184. t7=1./(1+l7.*sinh(k7).*sinh(k7));

  186. %透射系數(shù)接近1的能量
  187. %E=zeros(1,length(e));
  188. T7=zeros(1,length(e7));
  189. j7=1;
  190. for i7=1:length(e7)
  191.     if t7(1,i7)>0.9
  192.         T7(1,i7)=t7(1,i7);
  193.         E7(1,j7)=i7*0.0000001;
  194.         j7=j7+1;
  195.     end
  196. end

  198. %第8量子勢壘V8(82.4069Hz)
  199. e8=0:0.0000001:0.0234;
  200. v8=2.064;
  201. a8=a0+dz*8;

  203. k8=2*a8/hbar*sqrt(2*m*(v8-e8));
  204. l8=v8^2./(4*e8.*(v8-e8));
  205. t8=1./(1+l8.*sinh(k8).*sinh(k8));

  207. %透射系數(shù)接近1的能量
  208. %E=zeros(1,length(e));
  209. T8=zeros(1,length(e8));
  210. j8=1;
  211. for i8=1:length(e8)
  212.     if t8(1,i8)>0.9
  213.         T8(1,i8)=t8(1,i8);
  214.         E8(1,j8)=i8*0.0000001;
  215.         j8=j8+1;
  216.     end
  217. end

  219. %第9量子勢壘V9(87.3071Hz)
  220. e9=0:0.0000001:0.0248;
  221. v9=2.172;
  222. a9=a0+dz*9;

  224. k9=2*a9/hbar*sqrt(2*m*(v9-e9));
  225. l9=v9^2./(4*e9.*(v9-e9));
  226. t9=1./(1+l9.*sinh(k9).*sinh(k9));

  228. %透射系數(shù)接近1的能量
  229. %E=zeros(1,length(e));
  230. T9=zeros(1,length(e9));
  231. j9=1;
  232. for i9=1:length(e9)
  233.     if t9(1,i9)>0.9
  234.         T9(1,i9)=t9(1,i9);
  235.         E9(1,j9)=i9*0.0000001;
  236.         j9=j9+1;
  237.     end
  238. end

  240. %第10量子勢壘V10(92.4986Hz)
  241. e10=0:0.0000001:0.0263;
  242. v10=2.280;
  243. a10=a0+dz*10;

  245. k10=2*a10/hbar*sqrt(2*m*(v10-e10));
  246. l10=v10^2./(4*e10.*(v10-e10));
  247. t10=1./(1+l10.*sinh(k10).*sinh(k10));

  249. %透射系數(shù)接近1的能量
  250. %E=zeros(1,length(e));
  251. T10=zeros(1,length(e10));
  252. j10=1;
  253. for i10=1:length(e10)
  254.     if t10(1,i10)>0.9
  255.         T10(1,i10)=t10(1,i10);
  256.         E10(1,j10)=i10*0.0000001;
  257.         j10=j10+1;
  258.     end
  259. end

  261. %第11量子勢壘V11(97.9989Hz)
  262. e11=0:0.0000001:0.0278;
  263. v11=2.388;
  264. a11=a0+dz*11;

  266. k11=2*a11/hbar*sqrt(2*m*(v11-e11));
  267. l11=v11^2./(4*e11.*(v11-e11));
  268. t11=1./(1+l11.*sinh(k11).*sinh(k11));

  270. %透射系數(shù)接近1的能量
  271. %E=zeros(1,length(e));
  272. T11=zeros(1,length(e11));
  273. j11=1;
  274. for i11=1:length(e11)
  275.     if t11(1,i11)>0.9
  276.         T11(1,i11)=t11(1,i11);
  277.         E11(1,j11)=i11*0.0000001;
  278.         j11=j11+1;
  279.     end
  280. end

  282. %第12量子勢壘V12(103.826Hz)
  283. e12=0:0.0000001:0.0295;
  284. v12=2.496;
  285. a12=a0+dz*12;

  287. k12=2*a12/hbar*sqrt(2*m*(v12-e12));
  288. l12=v12^2./(4*e12.*(v12-e12));
  289. t12=1./(1+l12.*sinh(k12).*sinh(k12));

  291. %透射系數(shù)接近1的能量
  292. %E=zeros(1,length(e));
  293. T12=zeros(1,length(e12));
  294. j12=1;
  295. for i12=1:length(e12)
  296.     if t12(1,i12)>0.9
  297.         T12(1,i12)=t12(1,i12);
  298.         E12(1,j12)=i12*0.0000001;
  299.         j12=j12+1;
  300.     end
  301. end

  303. %第13量子勢壘V13(110Hz)
  304. e13=0:0.0000001:0.0313;
  305. v13=2.604;
  306. a13=a0+dz*13;

  308. k13=2*a13/hbar*sqrt(2*m*(v13-e13));
  309. l13=v13^2./(4*e13.*(v13-e13));
  310. t13=1./(1+l13.*sinh(k13).*sinh(k13));

  312. %透射系數(shù)接近1的能量
  313. %E=zeros(1,length(e));
  314. T13=zeros(1,length(e13));
  315. j13=1;
  316. for i13=1:length(e13)
  317.     if t13(1,i13)>0.9
  318.         T13(1,i13)=t13(1,i13);
  319.         E13(1,j13)=i13*0.0000001;
  320.         j13=j13+1;
  321.     end
  322. end

  324. %第14量子勢壘V14(116.541Hz)
  325. e14=0:0.0000001:0.0331;
  326. v14=2.712;
  327. a14=a0+dz*14;

  329. k14=2*a14/hbar*sqrt(2*m*(v14-e14));
  330. l14=v14^2./(4*e14.*(v14-e14));
  331. t14=1./(1+l14.*sinh(k14).*sinh(k14));

  333. %透射系數(shù)接近1的能量
  334. %E=zeros(1,length(e));
  335. T14=zeros(1,length(e14));
  336. j14=1;
  337. for i14=1:length(e14)
  338.     if t14(1,i14)>0.9
  339.         T14(1,i14)=t14(1,i14);
  340.         E14(1,j14)=i14*0.0000001;
  341.         j14=j14+1;
  342.     end
  343. end

  345. %第15量子勢壘V15(123.471Hz)
  346. e15=0:0.0000001:0.0351;
  347. v15=2.82;
  348. a15=a0+dz*15;

  350. k15=2*a15/hbar*sqrt(2*m*(v15-e15));
  351. l15=v15^2./(4*e15.*(v15-e15));
  352. t15=1./(1+l15.*sinh(k15).*sinh(k15));

  354. %透射系數(shù)接近1的能量
  355. %E=zeros(1,length(e));
  356. T15=zeros(1,length(e15));
  357. j15=1;
  358. for i15=1:length(e15)
  359.     if t15(1,i15)>0.9
  360.         T15(1,i15)=t15(1,i15);
  361.         E15(1,j15)=i15*0.0000001;
  362.         j15=j15+1;
  363.     end
  364. end

  366. %第16量子勢壘V16(130.813Hz)
  367. e16=0:0.0000001:0.0372;
  368. v16=2.928;
  369. a16=a0+dz*16;

  371. k16=2*a16/hbar*sqrt(2*m*(v16-e16));
  372. l16=v16^2./(4*e16.*(v16-e16));
  373. t16=1./(1+l16.*sinh(k16).*sinh(k16));

  375. %透射系數(shù)接近1的能量
  376. %E=zeros(1,length(e));
  377. T16=zeros(1,length(e16));
  378. j16=1;
  379. for i16=1:length(e16)
  380.     if t16(1,i16)>0.9
  381.         T16(1,i16)=t16(1,i16);
  382.         E16(1,j16)=i16*0.0000001;
  383.         j16=j16+1;
  384.     end
  385. end

  387. %第17量子勢壘V17(138.591Hz)
  388. e17=0:0.0000001:0.0394;
  389. v17=3.036;
  390. a17=a0+dz*17;

  392. k17=2*a17/hbar*sqrt(2*m*(v17-e17));
  393. l17=v17^2./(4*e17.*(v17-e17));
  394. t17=1./(1+l17.*sinh(k17).*sinh(k17));

  396. %透射系數(shù)接近1的能量
  397. %E=zeros(1,length(e));
  398. T17=zeros(1,length(e17));
  399. j17=1;
  400. for i17=1:length(e17)
  401.     if t17(1,i17)>0.9
  402.         T17(1,i17)=t17(1,i17);
  403.         E17(1,j17)=i17*0.0000001;
  404.         j17=j17+1;
  405.     end
  406. end

  408. %第18量子勢壘V18(146.832Hz)
  409. e18=0:0.0000001:0.0417;
  410. v18=3.144;
  411. a18=a0+dz*18;

  413. k18=2*a18/hbar*sqrt(2*m*(v18-e18));
  414. l18=v18^2./(4*e18.*(v18-e18));
  415. t18=1./(1+l18.*sinh(k18).*sinh(k18));

  417. %透射系數(shù)接近1的能量
  418. %E=zeros(1,length(e));
  419. T18=zeros(1,length(e18));
  420. j18=1;
  421. for i18=1:length(e18)
  422.     if t18(1,i18)>0.9
  423.         T18(1,i18)=t18(1,i18);
  424.         E18(1,j18)=i18*0.0000001;
  425.         j18=j18+1;
  426.     end
  427. end

  429. %第19量子勢壘V19(155.563Hz)
  430. e19=0:0.0000001:0.0442;
  431. v19=3.252;
  432. a19=a0+dz*19;

  434. k19=2*a19/hbar*sqrt(2*m*(v19-e19));
  435. l19=v19^2./(4*e19.*(v19-e19));
  436. t19=1./(1+l19.*sinh(k19).*sinh(k19));

  438. %透射系數(shù)接近1的能量
  439. %E=zeros(1,length(e));
  440. T19=zeros(1,length(e19));
  441. j19=1;
  442. for i19=1:length(e19)
  443.     if t19(1,i19)>0.9
  444.         T19(1,i19)=t19(1,i19);
  445.         E19(1,j19)=i19*0.0000001;
  446.         j19=j19+1;
  447.     end
  448. end

  450. %第20量子勢壘V20(164.814Hz)
  451. e20=0:0.0000001:0.0468;
  452. v20=3.26;
  453. a20=a0+dz*20;

  455. k20=2*a20/hbar*sqrt(2*m*(v20-e20));
  456. l20=v20^2./(4*e20.*(v20-e20));
  457. t20=1./(1+l20.*sinh(k20).*sinh(k20));

  459. %透射系數(shù)接近1的能量
  460. %E=zeros(1,length(e));
  461. T20=zeros(1,length(e20));
  462. j20=1;
  463. for i20=1:length(e20)
  464.     if t20(1,i20)>0.9
  465.         T20(1,i20)=t20(1,i20);
  466.         E20(1,j20)=i20*0.0000001;
  467.         j20=j20+1;
  468.     end
  469. end

  471. %第21量子勢壘V21(174.614Hz)
  472. e21=0:0.0000001:0.0496;
  473. v21=3.468;
  474. a21=a0+dz*21;

  476. k21=2*a21/hbar*sqrt(2*m*(v21-e21));
  477. l21=v21^2./(4*e21.*(v21-e21));
  478. t21=1./(1+l21.*sinh(k21).*sinh(k21));

  480. %透射系數(shù)接近1的能量
  481. %E=zeros(1,length(e));
  482. T21=zeros(1,length(e21));
  483. j21=1;
  484. for i21=1:length(e21)
  485.     if t21(1,i21)>0.9
  486.         T21(1,i21)=t21(1,i21);
  487.         E21(1,j21)=i21*0.0000001;
  488.         j21=j21+1;
  489.     end
  490. end

  492. %第22量子勢壘V22(184.997Hz)
  493. e22=0:0.0000001:0.0526;
  494. v22=3.576;
  495. a22=a0+dz*22;

  497. k22=2*a22/hbar*sqrt(2*m*(v22-e22));
  498. l22=v22^2./(4*e22.*(v22-e22));
  499. t22=1./(1+l22.*sinh(k22).*sinh(k22));

  501. %透射系數(shù)接近1的能量
  502. %E=zeros(1,length(e));
  503. T22=zeros(1,length(e22));
  504. j22=1;
  505. for i22=1:length(e22)
  506.     if t22(1,i22)>0.9
  507.         T22(1,i22)=t22(1,i22);
  508.         E22(1,j22)=i22*0.0000001;
  509.         j22=j22+1;
  510.     end
  511. end

  513. %第23量子勢壘V23(195.998Hz)
  514. e23=0:0.0000001:0.0557;
  515. v23=3.684;
  516. a23=a0+dz*23;

  518. k23=2*a23/hbar*sqrt(2*m*(v23-e23));
  519. l23=v23^2./(4*e23.*(v23-e23));
  520. t23=1./(1+l23.*sinh(k23).*sinh(k23));

  522. %透射系數(shù)接近1的能量
  523. %E=zeros(1,length(e));
  524. T23=zeros(1,length(e23));
  525. j23=1;
  526. for i23=1:length(e23)
  527.     if t23(1,i23)>0.9
  528.         T23(1,i23)=t23(1,i23);
  529.         E23(1,j23)=i23*0.0000001;
  530.         j23=j23+1;
  531.     end
  532. end

  534. %第24量子勢壘V24(207.652Hz)
  535. e24=0:0.0000001:0.059;
  536. v24=3.792;
  537. a24=a0+dz*24;

  539. k24=2*a24/hbar*sqrt(2*m*(v24-e24));
  540. l24=v24^2./(4*e24.*(v24-e24));
  541. t24=1./(1+l24.*sinh(k24).*sinh(k24));

  543. %透射系數(shù)接近1的能量
  544. %E=zeros(1,length(e));
  545. T24=zeros(1,length(e24));
  546. j24=1;
  547. for i24=1:length(e24)
  548.     if t24(1,i24)>0.9
  549.         T24(1,i24)=t24(1,i24);
  550.         E24(1,j24)=i24*0.0000001;
  551.         j24=j24+1;
  552.     end
  553. end

  555. %第25量子勢壘V25(220Hz)
  556. e25=0:0.0000001:0.0625;
  557. v25=3.9;
  558. a25=a0+dz*25;

  560. k25=2*a25/hbar*sqrt(2*m*(v25-e25));
  561. l25=v25^2./(4*e25.*(v25-e25));
  562. t25=1./(1+l25.*sinh(k25).*sinh(k25));

  564. %透射系數(shù)接近1的能量
  565. %E=zeros(1,length(e));
  566. T25=zeros(1,length(e25));
  567. j25=1;
  568. for i25=1:length(e25)
  569.     if t25(1,i25)>0.9
  570.         T25(1,i25)=t25(1,i25);
  571.         E25(1,j25)=i25*0.0000001;
  572.         j25=j25+1;
  573.     end
  574. end

  576. %第26量子勢壘V26(233.082Hz)
  577. e26=0:0.0000001:0.0662;
  578. v26=4.008;
  579. a26=a0+dz*26;

  581. k26=2*a26/hbar*sqrt(2*m*(v26-e26));
  582. l26=v26^2./(4*e26.*(v26-e26));
  583. t26=1./(1+l26.*sinh(k26).*sinh(k26));

  585. %透射系數(shù)接近1的能量
  586. %E=zeros(1,length(e));
  587. T26=zeros(1,length(e26));
  588. j26=1;
  589. for i26=1:length(e26)
  590.     if t26(1,i26)>0.9
  591.         T26(1,i26)=t26(1,i26);
  592.         E26(1,j26)=i26*0.0000001;
  593.         j26=j26+1;
  594.     end
  595. end

  597. %第27量子勢壘V27(246.942Hz)
  598. e27=0:0.0000001:0.0702;
  599. v27=4.116;
  600. a27=a0+dz*27;

  602. k27=2*a27/hbar*sqrt(2*m*(v27-e27));
  603. l27=v27^2./(4*e27.*(v27-e27));
  604. t27=1./(1+l27.*sinh(k27).*sinh(k27));

  606. %透射系數(shù)接近1的能量
  607. %E=zeros(1,length(e));
  608. T27=zeros(1,length(e27));
  609. j27=1;
  610. for i27=1:length(e27)
  611.     if t27(1,i27)>0.9
  612.         T27(1,i27)=t27(1,i27);
  613.         E27(1,j27)=i27*0.0000001;
  614.         j27=j27+1;
  615.     end
  616. end

  618. %第28量子勢壘V28(Hz)
  619. e28=0:0.0000001:0.0743;
  620. v28=4.224;
  621. a28=a0+dz*28;

  623. k28=2*a28/hbar*sqrt(2*m*(v28-e28));
  624. l28=v28^2./(4*e28.*(v28-e28));
  625. t28=1./(1+l28.*sinh(k28).*sinh(k28));

  627. %透射系數(shù)接近1的能量
  628. %E=zeros(1,length(e));
  629. T28=zeros(1,length(e28));
  630. j28=1;
  631. for i28=1:length(e28)
  632.     if t28(1,i28)>0.9
  633.         T28(1,i28)=t28(1,i28);
  634.         E28(1,j28)=i28*0.0000001;
  635.         j28=j28+1;
  636.     end
  637. end

  639. %第29量子勢壘V29(277.183Hz)
  640. e29=0:0.0000001:0.0787;
  641. v29=4.332;
  642. a29=a0+dz*29;

  644. k29=2*a29/hbar*sqrt(2*m*(v29-e29));
  645. l29=v29^2./(4*e29.*(v29-e29));
  646. t29=1./(1+l29.*sinh(k29).*sinh(k29));

  648. %透射系數(shù)接近1的能量
  649. %E=zeros(1,length(e));
  650. T29=zeros(1,length(e29));
  651. j29=1;
  652. for i29=1:length(e29)
  653.     if t29(1,i29)>0.9
  654.         T29(1,i29)=t29(1,i29);
  655.         E29(1,j29)=i29*0.0000001;
  656.         j29=j29+1;
  657.     end
  658. end

  660. %第30量子勢壘V30(293.665Hz)
  661. e30=0:0.0000001:0.0834;
  662. v30=4.44;
  663. a30=a0+dz*30;

  665. k30=2*a30/hbar*sqrt(2*m*(v30-e30));
  666. l30=v30^2./(4*e30.*(v30-e30));
  667. t30=1./(1+l30.*sinh(k30).*sinh(k30));

  669. %透射系數(shù)接近1的能量
  670. %E=zeros(1,length(e));
  671. T30=zeros(1,length(e30));
  672. j30=1;
  673. for i30=1:length(e30)
  674.     if t30(1,i30)>0.9
  675.         T30(1,i30)=t30(1,i30);
  676.         E30(1,j30)=i30*0.0000001;
  677.         j30=j30+1;
  678.     end
  679. end

  681. %第31量子勢壘V31(311.127Hz)
  682. e31=0:0.0000001:0.0884;
  683. v31=4.548;
  684. a31=a0+dz*31;

  686. k31=2*a31/hbar*sqrt(2*m*(v31-e31));
  687. l31=v31^2./(4*e31.*(v31-e31));
  688. t31=1./(1+l31.*sinh(k31).*sinh(k31));

  690. %透射系數(shù)接近1的能量
  691. %E=zeros(1,length(e));
  692. T31=zeros(1,length(e31));
  693. j31=1;
  694. for i31=1:length(e31)
  695.     if t31(1,i31)>0.9
  696.         T31(1,i31)=t31(1,i31);
  697.         E31(1,j31)=i31*0.0000001;
  698.         j31=j31+1;
  699.     end
  700. end

  702. %第32量子勢壘V32(329.628Hz)
  703. e32=0:0.0000001:0.0936;
  704. v32=4.656;
  705. a32=a0+dz*32;

  707. k32=2*a32/hbar*sqrt(2*m*(v32-e32));
  708. l32=v32^2./(4*e32.*(v32-e32));
  709. t32=1./(1+l32.*sinh(k32).*sinh(k32));

  711. %透射系數(shù)接近1的能量
  712. %E=zeros(1,length(e));
  713. T32=zeros(1,length(e32));
  714. j32=1;
  715. for i32=1:length(e32)
  716.     if t32(1,i32)>0.9
  717.         T32(1,i32)=t32(1,i32);
  718.         E32(1,j32)=i32*0.0000001;
  719.         j32=j32+1;
  720.     end
  721. end

  723. %第33量子勢壘V33(349.228Hz)
  724. e33=0:0.0000001:0.0992;
  725. v33=4.764;
  726. a33=a0+dz*33;

  728. k33=2*a33/hbar*sqrt(2*m*(v33-e33));
  729. l33=v33^2./(4*e33.*(v33-e33));
  730. t33=1./(1+l33.*sinh(k33).*sinh(k33));

  732. %透射系數(shù)接近1的能量
  733. %E=zeros(1,length(e));
  734. T33=zeros(1,length(e33));
  735. j33=1;
  736. for i33=1:length(e33)
  737.     if t33(1,i33)>0.9
  738.         T33(1,i33)=t33(1,i33);
  739.         E33(1,j33)=i33*0.0000001;
  740.         j33=j33+1;
  741.     end
  742. end

  744. %第34量子勢壘V34(369.994Hz)
  745. e34=0:0.000001:0.1051;
  746. v34=4.872;
  747. a34=a0+dz*34;

  749. k34=2*a34/hbar*sqrt(2*m*(v34-e34));
  750. l34=v34^2./(4*e34.*(v34-e34));
  751. t34=1./(1+l34.*sinh(k34).*sinh(k34));

  753. %透射系數(shù)接近1的能量
  754. %E=zeros(1,length(e));
  755. T34=zeros(1,length(e34));
  756. j34=1;
  757. for i34=1:length(e34)
  758.     if t34(1,i34)>0.9
  759.         T34(1,i34)=t34(1,i34);
  760.         E34(1,j34)=i34*0.000001;
  761.         j34=j34+1;
  762.     end
  763. end

  765. %第35量子勢壘V35(391.995Hz)
  766. e35=0:0.000001:0.1114;
  767. v35=4.98;
  768. a35=a0+dz*35;

  770. k35=2*a35/hbar*sqrt(2*m*(v35-e35));
  771. l35=v35^2./(4*e35.*(v35-e35));
  772. t35=1./(1+l35.*sinh(k35).*sinh(k35));

  774. %透射系數(shù)接近1的能量
  775. %E=zeros(1,length(e));
  776. T35=zeros(1,length(e35));
  777. j35=1;
  778. for i35=1:length(e35)
  779.     if t35(1,i35)>0.9
  780.         T35(1,i35)=t35(1,i35);
  781.         E35(1,j35)=i35*0.000001;
  782.         j35=j35+1;
  783.     end
  784. end

  786. %第36量子勢壘V36(415.305Hz)
  787. e36=0:0.000001:0.1180;
  788. v36=5.088;
  789. a36=a0+dz*36;

  791. k36=2*a36/hbar*sqrt(2*m*(v36-e36));
  792. l36=v36^2./(4*e36.*(v36-e36));
  793. t36=1./(1+l36.*sinh(k36).*sinh(k36));

  795. %透射系數(shù)接近1的能量
  796. %E=zeros(1,length(e));
  797. T36=zeros(1,length(e36));
  798. j36=1;
  799. for i36=1:length(e36)
  800.     if t36(1,i36)>0.9
  801.         T36(1,i36)=t36(1,i36);
  802.         E36(1,j36)=i36*0.000001;
  803.         j36=j36+1;
  804.     end
  805. end

  807. %第37量子勢壘V37(440Hz)
  808. e37=0:0.000001:0.1250;
  809. v37=5.196;
  810. a37=a0+dz*37;

  812. k37=2*a37/hbar*sqrt(2*m*(v37-e37));
  813. l37=v37^2./(4*e37.*(v37-e37));
  814. t37=1./(1+l37.*sinh(k37).*sinh(k37));

  816. %透射系數(shù)接近1的能量
  817. %E=zeros(1,length(e));
  818. T37=zeros(1,length(e37));
  819. j37=1;
  820. for i37=1:length(e37)
  821.     if t37(1,i37)>0.9
  822.         T37(1,i37)=t37(1,i37);
  823.         E37(1,j37)=i37*0.000001;
  824.         j37=j37+1;
  825.     end
  826. end

  828. %第38量子勢壘V38(466.164Hz)
  829. e38=0:0.000001:0.1324;
  830. v38=5.304;
  831. a38=a0+dz*38;

  833. k38=2*a38/hbar*sqrt(2*m*(v38-e38));
  834. l38=v38^2./(4*e38.*(v38-e38));
  835. t38=1./(1+l38.*sinh(k38).*sinh(k38));

  837. %透射系數(shù)接近1的能量
  838. %E=zeros(1,length(e));
  839. T38=zeros(1,length(e38));
  840. j38=1;
  841. for i38=1:length(e38)
  842.     if t38(1,i38)>0.9
  843.         T38(1,i38)=t38(1,i38);
  844.         E38(1,j38)=i38*0.000001;
  845.         j38=j38+1;
  846.     end
  847. end

  849. %第39量子勢壘V39(493.883Hz)
  850. e39=0:0.000001:0.1403;
  851. v39=5.412;
  852. a39=a0+dz*39;

  854. k39=2*a39/hbar*sqrt(2*m*(v39-e39));
  855. l39=v39^2./(4*e39.*(v39-e39));
  856. t39=1./(1+l39.*sinh(k39).*sinh(k39));

  858. %透射系數(shù)接近1的能量
  859. %E=zeros(1,length(e));
  860. T39=zeros(1,length(e39));
  861. j39=1;
  862. for i39=1:length(e39)
  863.     if t39(1,i39)>0.9
  864.         T39(1,i39)=t39(1,i39);
  865.         E39(1,j39)=i39*0.000001;
  866.         j39=j39+1;
  867.     end
  868. end

  870. %第40量子勢壘V40(523.251Hz)
  871. e40=0:0.000001:0.1487;
  872. v40=5.52;
  873. a40=a0+dz*40;

  875. k40=2*a40/hbar*sqrt(2*m*(v40-e40));
  876. l40=v40^2./(4*e40.*(v40-e40));
  877. t40=1./(1+l40.*sinh(k40).*sinh(k40));

  879. %透射系數(shù)接近1的能量
  880. %E=zeros(1,length(e));
  881. T40=zeros(1,length(e40));
  882. j40=1;
  883. for i40=1:length(e40)
  884.     if t40(1,i40)>0.9
  885.         T40(1,i40)=t40(1,i40);
  886.         E40(1,j40)=i40*0.000001;
  887.         j40=j40+1;
  888.     end
  889. end

  891. %第41量子勢壘V41(554.365Hz)
  892. e41=0:0.000001:0.1575;
  893. v41=5.628;
  894. a41=a0+dz*41;

  896. k41=2*a41/hbar*sqrt(2*m*(v41-e41));
  897. l41=v41^2./(4*e41.*(v41-e41));
  898. t41=1./(1+l41.*sinh(k41).*sinh(k41));

  900. %透射系數(shù)接近1的能量
  901. %E=zeros(1,length(e));
  902. T41=zeros(1,length(e41));
  903. j41=1;
  904. for i41=1:length(e41)
  905.     if t41(1,i41)>0.9
  906.         T41(1,i41)=t41(1,i41);
  907.         E41(1,j41)=i41*0.000001;
  908.         j41=j41+1;
  909.     end
  910. end

  912. %第42量子勢壘V42(587.33Hz)
  913. e42=0:0.000001:0.1669;
  914. v42=5.736;
  915. a42=a0+dz*42;

  917. k42=2*a42/hbar*sqrt(2*m*(v42-e42));
  918. l42=v42^2./(4*e42.*(v42-e42));
  919. t42=1./(1+l42.*sinh(k42).*sinh(k42));

  921. %透射系數(shù)接近1的能量
  922. %E=zeros(1,length(e));
  923. T42=zeros(1,length(e42));
  924. j42=1;
  925. for i42=1:length(e42)
  926.     if t42(1,i42)>0.9
  927.         T42(1,i42)=t42(1,i42);
  928.         E42(1,j42)=i42*0.000001;
  929.         j42=j42+1;
  930.     end
  931. end

  933. %第43量子勢壘V43(622.254Hz)
  934. e43=0:0.000001:0.1768;
  935. v43=5.844;
  936. a43=a0+dz*43;

  938. k43=2*a43/hbar*sqrt(2*m*(v43-e43));
  939. l43=v43^2./(4*e43.*(v43-e43));
  940. t43=1./(1+l43.*sinh(k43).*sinh(k43));

  942. %透射系數(shù)接近1的能量
  943. %E=zeros(1,length(e));
  944. T43=zeros(1,length(e43));
  945. j43=1;
  946. for i43=1:length(e43)
  947.     if t43(1,i43)>0.9
  948.         T43(1,i43)=t43(1,i43);
  949.         E43(1,j43)=i43*0.000001;
  950.         j43=j43+1;
  951.     end
  952. end

  954. %第44量子勢壘V44(659.255Hz)
  955. e44=0:0.000001:0.1873;
  956. v44=5.952;
  957. a44=a0+dz*44;

  959. k44=2*a44/hbar*sqrt(2*m*(v44-e44));
  960. l44=v44^2./(4*e44.*(v44-e44));
  961. t44=1./(1+l44.*sinh(k44).*sinh(k44));

  963. %透射系數(shù)接近1的能量
  964. %E=zeros(1,length(e));
  965. T44=zeros(1,length(e44));
  966. j44=1;
  967. for i44=1:length(e44)
  968.     if t44(1,i44)>0.9
  969.         T44(1,i44)=t44(1,i44);
  970.         E44(1,j44)=i44*0.000001;
  971.         j44=j44+1;
  972.     end
  973. end

  975. %第45量子勢壘V45(698.456Hz)
  976. e45=0:0.000001:0.1984;
  977. v45=6.06;
  978. a45=a0+dz*45;

  980. k45=2*a45/hbar*sqrt(2*m*(v45-e45));
  981. l45=v45^2./(4*e45.*(v45-e45));
  982. t45=1./(1+l45.*sinh(k45).*sinh(k45));

  984. %透射系數(shù)接近1的能量
  985. %E=zeros(1,length(e));
  986. T45=zeros(1,length(e45));
  987. j45=1;
  988. for i45=1:length(e45)
  989.     if t45(1,i45)>0.9
  990.         T45(1,i45)=t45(1,i45);
  991.         E45(1,j45)=i45*0.000001;
  992.         j45=j45+1;
  993.     end
  994. end

  996. %第46量子勢壘V46(739.989Hz)
  997. e46=0:0.000001:0.2102;
  998. v46=6.168;
  999. a46=a0+dz*46;

  1001. k46=2*a46/hbar*sqrt(2*m*(v46-e46));
  1002. l46=v46^2./(4*e46.*(v46-e46));
  1003. t46=1./(1+l46.*sinh(k46).*sinh(k46));

  1005. %透射系數(shù)接近1的能量
  1006. %E=zeros(1,length(e));
  1007. T46=zeros(1,length(e46));
  1008. j46=1;
  1009. for i46=1:length(e46)
  1010.     if t46(1,i46)>0.9
  1011.         T46(1,i46)=t46(1,i46);
  1012.         E46(1,j46)=i46*0.000001;
  1013.         j46=j46+1;
  1014.     end
  1015. end

  1017. %第47量子勢壘V47(783.991Hz)
  1018. e47=0:0.000001:0.2227;
  1019. v47=6.276;
  1020. a47=a0+dz*47;

  1022. k47=2*a47/hbar*sqrt(2*m*(v47-e47));
  1023. l47=v47^2./(4*e47.*(v47-e47));
  1024. t47=1./(1+l47.*sinh(k47).*sinh(k47));

  1026. %透射系數(shù)接近1的能量
  1027. %E=zeros(1,length(e));
  1028. T47=zeros(1,length(e47));
  1029. j47=1;
  1030. for i47=1:length(e47)
  1031.     if t47(1,i47)>0.9
  1032.         T47(1,i47)=t47(1,i47);
  1033.         E47(1,j47)=i47*0.000001;
  1034.         j47=j47+1;
  1035.     end
  1036. end

  1038. %第48量子勢壘V48(830.609Hz)
  1039. e48=0:0.000001:0.2360;
  1040. v48=6.384;
  1041. a48=a0+dz*48;

  1043. k48=2*a48/hbar*sqrt(2*m*(v48-e48));
  1044. l48=v48^2./(4*e48.*(v48-e48));
  1045. t48=1./(1+l48.*sinh(k48).*sinh(k48));

  1047. %透射系數(shù)接近1的能量
  1048. %E=zeros(1,length(e));
  1049. T48=zeros(1,length(e48));
  1050. j48=1;
  1051. for i48=1:length(e48)
  1052.     if t48(1,i48)>0.9
  1053.         T48(1,i48)=t48(1,i48);
  1054.         E48(1,j48)=i48*0.000001;
  1055.         j48=j48+1;
  1056.     end
  1057. end

  1059. %第49量子勢壘V49(880Hz)
  1060. e49=0:0.000001:0.2500;
  1061. v49=6.492;
  1062. a49=a0+dz*49;

  1064. k49=2*a49/hbar*sqrt(2*m*(v49-e49));
  1065. l49=v49^2./(4*e49.*(v49-e49));
  1066. t49=1./(1+l49.*sinh(k49).*sinh(k49));

  1068. %透射系數(shù)接近1的能量
  1069. %E=zeros(1,length(e));
  1070. T49=zeros(1,length(e49));
  1071. j49=1;
  1072. for i49=1:length(e49)
  1073.     if t49(1,i49)>0.9
  1074.         T49(1,i49)=t49(1,i49);
  1075.         E49(1,j49)=i49*0.000001;
  1076.         j49=j49+1;
  1077.     end
  1078. end

  1080. %第50量子勢壘V50(932.328Hz)
  1081. e50=0:0.000001:0.2649;
  1082. v50=6.6;
  1083. a50=a0+dz*50;

  1085. k50=2*a50/hbar*sqrt(2*m*(v50-e50));
  1086. l50=v50^2./(4*e50.*(v50-e50));
  1087. t50=1./(1+l50.*sinh(k50).*sinh(k50));

  1089. %透射系數(shù)接近1的能量
  1090. %E=zeros(1,length(e));
  1091. T50=zeros(1,length(e50));
  1092. j50=1;
  1093. for i50=1:length(e50)
  1094.     if t50(1,i50)>0.9
  1095.         T50(1,i50)=t50(1,i50);
  1096.         E50(1,j50)=i50*0.000001;
  1097.         j50=j50+1;
  1098.     end
  1099. end

  1101. %第51量子勢壘V51(987.767Hz)
  1102. e51=0:0.000001:0.2806;
  1103. v51=6.708;
  1104. a51=a0+dz*51;

  1106. k51=2*a51/hbar*sqrt(2*m*(v51-e51));
  1107. l51=v51^2./(4*e51.*(v51-e51));
  1108. t51=1./(1+l51.*sinh(k51).*sinh(k51));

  1110. %透射系數(shù)接近1的能量
  1111. %E=zeros(1,length(e));
  1112. T51=zeros(1,length(e51));
  1113. j51=1;
  1114. for i51=1:length(e51)
  1115.     if t51(1,i51)>0.9
  1116.         T51(1,i51)=t51(1,i51);
  1117.         E51(1,j51)=i51*0.000001;
  1118.         j51=j51+1;
  1119.     end
  1120. end

  1122. %第52量子勢壘V52(1046.5Hz)
  1123. e52=0:0.000001:0.2973;
  1124. v52=6.816;
  1125. a52=a0+dz*52;

  1127. k52=2*a52/hbar*sqrt(2*m*(v52-e52));
  1128. l52=v52^2./(4*e52.*(v52-e52));
  1129. t52=1./(1+l52.*sinh(k52).*sinh(k52));

  1131. %透射系數(shù)接近1的能量
  1132. %E=zeros(1,length(e));
  1133. T52=zeros(1,length(e52));
  1134. j52=1;
  1135. for i52=1:length(e52)
  1136.     if t52(1,i52)>0.9
  1137.         T52(1,i52)=t52(1,i52);
  1138.         E52(1,j52)=i52*0.000001;
  1139.         j52=j52+1;
  1140.     end
  1141. end

  1143. %第53量子勢壘V53(1108.73Hz)
  1144. e53=0:0.000001:0.3150;
  1145. v53=6.924;
  1146. a53=a0+dz*53;

  1148. k53=2*a53/hbar*sqrt(2*m*(v53-e53));
  1149. l53=v53^2./(4*e53.*(v53-e53));
  1150. t53=1./(1+l53.*sinh(k53).*sinh(k53));

  1152. %透射系數(shù)接近1的能量
  1153. %E=zeros(1,length(e));
  1154. T53=zeros(1,length(e53));
  1155. j53=1;
  1156. for i53=1:length(e53)
  1157.     if t53(1,i53)>0.9
  1158.         T53(1,i53)=t53(1,i53);
  1159.         E53(1,j53)=i53*0.000001;
  1160.         j53=j53+1;
  1161.     end
  1162. end

  1164. %第54量子勢壘V54(1174.66Hz)
  1165. e54=0:0.000001:0.6;
  1166. v54=7.032;
  1167. a54=a0+dz*54;

  1169. k54=2*a54/hbar*sqrt(2*m*(v54-e54));
  1170. l54=v54^2./(4*e54.*(v54-e54));
  1171. t54=1./(1+l54.*sinh(k54).*sinh(k54));

  1173. %透射系數(shù)接近1的能量
  1174. %E=zeros(1,length(e));
  1175. T54=zeros(1,length(e54));
  1176. j54=1;
  1177. for i54=1:length(e54)
  1178.     if t54(1,i54)>0.9
  1179.         T54(1,i54)=t54(1,i54);
  1180.         E54(1,j54)=i54*0.000001;
  1181.         j54=j54+1;
  1182.     end
  1183. end

  1185. %第55量子勢了V55(1244.51Hz)
  1186. e55=0:0.000001:0.6;
  1187. v55=7.14;
  1188. a55=a0+dz*55;

  1190. k55=2*a55/hbar*sqrt(2*m*(v55-e55));
  1191. l55=v55^2./(4*e55.*(v55-e55));
  1192. t55=1./(1+l55.*sinh(k55).*sinh(k55));

  1194. %透射系數(shù)接近1的能量
  1195. %E=zeros(1,length(e));
  1196. T55=zeros(1,length(e55));
  1197. j55=1;
  1198. for i55=1:length(e55)
  1199.     if t55(1,i55)>0.9
  1200.         T55(1,i55)=t55(1,i55);
  1201.         E55(1,j55)=i55*0.000001;
  1202.         j55=j55+1;
  1203.     end
  1204. end

  1206. %第56量子勢壘V56(1318.51Hz)
  1207. e56=0:0.000001:0.7;
  1208. v56=7.248;
  1209. a56=a0+dz*56;

  1211. k56=2*a56/hbar*sqrt(2*m*(v56-e56));
  1212. l56=v56^2./(4*e56.*(v56-e56));
  1213. t56=1./(1+l56.*sinh(k56).*sinh(k56));

  1215. %透射系數(shù)接近1的能量
  1216. %E=zeros(1,length(e));
  1217. T56=zeros(1,length(e56));
  1218. j56=1;
  1219. for i56=1:length(e56)
  1220.     if t56(1,i56)>0.9
  1221.         T56(1,i56)=t56(1,i56);
  1222.         E56(1,j56)=i56*0.000001;
  1223.         j56=j56+1;
  1224.     end
  1225. end

  1227. %第57量子勢壘V57(1396.91Hz)
  1228. e57=0:0.000001:0.7;
  1229. v57=7.356;
  1230. a57=a0+dz*57;

  1232. k57=2*a57/hbar*sqrt(2*m*(v57-e57));
  1233. l57=v57^2./(4*e57.*(v57-e57));
  1234. t57=1./(1+l57.*sinh(k57).*sinh(k57));

  1236. %透射系數(shù)接近1的能量
  1237. %E=zeros(1,length(e));
  1238. T57=zeros(1,length(e57));
  1239. j57=1;
  1240. for i57=1:length(e57)
  1241.     if t57(1,i57)>0.9
  1242.         T57(1,i57)=t57(1,i57);
  1243.         E57(1,j57)=i57*0.000001;
  1244.         j57=j57+1;
  1245.     end
  1246. end

  1248. %第58量子勢壘V58(1479.98Hz)
  1249. e58=0:0.000001:0.7;
  1250. v58=7.464;
  1251. a58=a0+dz*58;

  1253. k58=2*a58/hbar*sqrt(2*m*(v58-e58));
  1254. l58=v58^2./(4*e58.*(v58-e58));
  1255. t58=1./(1+l58.*sinh(k58).*sinh(k58));

  1257. %透射系數(shù)接近1的能量
  1258. %E=zeros(1,length(e));
  1259. T58=zeros(1,length(e58));
  1260. j58=1;
  1261. for i58=1:length(e58)
  1262.     if t58(1,i58)>0.9
  1263.         T58(1,i58)=t58(1,i58);
  1264.         E58(1,j58)=i58*0.000001;
  1265.         j58=j58+1;
  1266.     end
  1267. end

  1269. %第59量子勢壘V59(1567.98)
  1270. e59=0:0.000001:0.7;
  1271. v59=7.572;
  1272. a59=a0+dz*59;

  1274. k59=2*a59/hbar*sqrt(2*m*(v59-e59));
  1275. l59=v59^2./(4*e59.*(v59-e59));
  1276. t59=1./(1+l59.*sinh(k59).*sinh(k59));

  1278. %透射系數(shù)接近1的能量
  1279. %E=zeros(1,length(e));
  1280. T59=zeros(1,length(e59));
  1281. j59=1;
  1282. for i59=1:length(e59)
  1283.     if t59(1,i59)>0.9
  1284.         T59(1,i59)=t59(1,i59);
  1285.         E59(1,j59)=i59*0.000001;
  1286.         j59=j59+1;
  1287.     end
  1288. end

  1290. % %第60量子勢壘V60(1661.22Hz)
  1291. % e60=0:0.000001:0.9;
  1292. % v60=7.68;
  1293. % a60=a0+dz*60;
  1294. %
  1295. % k60=2*a60/hbar*sqrt(2*m*(v60-e60));
  1296. % l60=v60^2./(4*e60.*(v60-e60));
  1297. % t60=1./(1+l60.*sinh(k60).*sinh(k60));
  1298. %
  1299. % %透射系數(shù)接近1的能量
  1300. % %E=zeros(1,length(e));
  1301. % T60=zeros(1,length(e60));
  1302. % j60=1;
  1303. % for i60=1:length(e60)
  1304. %     if t60(1,i60)>0.9
  1305. %         T60(1,i60)=t60(1,i60);
  1306. %         E60(1,j60)=i60*0.000001;
  1307. %         j60=j60+1;
  1308. %     end
  1309. % end
  1310. %
  1311. % %第61量子勢壘V61(1760Hz)
  1312. % e61=0:0.000001:0.9;
  1313. % v61=7.788;
  1314. % a61=a0+dz*61;
  1315. %
  1316. % k61=2*a61/hbar*sqrt(2*m*(v61-e61));
  1317. % l61=v61^2./(4*e61.*(v61-e61));
  1318. % t61=1./(1+l61.*sinh(k61).*sinh(k61));
  1319. %
  1320. % %透射系數(shù)接近1的能量
  1321. % %E=zeros(1,length(e));
  1322. % T61=zeros(1,length(e61));
  1323. % j61=1;
  1324. % for i61=1:length(e61)
  1325. %     if t61(1,i61)>0.9
  1326. %         T61(1,i61)=t61(1,i61);
  1327. %         E61(1,j61)=i61*0.000001;
  1328. %         j61=j61+1;
  1329. %     end
  1330. % end
  1331. %
  1332. % %第62量子勢壘V62(1864.66Hz)
  1333. % e62=0:0.000001:0.9;
  1334. % v62=7.896;
  1335. % a62=a0+dz*62;
  1336. %
  1337. % k62=2*a62/hbar*sqrt(2*m*(v62-e62));
  1338. % l62=v62^2./(4*e62.*(v62-e62));
  1339. % t62=1./(1+l62.*sinh(k62).*sinh(k62));
  1340. %
  1341. % %透射系數(shù)接近1的能量
  1342. % %E=zeros(1,length(e));
  1343. % T62=zeros(1,length(e62));
  1344. % j62=1;
  1345. % for i62=1:length(e62)
  1346. %     if t62(1,i62)>0.9
  1347. %         T62(1,i62)=t62(1,i62);
  1348. %         E62(1,j62)=i62*0.000001;
  1349. %         j62=j62+1;
  1350. %     end
  1351. % end
  1352. %
  1353. % %第63量子勢壘V63(1975.53Hz)
  1354. % e63=0:0.000001:0.9;
  1355. % v63=8.004;
  1356. % a63=a0+dz*63;
  1357. %
  1358. % k63=2*a63/hbar*sqrt(2*m*(v63-e63));
  1359. % l63=v63^2./(4*e63.*(v63-e63));
  1360. % t63=1./(1+l63.*sinh(k63).*sinh(k63));
  1361. %
  1362. % %透射系數(shù)接近1的能量
  1363. % %E=zeros(1,length(e));
  1364. % T63=zeros(1,length(e63));
  1365. % j63=1;
  1366. % for i63=1:length(e63)
  1367. %     if t63(1,i63)>0.9
  1368. %         T63(1,i63)=t63(1,i63);
  1369. %         E63(1,j63)=i63*0.000001;
  1370. %         j63=j63+1;
  1371. %     end
  1372. % end
  1373. %
  1374. % %第64量子勢壘V64(2093Hz)
  1375. % e64=0:0.000001:0.9;
  1376. % v64=8.112;
  1377. % a64=a0+dz*64;
  1378. %
  1379. % k64=2*a64/hbar*sqrt(2*m*(v64-e64));
  1380. % l64=v64^2./(4*e64.*(v64-e64));
  1381. % t64=1./(1+l64.*sinh(k64).*sinh(k64));
  1382. %
  1383. % %透射系數(shù)接近1的能量
  1384. % %E=zeros(1,length(e));
  1385. % T64=zeros(1,length(e64));
  1386. % j64=1;
  1387. % for i64=1:length(e64)
  1388. %     if t64(1,i64)>0.9
  1389. %         T64(1,i64)=t64(1,i64);
  1390. %         E64(1,j64)=i64*0.000001;
  1391. %         j64=j64+1;
  1392. %     end
  1393. % end
  1394. %
  1395. % %第65量子勢壘V65(2217.46Hz)
  1396. % e65=0:0.000001:0.9;
  1397. % v65=8.22;
  1398. % a65=a0+dz*65;
  1399. %
  1400. % k65=2*a65/hbar*sqrt(2*m*(v65-e65));
  1401. % l65=v65^2./(4*e65.*(v65-e65));
  1402. % t65=1./(1+l65.*sinh(k65).*sinh(k65));
  1403. %
  1404. % %透射系數(shù)接近1的能量
  1405. % %E=zeros(1,length(e));
  1406. % T65=zeros(1,length(e65));
  1407. % j65=1;
  1408. % for i65=1:length(e65)
  1409. %     if t65(1,i65)>0.9
  1410. %         T65(1,i65)=t65(1,i65);
  1411. %         E65(1,j65)=i65*0.000001;
  1412. %         j65=j65+1;
  1413. %     end
  1414. % end
  1415. %
  1416. % %第66量子勢壘V66(2349.32Hz)
  1417. % e66=0:0.000001:0.9;
  1418. % v66=8.328;
  1419. % a66=a0+dz*66;
  1420. %
  1421. % k66=2*a66/hbar*sqrt(2*m*(v66-e66));
  1422. % l66=v66^2./(4*e66.*(v66-e66));
  1423. % t66=1./(1+l66.*sinh(k66).*sinh(k66));
  1424. %
  1425. % %透射系數(shù)接近1的能量
  1426. % %E=zeros(1,length(e));
  1427. % T66=zeros(1,length(e66));
  1428. % j66=1;
  1429. % for i66=1:length(e66)
  1430. %     if t66(1,i66)>0.9
  1431. %         T66(1,i66)=t66(1,i66);
  1432. %         E66(1,j66)=i66*0.000001;
  1433. %         j66=j66+1;
  1434. %     end
  1435. % end
  1436. %
  1437. % %第67量子勢壘V67(2489.02Hz)
  1438. % e67=0:0.000001:0.9;
  1439. % v67=8.436;
  1440. % a67=a0+dz*67;
  1441. %
  1442. % k67=2*a67/hbar*sqrt(2*m*(v67-e67));
  1443. % l67=v67^2./(4*e67.*(v67-e67));
  1444. % t67=1./(1+l67.*sinh(k67).*sinh(k67));
  1445. %
  1446. % %透射系數(shù)接近1的能量
  1447. % %E=zeros(1,length(e));
  1448. % T67=zeros(1,length(e67));
  1449. % j67=1;
  1450. % for i67=1:length(e67)
  1451. %     if t67(1,i67)>0.9
  1452. %         T67(1,i67)=t67(1,i67);
  1453. %         E67(1,j67)=i67*0.000001;
  1454. %         j67=j67+1;
  1455. %     end
  1456. % end
  1457. %
  1458. % %第68量子勢壘V68(2637.02Hz)
  1459. % e68=0:0.000001:0.9;
  1460. % v68=8.544;
  1461. % a68=a0+dz*68;
  1462. %
  1463. % k68=2*a68/hbar*sqrt(2*m*(v68-e68));
  1464. % l68=v68^2./(4*e68.*(v68-e68));
  1465. % t68=1./(1+l68.*sinh(k68).*sinh(k68));
  1466. %
  1467. % %透射系數(shù)接近1的能量
  1468. % %E=zeros(1,length(e));
  1469. % T68=zeros(1,length(e68));
  1470. % j68=1;
  1471. % for i68=1:length(e68)
  1472. %     if t68(1,i68)>0.9
  1473. %         T68(1,i68)=t68(1,i68);
  1474. %         E68(1,j68)=i68*0.000001;
  1475. %         j68=j68+1;
  1476. %     end
  1477. % end
  1478. %
  1479. % %第69量子勢壘V69(2793.83Hz)
  1480. % e69=0:0.000001:0.9;
  1481. % v69=8.652;
  1482. % a69=a0+dz*69;
  1483. %
  1484. % k69=2*a69/hbar*sqrt(2*m*(v69-e69));
  1485. % l69=v69^2./(4*e69.*(v69-e69));
  1486. % t69=1./(1+l69.*sinh(k69).*sinh(k69));
  1487. %
  1488. % %透射系數(shù)接近1的能量
  1489. % %E=zeros(1,length(e));
  1490. % T69=zeros(1,length(e69));
  1491. % j69=1;
  1492. % for i69=1:length(e69)
  1493. %     if t69(1,i69)>0.9
  1494. %         T69(1,i69)=t69(1,i69);
  1495. %         E69(1,j69)=i69*0.000001;
  1496. %         j69=j69+1;
  1497. %     end
  1498. % end
  1499. %
  1500. % %第70量子勢壘V70(2959.96Hz)
  1501. % e70=0:0.000001:0.9;
  1502. % v70=8.76;
  1503. % a70=a0+dz*70;
  1504. %
  1505. % k70=2*a70/hbar*sqrt(2*m*(v70-e70));
  1506. % l70=v70^2./(4*e70.*(v70-e70));
  1507. % t70=1./(1+l70.*sinh(k70).*sinh(k70));
  1508. %
  1509. % %透射系數(shù)接近1的能量
  1510. % %E=zeros(1,length(e));
  1511. % T70=zeros(1,length(e70));
  1512. % j70=1;
  1513. % for i70=1:length(e70)
  1514. %     if t70(1,i70)>0.9
  1515. %         T70(1,i70)=t70(1,i70);
  1516. %         E70(1,j70)=i70*0.000001;
  1517. %         j70=j70+1;
  1518. %     end
  1519. % end
  1520. %
  1521. % %第71量子勢壘V71(3135.96Hz);
  1522. % e71=0:0.000001:0.99;
  1523. % v71=8.868;
  1524. % a71=a0+dz*71;
  1525. %
  1526. % k71=2*a71/hbar*sqrt(2*m*(v71-e71));
  1527. % l71=v71^2./(4*e71.*(v71-e71));
  1528. % t71=1./(1+l71.*sinh(k71).*sinh(k71));
  1529. %
  1530. % %透射系數(shù)接近1的能量
  1531. % %E=zeros(1,length(e));
  1532. % T71=zeros(1,length(e71));
  1533. % j71=1;
  1534. % for i71=1:length(e71)
  1535. %     if t71(1,i71)>0.9
  1536. %         T71(1,i71)=t71(1,i71);
  1537. %         E71(1,j71)=i71*0.000001;
  1538. %         j71=j71+1;
  1539. %     end
  1540. % end
  1541. %
  1542. % %第72量子勢壘V72(3322.44Hz)
  1543. % e72=0:0.000001:0.99;
  1544. % v72=8.976;
  1545. % a72=a0+dz*72;
  1546. %
  1547. % k72=2*a72/hbar*sqrt(2*m*(v72-e72));
  1548. % l72=v72^2./(4*e72.*(v72-e72));
  1549. % t72=1./(1+l72.*sinh(k72).*sinh(k72));
  1550. %
  1551. % % 透射系數(shù)接近1的能量
  1552. % %E=zeros(1,length(e));
  1553. % T72=zeros(1,length(e72));
  1554. % j72=1;
  1555. % for i72=1:length(e72)
  1556. %     if t72(1,i72)>0.9
  1557. %         T72(1,i72)=t72(1,i72);
  1558. %         E72(1,j72)=i72*0.000001;
  1559. %         j72=j72+1;
  1560. %     end
  1561. % end
復制代碼






歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1