fourier這個(gè)函數(shù)在后面就取消不用了,只有老版本的matlab還能使用,當(dāng)然也不是不能用，現(xiàn)在的FOURIER函數(shù)要求輸入是sym型才能用，help fourier會(huì)顯示未定義的函數(shù)，需要help sym/fourier才能看到用法，現(xiàn)在多用fft函數(shù)吧。

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

function [D_B,D_C]=extraction(txt,num)%[num,txt,raw]=xlsread('D:\HeRui\selfd_fun\datafortest.xlsx');%txt為原始時(shí)間數(shù)據(jù)；格式為cell;%C為與txt對(duì)應(yīng)的數(shù)據(jù)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中沒有現(xiàn)成的Lagrange插值函數(shù)，必須編寫一個(gè)M文件實(shí)現(xiàn)Lagrange插值。
編寫一個(gè)名為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