分析此代码运行结果:f= @(x) sin(x); % 定义函数f(x) t= @(i) i.^2; % 定义函数t(i) x1 = 0:1/100:2*pi; %生成离散点 i1 = 0:1/100:1; %生成离散点 f1=f(x1); t1=t(i1); sin_mean=mean(f1) sin_var=var(f1) e_mean=mean(t1) e_var=var(t1) noise1=randn(size(x1)); noise2=randn(size(i1)); F_Noise = f(x1) + rand(size(x1)).*f(x1);%添加均匀分布噪声 sin_rand_mean=mean(F_Noise) sin_rand_var=var(F_Noise) F_Noise1 = f(x1) + randn(size(x1)).*f(x1);%添加正态分布噪声 sin_randn_mean=mean(F_Noise1) sin_randn_var=var(F_Noise1) t_Noise = t(i1) + (2*rand(size(i1))-1).*t(i1);%添加均匀分布噪声 r_rand_mean=mean(t_Noise) r_rand_var=var(t_Noise) t_Noise1= t(i1) + randn(size(i1)).*t(i1);%添加正态分布噪声 r_randn_mean=mean(t_Noise1) r_randn_var=var(t_Noise1) a=[1]; %分子的系数 b=[2,1]; %分母的系数 sys=tf(a,b); %生成 RC 系统的传递函数 k=lsim(sys,F_Noise1,x1); %求出系统在特定输入的情况下的输出 k1=lsim(sys,t_Noise1,i1); figure(1); subplot(421);plot(x1,f1);title('正弦函数'); subplot(422);plot(x1,noise1);title('白噪声信号'); subplot(423);plot(x1,F_Noise);title('加上相对误差水平为1的均匀分布噪声的正弦函数') subplot(424);plot(x1,F_Noise1);title('加上相对误差水平为1的正态分布噪声的正弦函数') subplot(425);autocorr(f1);xlabel('滞后阶数lag');ylabel('数据值');title('正弦信号的自相关函数ACF图') subplot(426);autocorr(noise1);xlabel('滞后阶数lag');ylabel('数据值');title('白噪声信号的自相关函数ACF图') subplot(427);plot(k);xlabel('输出样本数');ylabel('幅度'); title('加正态分布白噪的正弦信号通过rc系统后输出信号')
时间: 2023-08-14 15:03:50 浏览: 406
sinc:函数 sinc(x)=sin(pi*x)/(pi*x)-matlab开发
此代码主要是对正弦函数和平方函数添加各种类型的噪声进行分析,并对 RC 系统进行模拟求解。具体运行结果如下:
1. 定义了函数 $f(x) = \sin(x)$ 和 $t(i)=i^2$,并生成了离散点 $x_1$ 和 $i_1$。
2. 分别计算了 $f(x_1)$ 和 $t(i_1)$ 的均值和方差,分别为 $sin\_mean$、$sin\_var$、$e\_mean$ 和 $e\_var$。
3. 生成了两个相同大小的白噪声信号 $noise1$ 和 $noise2$。
4. 对 $f(x_1)$ 和 $t(i_1)$ 添加了相对误差水平为 $1$ 的均匀分布噪声和正态分布噪声,分别保存在 $F\_Noise$、$F\_Noise1$、$t\_Noise$ 和 $t\_Noise1$ 中。
5. 计算了添加噪声后的 $F\_Noise$ 和 $F\_Noise1$ 的均值和方差,分别为 $sin\_rand\_mean$、$sin\_rand\_var$ 和 $sin\_randn\_mean$、$sin\_randn\_var$。
6. 计算了添加噪声后的 $t\_Noise$ 和 $t\_Noise1$ 的均值和方差,分别为 $r\_rand\_mean$、$r\_rand\_var$ 和 $r\_randn\_mean$、$r\_randn\_var$。
7. 定义了 RC 系统的传递函数 $sys=tf(a,b)$,其中 $a$ 和 $b$ 分别为分子和分母的系数。
8. 使用 lsim 函数模拟了输入信号 $F\_Noise1$ 和 $t\_Noise1$ 在 RC 系统下的输出信号,分别保存在 $k$ 和 $k1$ 中。
9. 绘制了多个子图,包括正弦函数、白噪声信号、添加均匀分布噪声和正态分布噪声后的正弦函数、正弦信号的自相关函数 ACF 图、白噪声信号的自相关函数 ACF 图以及通过 RC 系统后的正弦信号输出信号。
阅读全文