使用matlab代码生成一个中心频率为10KHZ,带宽为400HZ方差为1的窄带高斯随机过程，与幅度为2相位为30°的余弦信号合成信号的包络，相位，包络平方的概率分布

好的，我来回答你的问题。

首先，我们可以使用MATLAB中的`randn`函数生成一个方差为1的高斯随机过程，然后使用`bandpass`函数将其限制在10KHz中心频率和400Hz带宽内。具体代码如下：

fs = 100e3; % 采样率
fc = 10e3; % 中心频率
bw = 400; % 带宽
var = 1; % 方差

t = 0:1/fs:1; % 时间范围
x = sqrt(var)*randn(size(t)); % 方差为1的高斯随机过程
y = bandpass(x, [fc-bw/2, fc+bw/2], fs); % 限制在10KHz中心频率和400Hz带宽内

接下来，我们可以生成幅度为2、相位为30°的余弦信号，并将其与随机过程合成：

A = 2; % 幅度
phi = 30/180*pi; % 相位

z = A*cos(2*pi*fc*t + phi); % 余弦信号
s = y.*z; % 合成信号

然后，我们可以使用MATLAB的`hilbert`函数将信号转换为解析信号，进而计算其包络、相位和包络平方：

s_analytic = hilbert(s); % 转换为解析信号
envelope = abs(s_analytic); % 计算包络
phase = angle(s_analytic); % 计算相位
envelope_square = abs(s_analytic).^2; % 计算包络平方

最后，我们可以使用MATLAB的`histogram`函数绘制包络、相位和包络平方的概率分布图：

figure;
subplot(3,1,1);
histogram(envelope);
title('Envelope Probability Distribution');
subplot(3,1,2);
histogram(phase);
title('Phase Probability Distribution');
subplot(3,1,3);
histogram(envelope_square);
title('Envelope Square Probability Distribution');

这样就可以得到包络、相位和包络平方的概率分布图了。