c++调用matlab编程

function newbuf=improcess(filename,threshold)
imagedata=mat2gray(imread(filename));  
imagedata=rgb2gray(imagedata);
newbuf=imagedata;   
[M N]=size(imagedata);   
for i=1:1:M   
  for j=1:1:N   
      if imagedata(i,j)>threshold/255   
          newbuf(i,j)=1;   
      else   
          newbuf(i,j)=0;   
      end   
  end   
end   
imwrite(newbuf,'result.bmp');   
figure,imshow(newbuf);
return;