关注我追更更多通信仿真代码文章目录摘要1. 引言2. 系统模型与信道估计算法2.1 OFDM系统模型2.2 导频结构与LS估计2.3 DFT-LS估计2.4 卡尔曼滤波信道跟踪3. 仿真设计3.1 系统参数3.2 对比算法与性能指标3.3 仿真流程4. 仿真结果与分析4.1 仿真图4.2 部分代码5. 总结摘要正交频分复用OFDM技术凭借其频谱效率高、抗多径能力强等优势成为第四代及第五代移动通信系统的核心物理层技术。信道估计的准确性对于接收端相干解调和系统误码率性能至关重要。本文针对时变多径信道研究了一种基于卡尔曼滤波的导频辅助信道估计与跟踪算法并与传统样条插值LS估计及DFT降噪LS估计进行对比。通过建立一阶自回归AR(1)时变信道模型将信道频率响应建模为状态变量利用导频符号观测值进行递推更新实现信道的动态跟踪。仿真采用16QAM调制、导频间隔为8、符号间相关系数为0.99的慢时变信道。结果表明卡尔曼滤波在低信噪比区域显著降低估计均方误差MSE改善系统误码率BER性能优于两种LS基准算法验证了其在时变OFDM系统中的有效性和鲁棒性。1. 引言随着移动通信和高速数据传输的飞速发展OFDM已成为4G/5G及未来无线系统的核心技术。在OFDM系统中信道估计是接收端均衡和相干解调的前提其精度直接决定系统性能。常用的信道估计方法包括最小二乘LS估计、最小均方误差MMSE估计以及基于变换域的降噪算法。LS估计简单易行但对噪声敏感MMSE估计性能最优但需预知信道统计特性DFT-based LS通过对时域冲激响应截断来抑制噪声是一种实用的改进方案。常用的信道估计方法可分为几类基于导频的最小二乘LS估计计算简单但对噪声敏感最小均方误差MMSE估计性能优异但需要信道统计先验信息且复杂度高基于变换域的DFT-LS估计通过对时域冲激响应截断来抑制噪声是LS的一种有效改进。然而上述方法均为每符号独立处理忽略了信道在时间维度的相关性。在移动通信场景中信道因多普勒效应而随时间变化相邻OFDM符号间的信道响应具有较强相关性。利用这种相关性进行跟踪滤波有望进一步提升估计精度。卡尔曼滤波作为一种最优线性递推估计器能够融合历史状态和当前观测动态估计时变系统状态非常适合信道跟踪。文献[1]最早将卡尔曼滤波应用于OFDM信道估计并验证了其在慢时变信道下的有效性。本文在此基础上设计了一种基于导频和卡尔曼滤波的OFDM信道跟踪方案并引入传统LS和DFT-LS作为对比基准通过仿真全面评估三种算法的MSE和BER性能。2. 系统模型与信道估计算法2.1 OFDM系统模型2.2 导频结构与LS估计2.3 DFT-LS估计2.4 卡尔曼滤波信道跟踪3. 仿真设计3.1 系统参数参数数值FFT 点数 (N)256循环前缀长度 (N_g)16调制方式16QAM导频比例1/8间隔 8 子载波导频符号BPSK±1信道抽头数10时间相关系数 (\alpha)0.99OFDM 符号数3000SNR 范围dB-10 : 5 : 30信道模型AR(1) 时变多径信道估计对比算法传统LS插值、DFT-LS、卡尔曼滤波3.2 对比算法与性能指标3.3 仿真流程4. 仿真结果与分析4.1 仿真图可以看到:仿真结果清晰地展示了三种信道估计算法的性能差异。卡尔曼滤波在低信噪比区域优势最为突出这是因为其能够利用信道时间相关性对噪声进行有效平滑而两种最小二乘算法则因每符号独立估计而受噪声影响较大最差。综上所述三种算法的性能排序为卡尔曼滤波优于基于离散傅里叶变换的最小二乘优于传统插值最小二乘验证了卡尔曼滤波在慢时变信道下的有效性与鲁棒性。然而卡尔曼滤波的优势依赖于模型参数如时间相关系数和噪声协方差与实际信道特性的匹配程度在参数失配时性能增益可能下降而两种最小二乘算法虽性能相对有限但具有实现简单、无需先验信息的优势在实际系统设计中需根据场景需求权衡选择。4.2 部分代码%%%%OFDM信道估计与卡尔曼滤波跟踪仿真代码%%%%对比算法传统LS样条插值、DFT-LS时域截断降噪、卡尔曼滤波%%%%信道模型AR(1)时变多径信道%%%%调制方式16QAM clear all;clc;close all;%%系统参数设置nFFT256;%FFT点数 nSym3000;%OFDM符号数 CP16;%循环前缀长度 M16;%调制阶数16QAM EbN0dB-10:5:30;%信噪比范围dB pr1/8;%导频比例每隔8个子载波插入一个导频%计算每个OFDM符号的数据子载波数和比特数 nDataSC(1-pr)*nFFT;%数据子载波个数 nBitsPerSymnDataSC*log2(M);%每符号比特数%生成导频图案BPSK调制±1 pilotsrepmat((2.*(rand(1,nFFT*pr)0.5)-1),nSym,1);pInd1:(1/pr):nFFT;%导频子载波索引 dataIndsetdiff(1:nFFT,pInd);%数据子载波索引%预分配性能向量 mse_lszeros(1,numel(EbN0dB));%传统LS的MSE mse_dftzeros(1,numel(EbN0dB));%DFT-LS的MSE mse_kalmanzeros(1,numel(EbN0dB));%卡尔曼滤波的MSE ber_lszeros(1,numel(EbN0dB));%传统LS的BER ber_dftzeros(1,numel(EbN0dB));%DFT-LS的BER ber_kalmanzeros(1,numel(EbN0dB));%卡尔曼滤波的BER%%主仿真循环forii1:numel(EbN0dB)%%----------发射机----------%生成随机数据比特 ipBit1.*(rand(1,nSym*nBitsPerSym)0.5);%16QAM调制归一化平均功率为1 ipModqammod(ipBit., M, InputType, bit, UnitAveragePower,true);ipModreshape(ipMod,nDataSC,nSym).;%构建频域OFDM符号数据导频 xFzeros(nSym,nFFT);xF(:,dataInd)ipMod;xF(:,pInd)pilots;%IFFT变换到时域并功率归一化 xt1(nFFT/sqrt(nFFT))*ifft(xF.).;%添加循环前缀 xt[xt1(:,nFFT-CP1:nFFT)xt1];%%----------时变多径信道AR(1)模型----------nTap10;%信道抽头数 rho0.99;%相邻符号相关系数越接近1变化越慢 htzeros(nSym,nTap);%初始抽头复高斯归一化功率ht(1,:)1/sqrt(2)*1/sqrt(nTap)*(randn(1,nTap)1i*randn(1,nTap));%递推生成后续符号的信道冲激响应forkk2:nSymht(kk,:)rho*ht(kk-1,:)sqrt(1-rho^2)*1/sqrt(2)*1/sqrt(nTap)*...(randn(1,nTap)1i*randn(1,nTap));end%计算真实信道频率响应用于MSE计算 hFfft(ht,nFFT,2);%信号与信道卷积考虑多径 xhtzeros(nSym,nTapnFFTCP-1);forjj1:nSymxht(jj,:)conv(ht(jj,:),xt(jj,:));end%将多符号串联成一维向量 xtreshape(xht.,1,nSym*(nFFTCPnTap-1));%加性高斯白噪声AWGN snrEbN0dB(ii);sigmasqrt(1/(2*(10^(snr/10))));%噪声标准差 ntsigma*(randn(1,numel(xt))1i*randn(1,numel(xt)));ytxtnt;%%----------接收机----------%重排为符号矩阵去掉循环前缀和信道拖尾 ytreshape(yt., (nFFTCPnTap-1), nSym).;ytyt(:,CP1:CPnFFT);%FFT变换到频域 yF(sqrt(nFFT)/nFFT)*(fft(yt.)).;%预分配三种算法的信道估计存储 hF_ls_allzeros(nSym,nFFT);%传统LS hF_dft_allzeros(nSym,nFFT);%DFT-LS hF_kalman_allzeros(nSym,nFFT);%卡尔曼滤波%----------卡尔曼滤波参数----------Q(1-rho^2)*mean(abs(hF(:)).^2);%过程噪声协方差由信道变化率决定 Rsigma^2;%观测噪声协方差等于噪声功率 hF_estones(1,nFFT)1i*ones(1,nFFT);%初始状态估计 P_est1;%初始估计误差协方差%用于BER统计的比特向量 ipBitR_LS[];ipBitR_DFT[];ipBitR_Kalman[];%----------逐符号处理----------forkk1:nSym yFRyF(kk,:);%当前符号的频域接收信号%1.传统LS估计仅样条插值无时域降噪%导频处LS估计后用样条插值到所有子载波 yF_est_pil_lsspline(pInd,yFR(pInd)./pilots(1,:),1:nFFT);hF_LSyF_est_pil_ls;%直接作为估计结果hF_ls_all(kk,:)hF_LS;%2.DFT-LS估计插值IFFT截断降噪%导频处LS估计-样条插值-IFFT-保留前CP个样点-FFT yF_est_pil_dftspline(pInd,yFR(pInd)./pilots(1,:),1:nFFT);ht_est(nFFT/sqrt(nFFT))*ifft(yF_est_pil_dft.).;hF_DFT(sqrt(nFFT)/nFFT)*(fft(ht_est(:,1:CP)., nFFT)).;hF_dft_all(kk,:)hF_DFT;%----------均衡与解调----------%传统LS均衡 yF_eq_LSyFR./hF_LS;yMod_LSyF_eq_LS(dataInd);ipBitHat_LSqamdemod(yMod_LS., M, OutputType, bit, UnitAveragePower, true).;ipBitR_LS[ipBitR_LS ipBitHat_LS];%DFT-LS均衡 yF_eq_DFTyFR./hF_DFT;yMod_DFTyF_eq_DFT(dataInd);ipBitHat_DFTqamdemod(yMod_DFT., M, OutputType, bit, UnitAveragePower, true).;ipBitR_DFT[ipBitR_DFT ipBitHat_DFT];end%%----------性能指标计算----------%均方误差MSEmse_ls(ii)mean(mean(abs(hF_ls_all-hF).^2));mse_dft(ii)mean(mean(abs(hF_dft_all-hF).^2));mse_kalman(ii)mean(mean(abs(hF_kalman_all-hF).^2));%误码率BER nErr_lssum(ipBitR_LS~ipBit);nErr_dftsum(ipBitR_DFT~ipBit);nErr_kalmansum(ipBitR_Kalman~ipBit);ber_ls(ii)nErr_ls/(nSym*nBitsPerSym);ber_dft(ii)nErr_dft/(nSym*nBitsPerSym);ber_kalman(ii)nErr_kalman/(nSym*nBitsPerSym);end%%绘图%图1MSE对比 figure;semilogy(EbN0dB,mse_ls,s-,LineWidth,2,DisplayName,传统LS插值);hold on;semilogy(EbN0dB,mse_dft,d-,LineWidth,2,DisplayName,DFT-LS时域截断);semilogy(EbN0dB,mse_kalman,o-,LineWidth,2,DisplayName,卡尔曼滤波);grid on;xlabel(SNR (dB));ylabel(MSE);title(信道估计均方误差MSE对比16QAM,ρ0.99,导频比1/8);legend(Location,southwest);%图2BER对比 figure;semilogy(EbN0dB,ber_ls,s-,LineWidth,2,DisplayName,传统LS插值);hold on;semilogy(EbN0dB,ber_dft,d-,LineWidth,2,DisplayName,DFT-LS时域截断);semilogy(EbN0dB,ber_kalman,o-,LineWidth,2,DisplayName,卡尔曼滤波);grid on;xlabel(SNR (dB));ylabel(BER);title(误码率BER性能对比16QAM,ρ0.99,导频比1/8);legend(Location,southwest);disp(仿真完成);5. 总结本文针对时变OFDM系统仿真验证了一种基于卡尔曼滤波的信道估计与跟踪算法。通过与传统的LS插值估计和DFT-LS降噪估计进行对比结果表明卡尔曼滤波能够有效利用信道的时间相关性显著降低估计均方误差MSE并改善系统误码率BER性能验证了其在慢时变信道下的有效性与鲁棒性。参考文献[1] Zheng Y. A novel channel estimation and tracking method for wireless OFDM systems based on pilots and Kalman filtering[J]. IEEE Transactions on Consumer Electronics, 2003, 49(2): 275-283.完整仿真代码可见往期文章文末VX公众号包含往期博客所有代码所见即所得
关注我追更更多通信仿真代码文章目录摘要1. 引言2. 系统模型与信道估计算法2.1 OFDM系统模型2.2 导频结构与LS估计2.3 DFT-LS估计2.4 卡尔曼滤波信道跟踪3. 仿真设计3.1 系统参数3.2 对比算法与性能指标3.3 仿真流程4. 仿真结果与分析4.1 仿真图4.2 部分代码5. 总结摘要正交频分复用OFDM技术凭借其频谱效率高、抗多径能力强等优势成为第四代及第五代移动通信系统的核心物理层技术。信道估计的准确性对于接收端相干解调和系统误码率性能至关重要。本文针对时变多径信道研究了一种基于卡尔曼滤波的导频辅助信道估计与跟踪算法并与传统样条插值LS估计及DFT降噪LS估计进行对比。通过建立一阶自回归AR(1)时变信道模型将信道频率响应建模为状态变量利用导频符号观测值进行递推更新实现信道的动态跟踪。仿真采用16QAM调制、导频间隔为8、符号间相关系数为0.99的慢时变信道。结果表明卡尔曼滤波在低信噪比区域显著降低估计均方误差MSE改善系统误码率BER性能优于两种LS基准算法验证了其在时变OFDM系统中的有效性和鲁棒性。1. 引言随着移动通信和高速数据传输的飞速发展OFDM已成为4G/5G及未来无线系统的核心技术。在OFDM系统中信道估计是接收端均衡和相干解调的前提其精度直接决定系统性能。常用的信道估计方法包括最小二乘LS估计、最小均方误差MMSE估计以及基于变换域的降噪算法。LS估计简单易行但对噪声敏感MMSE估计性能最优但需预知信道统计特性DFT-based LS通过对时域冲激响应截断来抑制噪声是一种实用的改进方案。常用的信道估计方法可分为几类基于导频的最小二乘LS估计计算简单但对噪声敏感最小均方误差MMSE估计性能优异但需要信道统计先验信息且复杂度高基于变换域的DFT-LS估计通过对时域冲激响应截断来抑制噪声是LS的一种有效改进。然而上述方法均为每符号独立处理忽略了信道在时间维度的相关性。在移动通信场景中信道因多普勒效应而随时间变化相邻OFDM符号间的信道响应具有较强相关性。利用这种相关性进行跟踪滤波有望进一步提升估计精度。卡尔曼滤波作为一种最优线性递推估计器能够融合历史状态和当前观测动态估计时变系统状态非常适合信道跟踪。文献[1]最早将卡尔曼滤波应用于OFDM信道估计并验证了其在慢时变信道下的有效性。本文在此基础上设计了一种基于导频和卡尔曼滤波的OFDM信道跟踪方案并引入传统LS和DFT-LS作为对比基准通过仿真全面评估三种算法的MSE和BER性能。2. 系统模型与信道估计算法2.1 OFDM系统模型2.2 导频结构与LS估计2.3 DFT-LS估计2.4 卡尔曼滤波信道跟踪3. 仿真设计3.1 系统参数参数数值FFT 点数 (N)256循环前缀长度 (N_g)16调制方式16QAM导频比例1/8间隔 8 子载波导频符号BPSK±1信道抽头数10时间相关系数 (\alpha)0.99OFDM 符号数3000SNR 范围dB-10 : 5 : 30信道模型AR(1) 时变多径信道估计对比算法传统LS插值、DFT-LS、卡尔曼滤波3.2 对比算法与性能指标3.3 仿真流程4. 仿真结果与分析4.1 仿真图可以看到:仿真结果清晰地展示了三种信道估计算法的性能差异。卡尔曼滤波在低信噪比区域优势最为突出这是因为其能够利用信道时间相关性对噪声进行有效平滑而两种最小二乘算法则因每符号独立估计而受噪声影响较大最差。综上所述三种算法的性能排序为卡尔曼滤波优于基于离散傅里叶变换的最小二乘优于传统插值最小二乘验证了卡尔曼滤波在慢时变信道下的有效性与鲁棒性。然而卡尔曼滤波的优势依赖于模型参数如时间相关系数和噪声协方差与实际信道特性的匹配程度在参数失配时性能增益可能下降而两种最小二乘算法虽性能相对有限但具有实现简单、无需先验信息的优势在实际系统设计中需根据场景需求权衡选择。4.2 部分代码%%%%OFDM信道估计与卡尔曼滤波跟踪仿真代码%%%%对比算法传统LS样条插值、DFT-LS时域截断降噪、卡尔曼滤波%%%%信道模型AR(1)时变多径信道%%%%调制方式16QAM clear all;clc;close all;%%系统参数设置nFFT256;%FFT点数 nSym3000;%OFDM符号数 CP16;%循环前缀长度 M16;%调制阶数16QAM EbN0dB-10:5:30;%信噪比范围dB pr1/8;%导频比例每隔8个子载波插入一个导频%计算每个OFDM符号的数据子载波数和比特数 nDataSC(1-pr)*nFFT;%数据子载波个数 nBitsPerSymnDataSC*log2(M);%每符号比特数%生成导频图案BPSK调制±1 pilotsrepmat((2.*(rand(1,nFFT*pr)0.5)-1),nSym,1);pInd1:(1/pr):nFFT;%导频子载波索引 dataIndsetdiff(1:nFFT,pInd);%数据子载波索引%预分配性能向量 mse_lszeros(1,numel(EbN0dB));%传统LS的MSE mse_dftzeros(1,numel(EbN0dB));%DFT-LS的MSE mse_kalmanzeros(1,numel(EbN0dB));%卡尔曼滤波的MSE ber_lszeros(1,numel(EbN0dB));%传统LS的BER ber_dftzeros(1,numel(EbN0dB));%DFT-LS的BER ber_kalmanzeros(1,numel(EbN0dB));%卡尔曼滤波的BER%%主仿真循环forii1:numel(EbN0dB)%%----------发射机----------%生成随机数据比特 ipBit1.*(rand(1,nSym*nBitsPerSym)0.5);%16QAM调制归一化平均功率为1 ipModqammod(ipBit., M, InputType, bit, UnitAveragePower,true);ipModreshape(ipMod,nDataSC,nSym).;%构建频域OFDM符号数据导频 xFzeros(nSym,nFFT);xF(:,dataInd)ipMod;xF(:,pInd)pilots;%IFFT变换到时域并功率归一化 xt1(nFFT/sqrt(nFFT))*ifft(xF.).;%添加循环前缀 xt[xt1(:,nFFT-CP1:nFFT)xt1];%%----------时变多径信道AR(1)模型----------nTap10;%信道抽头数 rho0.99;%相邻符号相关系数越接近1变化越慢 htzeros(nSym,nTap);%初始抽头复高斯归一化功率ht(1,:)1/sqrt(2)*1/sqrt(nTap)*(randn(1,nTap)1i*randn(1,nTap));%递推生成后续符号的信道冲激响应forkk2:nSymht(kk,:)rho*ht(kk-1,:)sqrt(1-rho^2)*1/sqrt(2)*1/sqrt(nTap)*...(randn(1,nTap)1i*randn(1,nTap));end%计算真实信道频率响应用于MSE计算 hFfft(ht,nFFT,2);%信号与信道卷积考虑多径 xhtzeros(nSym,nTapnFFTCP-1);forjj1:nSymxht(jj,:)conv(ht(jj,:),xt(jj,:));end%将多符号串联成一维向量 xtreshape(xht.,1,nSym*(nFFTCPnTap-1));%加性高斯白噪声AWGN snrEbN0dB(ii);sigmasqrt(1/(2*(10^(snr/10))));%噪声标准差 ntsigma*(randn(1,numel(xt))1i*randn(1,numel(xt)));ytxtnt;%%----------接收机----------%重排为符号矩阵去掉循环前缀和信道拖尾 ytreshape(yt., (nFFTCPnTap-1), nSym).;ytyt(:,CP1:CPnFFT);%FFT变换到频域 yF(sqrt(nFFT)/nFFT)*(fft(yt.)).;%预分配三种算法的信道估计存储 hF_ls_allzeros(nSym,nFFT);%传统LS hF_dft_allzeros(nSym,nFFT);%DFT-LS hF_kalman_allzeros(nSym,nFFT);%卡尔曼滤波%----------卡尔曼滤波参数----------Q(1-rho^2)*mean(abs(hF(:)).^2);%过程噪声协方差由信道变化率决定 Rsigma^2;%观测噪声协方差等于噪声功率 hF_estones(1,nFFT)1i*ones(1,nFFT);%初始状态估计 P_est1;%初始估计误差协方差%用于BER统计的比特向量 ipBitR_LS[];ipBitR_DFT[];ipBitR_Kalman[];%----------逐符号处理----------forkk1:nSym yFRyF(kk,:);%当前符号的频域接收信号%1.传统LS估计仅样条插值无时域降噪%导频处LS估计后用样条插值到所有子载波 yF_est_pil_lsspline(pInd,yFR(pInd)./pilots(1,:),1:nFFT);hF_LSyF_est_pil_ls;%直接作为估计结果hF_ls_all(kk,:)hF_LS;%2.DFT-LS估计插值IFFT截断降噪%导频处LS估计-样条插值-IFFT-保留前CP个样点-FFT yF_est_pil_dftspline(pInd,yFR(pInd)./pilots(1,:),1:nFFT);ht_est(nFFT/sqrt(nFFT))*ifft(yF_est_pil_dft.).;hF_DFT(sqrt(nFFT)/nFFT)*(fft(ht_est(:,1:CP)., nFFT)).;hF_dft_all(kk,:)hF_DFT;%----------均衡与解调----------%传统LS均衡 yF_eq_LSyFR./hF_LS;yMod_LSyF_eq_LS(dataInd);ipBitHat_LSqamdemod(yMod_LS., M, OutputType, bit, UnitAveragePower, true).;ipBitR_LS[ipBitR_LS ipBitHat_LS];%DFT-LS均衡 yF_eq_DFTyFR./hF_DFT;yMod_DFTyF_eq_DFT(dataInd);ipBitHat_DFTqamdemod(yMod_DFT., M, OutputType, bit, UnitAveragePower, true).;ipBitR_DFT[ipBitR_DFT ipBitHat_DFT];end%%----------性能指标计算----------%均方误差MSEmse_ls(ii)mean(mean(abs(hF_ls_all-hF).^2));mse_dft(ii)mean(mean(abs(hF_dft_all-hF).^2));mse_kalman(ii)mean(mean(abs(hF_kalman_all-hF).^2));%误码率BER nErr_lssum(ipBitR_LS~ipBit);nErr_dftsum(ipBitR_DFT~ipBit);nErr_kalmansum(ipBitR_Kalman~ipBit);ber_ls(ii)nErr_ls/(nSym*nBitsPerSym);ber_dft(ii)nErr_dft/(nSym*nBitsPerSym);ber_kalman(ii)nErr_kalman/(nSym*nBitsPerSym);end%%绘图%图1MSE对比 figure;semilogy(EbN0dB,mse_ls,s-,LineWidth,2,DisplayName,传统LS插值);hold on;semilogy(EbN0dB,mse_dft,d-,LineWidth,2,DisplayName,DFT-LS时域截断);semilogy(EbN0dB,mse_kalman,o-,LineWidth,2,DisplayName,卡尔曼滤波);grid on;xlabel(SNR (dB));ylabel(MSE);title(信道估计均方误差MSE对比16QAM,ρ0.99,导频比1/8);legend(Location,southwest);%图2BER对比 figure;semilogy(EbN0dB,ber_ls,s-,LineWidth,2,DisplayName,传统LS插值);hold on;semilogy(EbN0dB,ber_dft,d-,LineWidth,2,DisplayName,DFT-LS时域截断);semilogy(EbN0dB,ber_kalman,o-,LineWidth,2,DisplayName,卡尔曼滤波);grid on;xlabel(SNR (dB));ylabel(BER);title(误码率BER性能对比16QAM,ρ0.99,导频比1/8);legend(Location,southwest);disp(仿真完成);5. 总结本文针对时变OFDM系统仿真验证了一种基于卡尔曼滤波的信道估计与跟踪算法。通过与传统的LS插值估计和DFT-LS降噪估计进行对比结果表明卡尔曼滤波能够有效利用信道的时间相关性显著降低估计均方误差MSE并改善系统误码率BER性能验证了其在慢时变信道下的有效性与鲁棒性。参考文献[1] Zheng Y. A novel channel estimation and tracking method for wireless OFDM systems based on pilots and Kalman filtering[J]. IEEE Transactions on Consumer Electronics, 2003, 49(2): 275-283.完整仿真代码可见往期文章文末VX公众号包含往期博客所有代码所见即所得
