matlab常用代码

1.将数组数据存到text文件中

fid=fopen('test1.txt','wt'); %写的方式打开文件（若不存在，建立文件）；
fprintf(fid,'%d ',Scan_Pha);  % %d 表示以整数形式写入数据，这正是我想要的；
fclose(fid);  %关闭文件；

2.将txt文件数据导入matlab中成为数组

aa = textread('Pha100to20kampzero1.txt');

3.由数据画bode图

figure;
color1=[1,0,0]; %[0.64706	0.16471	0.16471];
color2=[0,0.45,0.74];%[0	0 1];
color4=[1,0,1];
color3=[0,0,0];%[0.23529 0.70196 0.44314];

subplot(2,1,1);
% semilogx(bodepsim1(:,1),20*log10(bodepsim1(:,2)),'-','Color',color1,'LineWidth',1.2);  
semilogx(Scan_Fre,20*log10(Scan_Amp),'-','Color',color1,'LineWidth',1.2);  
hold on;
% set(gca,'XLim',[1 1000],'Fontname', 'Helvetica','FontSize',8);
% set(gca,'YLim',[20 100],'Fontname', 'Helvetica','FontSize',8);
% ylabel('Magnitude(dB)','Fontname', 'Helvetica','FontSize',10);
grid on

4.指数矩阵数值积分法

当碰到矩阵无法求逆时，不能求出其解析解

A=rand(3);
fun=@(x)expm(A*x);
Q = integral(fun,0,1,'ArrayValued',true);