Protected: 土壤含水率插值-采用时间序列插值
There is no excerpt because this is a protected post.
Tag: Matlab
matlab中动态变量命名及赋值
可以采用eval函数，但应记得清除eval，即clear eval，如： for i =1:9 A=random(‘Normal’,0,1,100,3); eval([‘ratio’,num2str(i),’= A’]); clear eval end
矩阵有些行需要重复
%% 矩阵a每一行重复n次 a为已知矩阵 m=size(a,1); % a有多少行？ b=[] n=5; for i=1:m b0=repmat(a(i,:),n,1); b=[b;b0]; end %% 矩阵a每m行重复n次，a为已知矩阵 b=[]; n=5; m=5; k=size(a,1)/m; %循环k次 […]
用矩阵散点做等高线图和梯度场图
matlab提供了可以根据矩阵做等高线的命令contour，和求梯度的函数gradient，可以根据梯度用quiver命令将梯度做成箭头图。但要注意一些参数的设置。 一般情况下，例如： x=-2:0.2:2; y=-2:0.2:2; [X,Y]=meshgrid(x,y); Z=X.^3+Y.^2; contour(X,Y,Z); [DX,DY]=gradient(Z,0.2,0.2); hold on; quiver(X,Y,DX,DY); 作图如下： 设置一下不同的坐标，纵坐标的范围为 [-20 20]，间距为2。 x=-2:0.2:2; y=-20:2:20; [X,Y] = meshgrid(x,y); […]
Matlab将元胞数组输出到txt或者dat文件中
原链接：http://cn.mathworks.com/help/matlab/import_export/write-to-delimited-data-files.html#br2ypq2-1 方法一：用fprint输出元胞数组 这个例子展示了怎么用fprint方法将元胞数组输出到TEXT 文件中。 创建一个元胞数组例子, 以便输出 C = {‘Atkins’,32,77.3,’M’;’Cheng’,30,99.8,’F’;’Lam’,31,80.2,’M’} C = 3×4 cell array ‘Atkins’ [32] [77.3000] ‘M’ ‘Cheng’ [30] […]
MATLABR2012b破解lisence过期的问题（转载）
昨天一早就被matlab无法正常使用搞得天昏地暗，重装卸载重装了很多次，还是反复出现需要激活的页面，弄了大半天最后发现原来是matlab的lisence文件的截止时间是2017年11月。修改了lisence文件后立即便可用，快被自己蠢哭了。 附上lisence文件代码： %%%%%%%%%%%%%%%%%%%%%%%%%%%% INCREMENT Aerospace_Blockset MLM 99 permanent uncounted \ A05070F00D1EB1F92326 VENDOR_STRING=QQ=47399897 HOSTID=ANY \ ck=216 SN=888888 TS_OK INCREMENT Aerospace_Toolbox MLM […]
[转载]Matlab拟合工具箱的几个误差参数说明
Matlab的拟合、优化和统计等工具箱常有下面几个名词： SSE(和方差、误差平方和)：The sum of squares due to error MSE(均方差、方差)：Mean squared error RMSE(均方根、标准差)：Root mean squared error R-square(确定系数)：Coefficient of determination Adjusted R-square：Degree-of-freedom […]
调用拟合好的函数
Matlab 2016要调用一个拟合好的函数非常简单： f = fittype(‘a*x^2+b*exp(n*x)’) %给出函数 f = General model: f(a,b,n,x) = a*x^2+b*exp(n*x) c = cfit(f,1,10.3,-1e2) %给出相应参数，返回函数值 c […]
[转载]matlab化简符号表达式
化简符号表达式计算机毕竟还是挺笨的, 经过一系列的符号计算后, 得到的结果可能只有它自己才能看懂， Matlab提供大量函数以用于符号表达式的化简. collect(f): 函数用途是合并多项式中相同的项, 如:syms x t f=(1+x)*t+x*t; collect(f) expand(f):展开多项式,syms x f=x*(x*(x-1)+3)+2; expand(f); horner(f) 对转换多项式为Horner形式, 这种形式的特点是乘法嵌套, 其有着不错的数值计算性质.syms x; […]
长短矩阵的对齐
短矩阵向长矩阵对齐 A=[]; % 长系列序号 B=[]; % 短系列序号 Z=cell(1,1); % 和短系列同长度的矩阵，可包含汉字 [Lia, Locb]=ismember(A,B,’rows’); % 把短系列的矩阵向长系列看齐，缺的值以[]填写 for i=1:749; if Locb(i)==0 […]