Matlab程序：非矩形网格图

※鱼鱼╰☆

1. function [x,y,z]=meshf(f_name,xa,xb,y_lower,y_upper,m,n)
   
2. %二元函数z=f(x,y)非矩形网格图
   
3. %用法：[x,y,z]=meshf('f_name',xa,xb,'y_lower','y_upper',m,n)
   
4. %     f_name: 表示二元函数z=f(x,y)的M函数名
   
5. %     xa,xb: x下端与上端
   
6. %     y_lower,y_upper:  y用x表示的上端函数与下端函数
   
7. %     m,n: 网格数（默认m=n=20）
   
8. %     [x,y,z]: 返回网格数据，可用mesh(x,y,z)作图，若无输出变量则自动作图
   
9. %例题 画z=exp(-x^2/2-y^2/2)/2/pi在x^2+y^2<=9内的图。
   
10. %先写三个M函数top_f.m, low_f.m, upp_f.m
    
11. %                    M函数top_f.m
    
12. %                 function z=f(x,y)
    
13. %                 z=exp(-x.^2/2-y.^2/2)/2/pi;
    
14. %                    M函数low_f.m
    
15. %                 function y=f(x)
    
16. %                 y=-sqrt(9-x.^2);
    
17. %                    M函数upp_f.m
    
18. %                 function y=f(x)
    
19. %                 y=sqrt(9-x.^2);
    
20. %现在用命令
    
21. %     meshf('top_f',-3,3,'low_f','upp_f')
    
22. %
    
23. % Purpose: Generate mesh of z=f(x,y)
    
24. % Synopsis: [x,y,z]=meshf('f_name',xa,xb,'y_lower','y_upper',m,n)
    
25. %                f_name: file name of 2-D function f(x,y)
    
26. %                xa,xb: lower and upper limit of x
    
27. %                y_lower,y_upper: file names of lower and upper
    
28. %                                limit of y (1-D function of x)
    
29. %                m,n: number of meshgrid (Default: m=n=20)
    
30. %                [x,y,z]: return mesh data used by MESH of SURF
    
31. 
    
32. % L.J.Hu 8-8-1998
    
33. 
    
34. if nargin<7,n=20;end
    
35. if nargin<6,m=20;end
    
36. 
    
37. x=linspace(xa,xb,m);
    
38. yc=feval(y_lower,x);yd=feval(y_upper,x);
    
39. 
    
40. x=ones(n,1)*x;
    
41. y=zeros(n,m);
    
42. for i=1:m
    
43.   y(:,i)=linspace(yc(i),yd(i),n)';
    
44. end
    
45. z=feval(f_name,x,y);
    
46. 
    
47. if nargout==0,mesh(x,y,z);xlabel('X');ylabel('Y');end
    

复制代码