ファブリペロー干渉計について

Question

yu 2019-12-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/497767-

编辑： yu 2019-12-27

采纳的回答： Akira Agata

はじめてコードを書いた初心者です

ファブリぺロー干渉計の入力に対する出力を書きたいと思い書きました。

ですが、あまり上手くいきません。

間違えている場合や別の書き方があるならぜひ教えてください。

分かりにくい質問かもしれませんがよろしくお願いします。

式の内容はコメントに追記します

A=constant

R=constant

n=constant

c=constant

%

% This program simulates an Fabry-Perot

% Create an output structure similar to the input

OutputPort1 = InputPort1;

% Length variation

Length = Parameter0;

% calculate the optical signal

if(InputPort1.TypeSignal =='Optical' )

% verify how many sampled signals are in the structure

[ls, cs] = size(InputPort1.Sampled);

if( ls > 0 )

% caculate the at each signal

for counter1=1:cs

OutputPort1.Sampled(1, counter1).Signal = InputPort1.Sampled(1, counter1).Signal * ((1-A-R).^2/(1-R).^2+4*R*sin.^2*(2*pi*n*Length*f/c));

end

4 个评论
显示 2更早的评论隐藏 2更早的评论

Akira Agata 2019-12-26

ファブリペロー干渉計の反射特性は、半透明膜の透過率と反射率にも依存するはずですが、この式にはそれらの項が含まれていないようです。念のため、いまいちど理論式をご確認ください。

yu 2019-12-26

コメント及びご助言ありがとうございます。

おっしゃられた通り、今までの式は恐らく反射膜と半透明膜が同一のものだった場合の式だったようです

P=(((Tx.^2+Rx.^2)*Ry-Rx).^2+4*Rx*Ry(Tx.^2+Rx.^2)*sin.^2*(φ/2)/(1-Rx*Ry) .^2+4*Rx*Ry*sin.^2(φ/2))*│Ei│.^2

φ＝2Lω/ｃ

Ei＝レーザー光振幅　L=反射膜と半透明膜間の距離（変数）Rx,Tx.=半透明膜の反射率及び透過率　Ry,Ty=反射膜の反射率及び透過率　ｃ＝光速　ω＝レーザーの角周波数

恐らくこちらの式が正しいようです。

この式からL（反射膜と半透明膜間の距離）で変化する反射特性をプロットしたいと考えます。

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Akira Agata 2019-12-26

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/497767-#answer_407841

编辑：Akira Agata 2019-12-27

初期パラメータを以下のように想定して計算したところ、波長によっては透過率（＝出力光電力／入力光電力）が 1 を超えるという結果になっています。記載頂いた式のどこかに誤りがあると思われますので、再度ご確認頂けないでしょうか。

% 設定パラメータ
Rx = 0.45;                % 半透明膜の反射率
Tx = 0.45;                % 半透明膜の透過率
Ry = 0.95;                % 反射膜の反射率
L = 100.0e-6;             % 反射膜と半透明膜間の距離 [m]
c = 3.0e8;                % 光速 [1/m]
lambda =...               % 波長範囲 [m]
  linspace(1.5e-6,1.6e-6,1000);
% 設定パラメータからの換算値
f = c./lambda;            % 周波数 [Hz]
omega = 2*pi*f;           % 各周波数 [rad/s]
phi = 2*L*omega/c;
% ファブリペロー干渉計の透過率(Tfp)を計算 (ただしTfp = P_out/P_inと定義)
Tfp =...
  ((Tx^2+Rx^2)*Ry-Rx)^2 +...
  (4*Rx*Ry*(Tx^2+Rx^2)*(sin(phi/2)).^2)./((1-Rx*Ry)^2) +...
  4*Rx*Ry*(sin(phi/2)).^2;
% 計算結果をプロット
figure
plot(lambda*1e6,Tfp)
xlabel('Wavelength [\mum]','FontSize',12)
ylabel('Transmittance','FontSize',12)

3 个评论
显示 1更早的评论隐藏 1更早的评论

yu 2019-12-27

自分の知識不足及び説明不足で手間取らせてしまい大変申し訳ございません。

│Ei│.^2の前の(1-Rx*Ry) .^2+4*Rx*Ry*sin.^2(φ/2))が全て分子でした。

そして前回Pの記述がなくとEiに誤りがありました。

P=反射光強度であり、Ei＝レーザー光の電場振幅でした。

