index fc65eeb..53ffbd2 100644 (file)
@@ -9,7 +9,7 @@
<synopsis>[v]=ell1mag(eps,m1,z)</synopsis>
</refsynopsisdiv>
<refsection>
-        <title>パラメータ</title>
+        <title>引数</title>
<variablelist>
<varlistentry>
<term>eps</term>
<title>説明</title>
<para>
楕円フィルタの二乗振幅を得る際に使用される関数.
-            通常は,<literal>m1=eps*eps/(a*a-1)</literal>です.
-            <literal>s=%sn(z,m1)</literal>として,
+            通常は,<literal>m1=eps*eps/(a*a-1)</literal>です.
+            <literal>s=ellipj(z,m1)</literal>として,
<literal>v=real(ones(z)./(ones(z)+eps*eps*s.*s))</literal>を返します.
</para>
</refsection>
<refsection>
<title>例</title>
-        <programlisting role="example"><![CDATA[
+        <programlisting role="example"><![CDATA[
deff('[alpha,BeTa]=alpha_beta(n,m,m1)',...
-'if 2*int(n/2)==n then, BeTa=K1; else, BeTa=0;end;...
-alpha=%k(1-m1)/%k(1-m);')
+'if 2*int(n/2)==n then, BeTa=K1; else, BeTa=0;end;'+...
+'alpha=%k(1-m1)/%k(1-m);')
epsilon=0.1;A=10;  //リップルパラメータ
m1=(epsilon*epsilon)/(A*A-1);n=5;omegac=6;
m=find_freq(epsilon,A,n);omegar = omegac/sqrt(m)
%k(1-m1)*%k(m)/(%k(m1)*%k(1-m))-n   //チェック...
[alpha,Beta]=alpha_beta(n,m,m1)
-alpha*%asn(1,m)-n*%k(m1)      //チェック
+alpha*delip(1,sqrt(m))-n*%k(m1)      //チェック
sample=0:0.01:20;
//等高線に正の実軸をマップ...
-z=alpha*%asn(sample/omegac,m)+Beta*ones(sample);
+z=alpha*delip(sample/omegac,sqrt(m))+Beta*ones(sample);
plot(sample,ell1mag(epsilon,m1,z))
]]></programlisting>
+        <scilab:image><![CDATA[
+deff('[alpha,BeTa]=alpha_beta(n,m,m1)',...
+'if 2*int(n/2)==n then, BeTa=K1; else, BeTa=0;end;'+...
+'alpha=%k(1-m1)/%k(1-m);')
+epsilon=0.1;A=10;  //ripple parameters
+m1=(epsilon*epsilon)/(A*A-1);n=5;omegac=6;
+m=find_freq(epsilon,A,n);omegar = omegac/sqrt(m)
+%k(1-m1)*%k(m)/(%k(m1)*%k(1-m))-n   //Check...
+[alpha,Beta]=alpha_beta(n,m,m1)
+alpha*delip(1,sqrt(m))-n*%k(m1)      //Check
+samples=0:0.01:20;
+//Now we map the positive real axis into the contour...
+z=alpha*delip(samples/omegac,sqrt(m))+Beta*ones(samples);
+plot(samples,ell1mag(epsilon,m1,z))
+]]>     </scilab:image>
</refsection>
<refsection role="see also">
<title>参照</title>
@@ -77,3 +92,4 @@ plot(sample,ell1mag(epsilon,m1,z))
</simplelist>
</refsection>
</refentry>
+