帮我看一个程序，给这个程序求出的矩阵赋予个名称

xyq001303

function imgrad=shiyan1a=imread('klj.jpg');imshow(a);[x y]=ginput(2);x=round(x);y=round(y);n1=abs(x(1)-x(2))+1;n2=abs(y(1)-y(2))+1;[n which_max]=max([n1,n1]);imgrad=ones(3,n);for i=1:3    b=a(:,:,i);    if y(1)<y(2)&&x(1)<x(2)        if which_max==1            xx=x(1):x(2);            yy=subf(x,y,xx);        else            yy=y(1):y(2);            xx=subf(x,y,yy);        end    elseif y(1)<y(2)&&x(1)>x(2)        if which_max==1            xx=x(2):x(1);            yy=subf(x,y,xx);        else            yy=y(1):y(2);            xx=subf(x,y,yy);        end    elseif y(1)>y(2)&&x(1)<x(2)        if which_max==1            xx=x(1):x(2);            yy=subf(x,y,xx);        else            yy=y(2):y(1);            xx=subf(x,y,yy);        end    elseif y(1)>y(2)&&x(1)>x(2)        if which_max==1            xx=x(2):x(1);            yy=subf(x,y,xx);        else            yy=y(2):y(1);            xx=subf(x,y,yy);        end    end        ind=sub2ind(size(b),yy,xx);    data=b(ind);    data=double(data);    imgrad(i,:)=gradient(data);endfunction coordinate=subf(x,y,input)k=(y(1)-y(2))/(x(1)-x(2));if abs(x(1)-x(2))>=abs(y(1)-y(2))    coordinate=k*(input-x(1))+y(1);    coordinate=round(coordinate);else    coordinate=(input-y(1))/k+x(1);    coordinate=round(coordinate);end 以上代码是用来给图片求颜色变化梯度的，求出来的是一个3行n列的矩阵，我想给这个矩阵赋予一个名称，方便后面对其进行操作，但是不知道该怎么添加代码，请高手解惑，谢谢！！！

messenger

真服了，怎么贴出来的代码没有换行呢？

如果想给求出来的矩阵赋一个名，只要在调用函数前面加=号就可以了。

比如，你将这段程序存为shiyan.m文件，在调用时用 x=shiyan，就把求出来的矩阵赋成 x 这个矩阵名了。


建议你找一本基础的Matlab书，看看函数调用方面的内容。

1. 
   
2. 
   
3. 
   
4. function imgrad=shiyan1
   
5. a=imread('klj.jpg');
   
6. imshow(a);
   
7. [x y]=ginput(2);
   
8. x=round(x);
   
9. y=round(y);
   
10. n1=abs(x(1)-x(2))+1;
    
11. n2=abs(y(1)-y(2))+1;
    
12. [n which_max]=max([n1,n1]);
    
13. imgrad=ones(3,n);
    
14. for i=1:3   
    
15. b=a(:,:,i);   
    
16. if y(1)<y(2)&&x(1)<x(2)        
    
17. if which_max==1           
    
18. xx=x(1):x(2);            
    
19. yy=subf(x,y,xx);        
    
20. else            
    
21. yy=y(1):
    
22. y(2);           
    
23. xx=subf(x,y,yy);        
    
24. end   
    
25. elseif
    
26. y(1)<y(2)&&x(1)>x(2)        
    
27. if which_max==1            
    
28. xx=x(2):x(1);            
    
29. yy=subf(x,y,xx);        
    
30. else            
    
31. yy=y(1):y(2);            
    
32. xx=subf(x,y,yy);        
    
33. end   
    
34. elseif y(1)>y(2)&&x(1)<x(2)        
    
35. if which_max==1            
    
36. xx=x(1):x(2);            
    
37. yy=subf(x,y,xx);        
    
38. else            
    
39. yy=y(2):y(1);            
    
40. xx=subf(x,y,yy);        
    
41. end   
    
42. elseif y(1)>y(2)&&x(1)>x(2)        
    
43. if which_max==1            
    
44. xx=x(2):x(1);            
    
45. yy=subf(x,y,xx);        
    
46. else            
    
47. yy=y(2):y(1);            
    
48. xx=subf(x,y,yy);        
    
49. end   
    
50. end        
    
51. ind=sub2ind(size(b),yy,xx);   
    
52. data=b(ind);   
    
53. data=double(data);   
    
54. imgrad(i,:)=gradient(data);
    
55. end
    
56. 
    
57. function coordinate=subf(x,y,input)
    
58. k=(y(1)-y(2))/(x(1)-x(2));
    
59. if abs(x(1)-x(2))>=abs(y(1)-y(2))   
    
60. coordinate=k*(input-x(1))+y(1);   
    
61. coordinate=round(coordinate);
    
62. else   
    
63. coordinate=(input-y(1))/k+x(1);   
    
64. coordinate=round(coordinate);
    
65. end
    
66. 
    

复制代码