Matlab:绘图误差分析

Matlab:绘图误差分析,matlab,error-handling,plot,signal-processing,Matlab,Error Handling,Plot,Signal Processing,我创建了一个函数,用高斯噪声表示浓度与时间的理想特性和浓度与时间的真实特性 下面是理想函数和实函数的曲线 我想绘制误差曲线,表示理想值和实际值之间的差异(带噪声) 特性必须以(噪声-X轴与误差-Y轴)表示,而不是以浓度与时间表示 下面是生成上述两条曲线的函数 function [c_t,c_t_noise] =Noise_ConstrainedK2(t,a1,a2,a3,b1,b2,b3,td,tmax,k1,k2,k3) K_1 = (k1*k2)/(k2+k3); K_2

我创建了一个函数,用高斯噪声表示浓度与时间的理想特性和浓度与时间的真实特性

下面是理想函数和实函数的曲线

我想绘制误差曲线,表示理想值和实际值之间的差异(带噪声)

特性必须以(噪声-X轴与误差-Y轴)表示,而不是以浓度与时间表示

下面是生成上述两条曲线的函数

function [c_t,c_t_noise] =Noise_ConstrainedK2(t,a1,a2,a3,b1,b2,b3,td,tmax,k1,k2,k3)



K_1   = (k1*k2)/(k2+k3);
K_2   = (k1*k3)/(k2+k3);
%DV_free= k1/(k2+k3);




c_t = zeros(size(t));

ind = (t > td) & (t < tmax);

  c_t(ind)= conv(((t(ind) - td) ./ (tmax - td) * (a1 + a2 + a3)),(K_1*exp(-(k2+k3)*t(ind)+K_2)),'same');



ind = (t >= tmax);

  c_t(ind)=conv((a1 * exp(-b1 * (t(ind) - tmax))+ a2 * exp(-b2 * (t(ind) - tmax))) + a3 * exp(-b3 * (t(ind) - tmax)),(K_1*exp(-(k2+k3)*t(ind)+K_2)),'same');

for mean=1:10
 for sigma=1:45

 c_t_noise = c_t+normrnd(mean,sigma,[1 500],length(c_t));

 end

end
dg = [0 0.5 0];
plot(t,c_t,'r');
hold on; 
plot(t,c_t_noise,'Color',dg);
hold off; 

axis([0 50 0 1900]);
xlabel('Time[mins]');
ylabel('concentration [Mbq]');
title('My signal');

end
函数[c_t,c_t_noise]=噪声受约束k2(t,a1,a2,a3,b1,b2,b3,td,tmax,k1,k2,k3)
K_1=(k1*k2)/(k2+k3);
K_2=(k1*k3)/(k2+k3);
%DV_free=k1/(k2+k3);
c_t=零(大小(t));
ind=(t>td)和(t=tmax);
c_t(ind)=conv((a1*exp(-b1*(t(ind)-tmax))+a2*exp(-b2*(t(ind)-tmax))+a3*exp(-b3*(t(ind)-tmax)),(K_1*exp(-k2+k3)*t(ind)+K_2)),‘相同’;
平均值=1:10
对于sigma=1:45
c_t_噪声=c_t+正常值(平均值,西格玛[1500],长度(c_t));
结束
结束
dg=[0.50];
图(t,c_t,'r');
等等
绘图(t、c、t、噪音、‘颜色’、dg);
拖延;
轴([0 50 0 1900]);
xlabel('时间[分钟]);
ylabel('浓度[Mbq]');
标题(“我的信号”);
结束

请建议我一种方法,我可以计算c_t和c_t_噪声之间的差异,并绘制它

-1:plot(t,(c_t-c_t_噪声))?@Schorsch我想用(噪声(变量Sigma)与误差)而不是浓度与时间来表示图,噪声与误差之间的差异是什么?
c\u t\u noise-c\u t
等于噪声还是误差?