% 設定パラメータ

Rx = 0.45; % 半透明膜の反射率

Tx = 0.45; % 半透明膜の透過率

Ry = 0.95; % 反射膜の反射率

L = 100.0e-6; % 反射膜と半透明膜間の距離 [m]

c = 3.0e8; % 光速 [1/m]

lambda =... % 波長範囲 [m]

linspace(1.5e-6,1.6e-6,1000);

% 設定パラメータからの換算値

f = c./lambda; % 周波数 [Hz]

omega = 2*pi*f; % 各周波数 [rad/s]

phi = 2*L*omega/c;

% ファブリペロー干渉計の透過率(Tfp)を計算 (ただしTfp = P_out/P_inと定義)

Tfp =...

((Tx^2+Rx^2)*Ry-Rx)^2 +...

(4*Rx*Ry*(Tx^2+Rx^2)*(sin(phi/2)).^2)./(((1-Rx*Ry)^2) +...

4*Rx*Ry*(sin(phi/2)).^2);

% 計算結果をプロット

figure

plot(lambda*1e6,Tfp)

xlabel('Wavelength [\mum]','FontSize',12)

ylabel('Transmittance','FontSize',12)

と直してみたところ1を超えない範囲になりました。

横軸をLの変化で見る場合は波長範囲を固定し、Lを変化させるようにすればよろしいのでしょうか

Akira Agata 2019-12-27

编辑：Akira Agata 2019-12-27

ご確認頂き、ありがとうございます。

ここまで来れば、あとはご指摘のとおり、波長範囲を固定してLを変化させれば問題ありません。

forループを使って計算しても良いのですが、せっかくですので行列演算に強いMATLABの特徴を活用してみてください。たとえば以下のようにすると、forループを使わずに計算することができます（lambdaとLを、それぞれ行ベクトルと列ベクトルにするところがポイントです）。

% 設定パラメータ
Rx = 0.45;                      % 半透明膜の反射率
Tx = 0.45;                      % 半透明膜の透過率
Ry = 0.95;                      % 反射膜の反射率
L = linspace(100.0e-6,105e-6)'; % 反射膜と半透明膜間の距離 [m]
c = 3.0e8;                      % 光速 [1/m]
lambda =...                     % 波長範囲 [m]
  linspace(1.5e-6,1.55e-6,1000);
% 設定パラメータからの換算値
f = c./lambda;            % 周波数 [Hz]
omega = 2*pi*f;           % 各周波数 [rad/s]
phi = 2*L*omega./c;
% ファブリペロー干渉計の透過率(Tfp)を計算 (ただしR = P_out/P_inと定義)
Tfp =...
  ((Tx^2+Rx^2)*Ry-Rx)^2 +...
  (4*Rx*Ry*(Tx^2+Rx^2)*(sin(phi/2)).^2)./(((1-Rx*Ry)^2) +...
  4*Rx*Ry*(sin(phi/2)).^2);
% 計算結果をプロット
figure
surf(lambda*1e6,L*1e6,Tfp,...
  'EdgeColor','none')
xlabel('Wavelength [\mum]',...
  'FontSize',12,...
  'Rotation',-20)
ylabel('Cavity length [\mum]',...
  'FontSize',12,...
  'Rotation',45)
zlabel('Transmittance','FontSize',12)
view(30,60)

yu 2019-12-27

本当に大変お世話になりました

forループを使わない計算なども大変参考になりました。

今後もコードの勉強に励んでいきたいと思います。

请先登录，再进行评论。

Answer 2

Hiroumi Mita 2019-12-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/497767-#answer_407439

構造体を使っているようですが、プログラムの基本は

部分から全体、簡単なものから複雑なものへ変化させることなので、基本から進めましょう。

#1. まず、この式の正誤を確認しましょう。

OutputPort1.Sampled(1, counter1).Signal = InputPort1.Sampled(1, counter1).Signal * ((1-A-R).^2/(1-R).^2+4*R*sin.^2*(2*pi*n*Length*f/c));

#2. あっているようであれば

適当な数値を各パラメータにあてはめて期待値がでるか確認しましょう。

3, 関数にします。

4．入出力を構造体にします。

このように、確実に階段を昇っていくことがプログラムのコツです。

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

yu 2019-12-24

ご回答ありがとうございます。

まだおぼつかないですが、言われた通り1つずつこなしていきます。

请先登录，再进行评论。