How to plot this function?
1 view (last 30 days)
Show older comments
Hello,
Could you please help me plot this function with Matlab ?
<http://img638.imageshack.us/img638/3205/fonction.png >
I'm a beginner. Your help will be valuable to me.
Many thanks
0 Comments
Accepted Answer
Matt Fig
on 15 Mar 2011
If you want to do it manually, try this:
f = @(z) quadv(@(x) (1/(3*sqrt(2*pi)))*exp(-(x-30).^2/18),z,100);
Z = 0:.25:50;
G = zeros(1,length(Z));
for ii = 1:length(Z)
G(ii) = f(Z(ii));
end
plot(Z,G)
xlabel('z')
ylabel('(1/\surd(2*pi))\inte^{-(x-30)^2/18}')
0 Comments
More Answers (4)
Matt Tearle
on 14 Mar 2011
Then plotting can be done by
x = linspace(xmin,xmax,npoints);
y = constant*erfc(x);
plot(x,y)
2 Comments
Matt Tearle
on 15 Mar 2011
Why would you want to? This function is defined via an improper integral, which means you'll have to do an infinite-domain integral for each value of z. This will be slow and will introduce all sorts of numerical errors. On the other hand, the error function is a well-known function -- its implementation in MATLAB is stable and tested. Why reinvent the wheel?
But if you really want to, you can use something like quad to do the integral.
See Also
Categories
Find more on Multirate Signal Processing in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!