fourier這個函數在后面就取消不用了,只有老版本的matlab還能使用,當然也不是不能用，現在的FOURIER函數要求輸入是sym型才能用，help fourier會顯示未定義的函數，需要help sym/fourier才能看到用法，現在多用fft函數吧。

求極限命令limit是對于有變量的算式操作的，對于數值是不能求極限的，如limit(sin(x),x,0)
結果是：ans = 0。
但是，limit(10)就沒意義，不能操作，會有錯誤提示。

function [D_B,D_C]=extraction(txt,num)%[num,txt,raw]=xlsread('D:\HeRui\selfd_fun\datafortest.xlsx');%txt為原始時間數據；格式為cell;%C為與txt對應的數據B=char(txt);D_B=B;D_C=num;for i=1:1:320if double(B(i,6))<57D_B(i,:)='0';end;for j=0:1:320-iif double(B(i+j,6))==57D_B(i+j,:)=B(i+j,:);end;if double(B(i+j,6))==49&&double(B(i+j,7))==48D_B(i+j,:)=B(i+j,:);endend;if double(B(i,7))>48D_B(i)='0';end;end;D_B

A=imread('50粒 小底盤.JPG');
B=rgb2gray(A);
figure;imshow(B);title('原始圖像');
B= double(B);
H=imnoise(B,'gaussian');

Matlab中沒有現成的Lagrange插值函數，必須編寫一個M文件實現Lagrange插值。
編寫一個名為lagrange.m的M文件：

function y=lagrange(x0,y0,x);
n=length(x0);m=length(x);
for i=1:m
z=x(i);
s=0.0;
for k=1:n
p=1.0;
for j=1:n
if j~=k
p=p*(z-x0(j))/(x0(k)-x0(j));
end
end
s=p*y0(k)+s;
end
y(i)=s;
end