算是matlab的一些积累叭（常用基本操作）

绘图

基本操作

%% 用来绘图
plot(L(1:18),xp)                            %输出经验频率曲线
xlabel('P(%)')                              %X轴名  P(%)
ylabel('Q(m^{3}/s)')                        %Y轴名  降水mm)
%set(gca,'ytick',[0:10:1400])
set(gca,'yMinorTick','on')
set(gca,'XTick',L(1:18) )                   %X轴按L()算法伸缩            
set(gca,'XTickLabel',p(1:18)*100)           %按p（1：18）标识伸缩后的X轴
grid on                                     %加网格线  
hold on                                     %保持图形    
title(strcat('Cv=',num2str(cv),';Cs/Cv=',num2str(k)))
plot(L(19:(18+n)),w,'*k')                   %绘制经验点 黑色的*
hold off                                    %取消图形保持  

颜色设置

渐变色散点图

% 绘制绳套曲线
ColorMaker=1:1:n;
plot(Q_e,W);
hold on scatter(Q_e,W,10,ColorMaker,'o','filled'); xlabel('Q/(m^3/s)'); ylabel('W(12h.m^3/s)');

坐标轴颜色设置

ylabel('up stream','color','b');%表示 up stram是坐标轴名称，color和b表示蓝色

循环语句

while

用途

依据某一判断一直重复进行，直到满足条件后再退出
用 ctrl+C可以强制退出

对外交互（？）

input

mm=input('是否满意？（1\\2）');
if mm==1    
   break
else
   disp('好嘞，继续运行... ...')
   continue

注意：

1. 换行符在字符串里用 \

读写文件

读取文件

filename = "D:\1-3大四上\2-水文预报\sj\上机实验2\马斯京根法.xls"; data=xlsread(filename,'A7:d24');

写出文件

1. 判断文件是否存在，并且输出，下段代码摘自某各函数定义中，该代码要求输入的文件名是字符串！！！！并且保存位置是在该运行文件夹下。

filename3=Rfile;
fid=fopen(filename3,'wt');
if fid==-1
   error('Error when opening training set file !');                         
end
fprintf(fid,'%f\n',R);                                      
fclose(fid);

常见函数

拟合函数

1. ployval（代入）&polyfit（求解）

r=polyfit(Q_e,W,1); K=r(1);									%蓄泄方程
x=xx(m);
%直线拟合R2
W=polyval(r,QQ); 
R2=1-sum((QQ-Q_e).^2)/sum((Q_e-mean(Q_e)).^2);

寻找极值点的函数(findpeaks)

%% 设置寻找极值点的函数
function [maxx,loc_max]=peak(s,mph) [maxx,loc_max]=findpeaks(s,'minpeakheight',mph);
 end

合并字符串（strcat）

 title(strcat('x=',num2str(x),'K=',num2str(round(r(1),3)),'误差max',num2str(round(delt_Q_max,3))));