MATLAB常用程序集合

IIR滤波程序

function y=iirfilter(b,a,x)

nf=length(a);

c=zeros(nf-1,nf-1);

for i=nf-1:-1:1

for j=1:i

c(i,j)=b(i-j+1);

d(i,j)=-a(i-j+1);

end

end

d(1,1)=0;

y(1)=b(1)*x(1);

for i=2:nf-1

ss=0;

for j=2:i

ss=ss+x(j)*c(i,j)+y(j-1)*d(i,j-1);

end

y(i)=ss+x(1)*c(i,1);

end

a(1)=0;

for i=nf:length(x)

sum=0;

for j=2:nf

sum=sum+b(j)*x(i-j+1)-a(j)*y(i-j+1);

end

y(i)=sum+b(1)*x(i);

end

读写二进制文件

写文件

file_out = fopen('tdip.b', 'wb');

fwrite(file_out,y,'double');%int8,int16,int24,double,float

读文件

[FileName,PathName] = uigetfile('*.*','Select the M-file');%选择需要处理的数据文件%

% % 打开方式为二进制打开，其实'r'就行，matlab是默认二进制形式打开文件的file_id = fopen(FileName, 'rb');

%设置输出文件名字

jump_distance = 0; % 这里是跳过的字节数

fseek(file_id, jump_distance, 'bof');

i=0;

week=0;

week0=headt(7);

while ~feof(file_id)

% while(i<3)

[data, ele_count] = fread(file_id, [4,fs/10] , 'int32' ) ;

End

批量处理文件

编写程序时经常会遇到处理大量数据文件的情况，数据文件的命名比较相似，以一系列编号区分，如“a1.mat，a2.mat，... ，a100.mat " (假设其中的数据名称也为a1...a100) 。为了可以批量处理这些文件，可以采用如下的方法：

for i=1:1:100

s=strcat('a', int2Str(i));

load(strcat(s,'.mat'));

x=eval(s); %将数据a1赋值给x，便于后继统一处理

% ... %统一处理程序

end

写txt文件

fid=fopen(‘1.txt’,’w t’);

fprintf(fid,‘这是一个测试文件：\n’)

fprintf(fid,’%f’,a)

m=[A nskip nfft fxb real(pxb) imag(pxb)];

fid = fopen('pxb.txt', 'at');

fprintf(fid,'振幅开始时间fft点数激励频率pxb:\n');

fprintf(fid,' %d %d %d %6.8f %6.8f %6.8f\n\n', m);

% fprintf(fid,'l1 l2 l3 l4:\n')

% fprintf(fid, ' %d %d %d %d %d %d %d %d\n ', l');

fclose(fid);

%a=’00:00:00’;

% fprintf(fid,'%s ',a);

用cell生成自动变量名并幅值

clc

clear

for i=1993:1996

matName{i-1992}=strcat('Y',num2str(i)); % str=[num2str(i-1992) '.csv'];

% matName{i-1992}=xlsread(str); matName{i-1992}=rand(1,i-1992);

% eval([matName{i-1992} ' = rand(1,i);']); end

matName{1}(1)

关于errorbar

误差线：体现数据点置信水平和偏差MATLAB中errorbar的使用

example

#上下误差相等

X = 0:pi/10:pi;

Y = sin(X);

E = 0:0.05:0.5;

figure

errorbar(X,Y,E,'r');

hold on

plot(X,Y,'-mo',...

'LineWidth',2,...

'MarkerEdgeColor','k',...

'MarkerFaceColor',[.49 1 .63],... 'MarkerSize',12);

xlabel('x');ylabel('sin(x)');

#上下误差不相等

X = 0:pi/10:pi;

Y = sin(X);

L = 0:0.05:0.5;

U=0.5:-0.05:0;

figure

errorbar(X,Y,L,U,'r');

hold on

plot(X,Y,'-mo',...

'LineWidth',2,...

'MarkerEdgeColor','k',...

'MarkerFaceColor',[.49 1 .63],... 'MarkerSize',12);

xlabel('x');ylabel('sin(x)');