解答如下
clear;
clc;
lamda=1550e-9;%1550nm
fs=500e6;%500MHz
T=400e-9;%400ns
w=sqrt(log(2))/(log(2)*pi*T*fs);
M=T*fs;
N=[100,1000,10000];
for i=1:length(N)
N1=N(i);
CNR1=linspace(-40,10,50);%载噪比其实就是信噪比
CNR=10.^(CNR1/10);
V_cr1=zeros(1,50);
qianzhixishu1=lamda*fs*w/(2*sqrt(N1*M));%前值系数，就是公式中积分大括号之前的部分
for index=1:50
houzhixishu=@(f)(f./w).^2./(1+(CNR(index)./(sqrt(2*pi).*w).*exp(-f.^2./(2.*w.^2)) ).^-1).^2;%后置系数，其实就是大括号内的积分部分
res=integral(houzhixishu,-1/2,1/2);
V_cr1(index)=qianzhixishu1*(res.^(-1/2));
end
semilogy(CNR1,V_cr1);hold on
end
title('V_c_r_l_b with CNR');
ylabel('V_c_r_l_b(m/s)');
xlabel('CNR(dB)');
ylim([1e-3,1e1])
grid on
参数一样的情况下，怀疑文献结果有点问题
