得到频域波形的坐标数组_PA入门之时域、频域中电压、电流、功率ADS仿真结果后处理...

仿真原理图

结果后处理-频域

结果后处理-时域

理论分析时域

一、控制面板

首先先看控制面板，m6定义交流电压峰值Vpk，m7定义了交流电流峰值Ipk，m4定义了直流电压Vdc，m5定义了直流电流Idc，VI Phase是电流和电压之间的相位差，除以了pi进行了归一化，pichoff mode决定电流小于0会不会夹断。通过写好公式然后再用slider画出即可。

二、电压电流时域波形图

首先第一个图，这个图表示时域的电流电压波形，以及标记了直流为0的那条线。

公式如下：

1. 定义频率，角频率，截至时间，生成行向量画图

ff=1E9，w=2*pi*ff，tstop=1/ff，t=[0::(tstop/200)::tstop]

2. 标记直流，生成一个1x2的数组

tz=[0::tstop*1E9]，Zero=[0::0]

3. 电压波形

V_t=vs(indep(m6)*sin(w*t)+indep(m4),t)，其中indep()是返回数据的自变量indep(m6)=vpk，indep(m4)=Vdc。vs(dependent, independent)表示讲两组变量关联，前面是因变量。用plot_vs(V_t, t*1E9)画出即可。

4. 电流波形

电流波形和电压波形相似，但是根据pinch off的值需要if来判断。同时ph为相位差

ph=m3

I_t_sine=vs(indep(m7)*sin(w*t+ph)+indep(m5),t)

I_t_pinchoff=if(I_t_sine>0)then(I_t_sine)else(0)endif

I_t=if(indep(m8)==0)then(I_t_sine)else(I_t_pinchoff)endif

三、功率图

1. 直流功耗

Pdc=indep(m4)*indep(m5)

2. 复合功耗和交流功耗

V_t_rf=vs(indep(m6)*sin(w*t),t)

I_t_rf_sine=vs(indep(m7)*sin(w*t+ph),t)

Idc_pinchoff=fs(I_t)[0]，把I_t做一个时域到频域的转换，然后取出dc值

I_t_rf_pinchoff=I_t-mag(Idc_pinchoff)，I_t-dc值就是有夹断情况的交流电流

I_t_rf=if(indep(m8)==0)then(I_t_rf_sine)else(I_t_rf_pinchoff)endif，做判断

P_diss_gen_composite=vs((V_t*I_t),t*1E9)，计算复合功耗

P_diss_gen_rf=vs((V_t_rf*I_t_rf),t*1E9)，纯rf交流功率

四、电压电流功耗的极坐标

V_f=fs(V_t,0,3E9)，V_f[1]，变换到频域，取出基波

I_f=fs(I_t,0,3E9)，I_f[1]，变换到频域，取出基波

pspec：This measurement gives a power frequency spectrum in harmonic balance analyses

P_f=(pspec(V_f,0,I_f))，P_f[1]，变换到频域，取出基波

五、负载线

plot_vs(I_t, V_t)

六、功耗与电压电流相位差

Phi=[0::(pi/100)::pi]

Power=vs((m6*m7)*cos(Phi), Phi)，注：Vpk*Ipk*cos(phi)，扫相位范围

Phx=if(indep(m3)>1)then(1-(indep(m3)-1))else(indep(m3))endif，让范围在0-pi内

Pphase=vs((m6*m7)*cos(m3), Phx)，在某个phi上的power，Vpk*Ipk*cos(VI_phase)

七、计算功耗和效率

Pdc1=mag(indep(m4)*Idc_pinchoff)，V_dc*I_dc

Pd_g_composite=integrate(V_t*I_t)/1E-9

Pd_pure_rf=integrate(V_t_rf*I_t_rf)/1E-9

delta=Pd_g_composite-Pd_pure_rf

Efficiency=real(-Pd_pure_rf/Pdc1)*100

因为是理想模型，所以存在物理不纯在的情况

Efficiency1=if(Efficiency<100&&Efficiency>0)then(rightstr(identify_value(int(Efficiency*100)/100),5))else(if(Efficiency<0)then(identify_value(0))else("Non-physical"))endif

仿真的频域和时域的对比

1. FREQ DOMAIN:SELECT ONLY THE RF SIGNAL

Pavg=pspec(V,0,I.i)

Vpeak=V[1]

Ipeak=I.i[1]

Pavg_rf_fund=real(0.5*Vpeak*conj(Ipeak))

2. COMBINED DC + RF POWER

P_t=(ts(V)*ts(I.i))

time=indep(ts(V))

Last_Time_Point=time[sweep_size(time)-1]

P_rf_and_dc_avg=integrate(P_t)/Last_Time_Point

3. ONLY THE RF POWER

因为仿的是HB，所以要用ts：frequency-to-time transform

Ipp=max(ts(I.i))-min(ts(I.i))

Vpp=max(ts(V))-min(ts(V))

P_avg_rf=(Vpp*Ipp)/8，Pavg=1/2*Vp*Ip*cos(phi)