Matlab – zhangxudong

November 15, 2024November 15, 2024learn

根据太阳辐射来计算日照时数DLc 和根据比湿计算相对湿度RHu

% 张旭东更新20240724 % DLc——平均日照时数[hour] % Rn——冠层表面净辐射[MJ/m^2/day] % Rns——净短波辐射[MJ/m^2/day] % Rnl——净长波辐射[MJ/m^2/day] % Rs——入射的太阳辐射[MJ/m^2/day] % Ra24——碧空太阳总辐射[MJ/m^2/day] % latitude——地理纬度转换为[度] % elevation——海拔高度，[m] % RHmean——平均相对湿度，% […]

April 30, 2024April 30, 2024learn

Protected: 土壤含水率插值-采用时间序列插值

There is no excerpt because this is a protected post.

June 23, 2023June 23, 2023Uncategorized

matlab中动态变量命名及赋值

可以采用eval函数，但应记得清除eval，即clear eval，如： for i =1:9 A=random(‘Normal’,0,1,100,3); eval([‘ratio’,num2str(i),’= A’]); clear eval end

November 10, 2020November 10, 2020learn

矩阵有些行需要重复

%% 矩阵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次 […]

October 30, 2018October 30, 2018learn

用矩阵散点做等高线图和梯度场图

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); […]

December 17, 2017December 17, 2017learn

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] […]

November 19, 2017April 19, 2018learn

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 […]

November 17, 2016November 29, 2016learn

[转载]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 […]

November 17, 2016November 17, 2016learn

调用拟合好的函数

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 […]

November 3, 2016November 3, 2016learn

[转载]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; […]