In[]:=
Aatilda=47;Aptilda=0.05;(*Aatilda=48;Aptilda=0.11;*)deltaa=10^(-0.05Aatilda)
Out[]=
0.00446684
In[]:=
deltap=(10^(0.05Aptilda)-1)/(10^(0.05Aptilda)+1)
Out[]=
0.00287822
In[]:=
delta=Min[deltap,deltaa]
Out[]=
0.00287822
In[]:=
Aa=-20Log10[delta]
Out[]=
50.8175
In[]:=
Ap=20Log10[(1+delta)/(1-delta)]
Out[]=
0.05
In[]:=
alpha=0.1102(Aa-8.7)
Out[]=
4.64135
In[]:=
Dv=(Aa-7.95)/14.36(*Dvalueequationdepends*)
Out[]=
2.9852
In[]:=
Nn=85;
In[]:=
Io[x_]:=1+btt[n_]:=alphaSqrt[1-(n/Nn)^2]
Infinity
∑
k=1
2
1
k!
k
x
2
In[]:=
Wk[n_]:=Io[btt[n]]/Io[alpha]
In[]:=
DiscretePlot[Wk[x],{x,-85,85,2},PlotLabel->"Kaiser Window \nα=4.64135\n N=85"]
Out[]=
In[]:=
omegap1=500;omegap2=1050;omegaa1=600;omegaa2=900;omegas=2800;T=2Pi/omegas
Out[]=
π
1400
In[]:=
bt1=omegaa1-omegap1;bt2=omegap2-omegaa2;btm=Min[bt1,bt2]omegaC1=omegap1+btm/2;omegaC2=omegap2-btm/2;H[w_]:=Piecewise[{{1,w<omegaC1},{1,w>omegaC2}}]
Out[]=
100
In[]:=
line1=Line[{{0,1},{500,1}}];line2=Line[{{500,1},{500,0}}];line3=Line[{{600,0},{600,1}}];line4=Line[{{900,0},{900,1}}];line6=Line[{{1050,0},{1050,1}}];line5=Line[{{1050,1},{1500,1}}];lineStyle={Thickness[0.002],Red,Dashed};Plot[H[w],{w,0,1500},Exclusions->None,Epilog->{Directive[lineStyle],line1,line2,line3,line4,line5,line6}]
Out[]=
In[]:=
hantoni[n_]:=Piecewise[{{1+(2(omegaC1-omegaC2))/omegas,n==0},{1/(nPi)(Sin[omegaC1nT]-Sin[omegaC2nT]),True}}]
In[]:=
DiscretePlot[hantoni[n],{n,0,100}]
Out[]=