-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathw_245_270.m
136 lines (136 loc) · 3.01 KB
/
w_245_270.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
n1=[2.67019550348525
2.61254915016153
2.55476136483365
2.49663222354213
2.43799955207791
2.37873387922333
2.31873620098661
2.25793842437397
2.19630653460753
2.13384685644047
2.07061594679781
2.00673484248535
1.94240826764574
1.87794924232453
1.81380794715344
1.75060196336845
1.68914030281478 ];
n2 = [2.68277101839115
2.62575491546820
2.56859358254407
2.51108659934529
2.45307077647103
2.39441489272308
2.33501715324231
2.27480520602112
2.21373874425822
2.15181499725080
2.08907760946066
2.02562955445242
1.96165085221026
1.89742145611209
1.83334880800886
1.76999770055333
1.70811647828944 ];
n3 = [2.69483806136894
2.63843514891231
2.58188411183487
2.52498411955321
2.46757112388420
2.40951238749609
2.35070367941317
2.29106892160114
2.23056224708228
2.16917284515792
2.10693288839289
2.04392930181668
1.98032001588301
1.91635521598013
1.85240351959459
1.78898116284358
1.72677950358653 ];
n4 =[2.70642719615137
2.65062114036387
2.59466492724553
2.53835738824716
2.48153374870663
2.42405997114718
2.36582970753519
2.30676359762843
2.24681089682922
2.18595367776141
2.12421393327399
2.06166429421232
1.99844291333728
1.93477322236579
1.87098862913107
1.80756076182924
1.74512763330745 ];
n5 = [2.71755910423336
2.66233390645350
2.60695723124035
2.55122765142830
2.49497977793245
2.43807844677522
2.38041545913565
2.32190856615602
2.26250271070677
2.20217363694856
2.14093428783813
2.07884451323715
2.01602473865598
1.95267432082047
1.88909467147048
1.82571656987581
1.76312835159412 ];
n6 = [2.72825865936640
2.67359823627010
2.61878568861242
2.56361937405671
2.50793338308880
2.45159158253764
2.39448415270415
2.33652630423481
2.27765916764700
2.21785285460938
2.15711227445257
2.09548583464902
2.03307809437941
1.97006669804942
1.90672412548622
1.84344377484734
1.78076791981932 ];
n1 = n1';
n2 = n2';
n3 = n3';
n4 = n4';
n5 = n5';
n6 = n6';
lambda = 1.2e-6:0.05e-6:2e-6;
C = 3e8;
%omega = 2*pi*C./lambda;
dn1dlambda = [diff(n1)./diff(lambda),0];
ddn1ddlambda = [diff(dn1dlambda)./diff(lambda),0];
D1 = 1e12*1e3*1e-9 .* -lambda ./ C .* ddn1ddlambda;
dn2dlambda = [diff(n2)./diff(lambda),0];
ddn2ddlambda = [diff(dn2dlambda)./diff(lambda),0];
D2 = 1e12*1e3*1e-9 .* -lambda ./ C .* ddn2ddlambda;
dn3dlambda = [diff(n3)./diff(lambda),0];
ddn3ddlambda = [diff(dn3dlambda)./diff(lambda),0];
D3 = 1e12*1e3*1e-9 .* -lambda ./ C .* ddn3ddlambda;
dn4dlambda = [diff(n4)./diff(lambda),0];
ddn4ddlambda = [diff(dn4dlambda)./diff(lambda),0];
D4 = 1e12*1e3*1e-9 .* -lambda ./ C .* ddn4ddlambda;
dn5dlambda = [diff(n5)./diff(lambda),0];
ddn5ddlambda = [diff(dn5dlambda)./diff(lambda),0];
D5 = 1e12*1e3*1e-9 .* -lambda ./ C .* ddn5ddlambda;
dn6dlambda = [diff(n6)./diff(lambda),0];
ddn6ddlambda = [diff(dn6dlambda)./diff(lambda),0];
D6 = 1e12*1e3*1e-9 .* -lambda ./ C .* ddn6ddlambda;
plot(lambda,D1,'r-'),hold on;
plot(lambda,D2,'g-'),hold on;
plot(lambda,D3,'b-'),hold on;
plot(lambda,D4,'c+-'),hold on;
plot(lambda,D5,'ko-'),hold on;
plot(lambda,D6,'m-'),hold on;
axis([1.2e-6,1.9e-6,-2000,2300])