Matlab: Plotting error Analysis -

-

i have created function expresses ideal characteristics of concentration vs time , real characteristics of concentration vs time gaussian noise.

following curves of both ideal , real function

ideal real

i plot error curve expresses difference between ideal , real values(with noise)

the characteristics must expressed in( noise - x axis vs error -y axis) instead of concentration vs time

following function generates 2 curves expressed above

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');  mean=1:10  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

kindly suggest me way in can possibly compute difference between c_t , c_t_noise , plot it


Comments

Post a Comment

Popular posts from this blog

windows - Single EXE to Install Python Standalone Executable for Easy Distribution -

-
i used pyinstaller create standalone exe file pyqt project. did not use single file option because exe created single file option takes long open. therefore, pyinstaller created dist folder contains bunch of files including program.exe file. create setup.exe file place contents of dist folder in program directory , create shortcut in start menu , on desktop. want super simple user. maybe setup.exe files when download program cnet. found inno-setup, looks promising. however, not know if there special considerations because program standalone program create python/pyqt program. thanks! have experience task? there program using task not know about? inno-setup or nsis easiest use. tell them files include , put them , create nice installer you. wrote short tutorial on experiences using innosetup might find helpful: http://www.blog.pythonlibrary.org/2008/08/27/packaging-wxpymail-for-distribution/ note tutorial based around wxpython app wrote, concepts same.
Read more

c# - Access objects in UserControl from MainWindow in WPF -

-
i have simple xaml, how can change text property of textblock in usercontrol1 mainwindow ? <window x:class="refactorxaml.mainwindow" xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation" xmlns:x="http://schemas.microsoft.com/winfx/2006/xaml" xmlns:refactorxaml="clr-namespace:refactorxaml" title="mainwindow" height="350" width="525"> <grid> <grid.rowdefinitions> <rowdefinition height="100"></rowdefinition> <rowdefinition height="100"></rowdefinition> </grid.rowdefinitions> <grid> <refactorxaml:usercontrol1></refactorxaml:usercontrol1> </grid> <grid grid.row="1"> <button>change text</button> </grid> </grid> </window> <usercon...
Read more

javascript - How to name a jQuery function to make a browser's back button work? -

-
i'm trying write function ajaxyfy web site, works perfect, except function button. that's have now: function loading() { if (typeof history.pushstate !== "undefined") { var historycount = 0; $('.menu-item, .logo').on('click','a',function(e){ e.preventdefault(); $(this).data('clicked', true); if ($(this).is('.logo a')) { var homepage = $(this).attr('href'); function1(homepage); history.pushstate(null, null, homepage); } else if ($(this).is('.projects a')) { var projects = $(this).attr('href'); function2(projects); history.pushstate(null, null, projects); } else { var pages = $(this).attr('href'); function3(pages); history.pushstate(null, null, pages); ...
Read more