函数绘图---例题

森之张卫东

创建一个函数的函数，它能够画出所有只有一个自变量的MATLAB函数的图象，自变量的范围是用户指定的始值和终值。

答案：

这个函数有两个输入参数，第一个是要画的函数的函数名，第二个是两元素向量，它指明了画图的取值范围。

1.陈述问题

创建一个函数的函数，它能够画出所有只有一个自变量的MATLAB函数的图象，自变量的范围由用户指定。

2.定义输入输出函数的输入有两个

(1)包含有函数名的字符串

(2)包含有起始值和终值的2元素向量函数的输出是所要画的图象

3.设计算法这个函数可以分为4大步：

Check for a legal number of arguments

Check that the second argument has two elements

Calculate the value of the function between the start and stoppoints

Plot and label the function

第三四大步的伪代码如下

n_steps ← 100

step_size ← (xlim(2) – xlim(1)) / nsteps

x ← xlim(1):step_size:xlim(2)

y ← feval(fun, x)

plot(x, y)

title(['bf Plot of function ' fun ' (x)'])

xlabel('\bfx;)

ylabel(['bf ' fun ' (x)'])

森之张卫东

2. function quickplot(fun,xlim)
   
3. %QUICKPLOT Generate quick plot of a function
   
4. % Function QUICKPLOT generates a quick plot
   
5. % of a function contained in a external mfile,
   
6. % between user-specified x limits.
   
7. % Define variables:
   
8. % fun                   --Function to plot
   
9. % msg                   --Error message
   
10. % n_steps               --Number of steps to plot
    
11. % step_size             --Step size
    
12. % x                     --X-values to plot
    
13. % y                     --Y-values to plot
    
14. % xlim                  --Plot x limits
    
15. % Record of revisions:
    
16. % Date Programmer Description of change
    
17. % ==== ========== =====================
    
18. % 12/17/98 S. J. Chapman Original code
    
19. % Check for a legal number of input arguments.
    
20. msg = nargchk(2,2,nargin);
    
21. error(msg);
    
22. % Check the second argument to see if it has two
    
23. % elements. Note that this double test allows the
    
24. % argument to be either a row or a column vector.
    
25. if ( size(xlim,1) == 1 & size(xlim,2) == 2 ) | ...
    
26. ( size(xlim,1) == 2 & size(xlim,2) == 1 )
    
27.     % Ok        --continue processing.
    
28.     n_steps = 100;
    
29.     step_size = (xlim(2) - xlim(1)) / n_steps;
    
30.     x = xlim(1):step_size:xlim(2);
    
31.     y = feval(fun,x);
    
32.     plot(x,y);
    
33.     title(['\bfPlot of function ' fun '(x)']);
    
34.     xlabel('\bfx');
    
35.     ylabel(['\bf' fun '(x)']);
    
36. else
    
37.     % Else wrong number of elements in xlim.
    
38.     error('Incorrect number of elements in xlim.');
    
39. end
    

function quickplot(fun,xlim)<br /> %QUICKPLOT Generate quick plot of a function<br /> % Function QUICKPLOT generates a quick plot<br /> % of a function contained in a external mfile,<br /> % between user-specified x limits.<br /> % Define variables:<br /> % fun                   --Function to plot<br /> % msg                   --Error message<br /> % n_steps               --Number of steps to plot<br /> % step_size             --Step size<br /> % x                     --X-values to plot<br /> % y                     --Y-values to plot<br /> % xlim                  --Plot x limits<br /> % Record of revisions:<br /> % Date Programmer Description of change<br /> % ==== ========== =====================<br /> % 12/17/98 S. J. Chapman Original code<br /> % Check for a legal number of input arguments.<br /> msg = nargchk(2,2,nargin);<br /> error(msg);<br /> % Check the second argument to see if it has two<br /> % elements. Note that this double test allows the<br /> % argument to be either a row or a column vector.<br /> if ( size(xlim,1) == 1 & size(xlim,2) == 2 ) | ...<br /> ( size(xlim,1) == 2 & size(xlim,2) == 1 )<br />     % Ok        --continue processing.<br />     n_steps = 100;<br />     step_size = (xlim(2) - xlim(1)) / n_steps;<br />     x = xlim(1):step_size:xlim(2);<br />     y = feval(fun,x);<br />     plot(x,y);<br />     title(['\bfPlot of function ' fun '(x)']);<br />     xlabel('\bfx');<br />     ylabel(['\bf' fun '(x)']);<br /> else<br />     % Else wrong number of elements in xlim.<br />     error('Incorrect number of elements in xlim.');<br /> end<br />

森之张卫东

附件：详细内容。