How to plot this function?

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

 采纳的回答

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}')

更多回答(4 个)

Take a look at erfc. A simple change of coordinates should get you there.
Then plotting can be done by
x = linspace(xmin,xmax,npoints);
y = constant*erfc(x);
plot(x,y)

2 个评论

Thank you for your answer.
But can't I put my function directly (without going through a pre-determined function)?
Thanks
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.

请先登录,再进行评论。

Michaël
Michaël 2011-3-14

0 个投票

Thank you for your answer.
But can't I put my function directly (without going through a pre-determined function)?
Thanks

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by