ATL

　　1、std函数是用来计算标准偏差的一个函数，由于其有不同的参数，我们就用下面的例子进行介绍：　　A =　　1 2 3　　1 1 1　　2、std(A):　　std（A）函数求解的是最常见的标准差，此时除以的是N-1。　　注意：此函数命令不能对矩阵求整体的标准差，只能按照行或者列进行逐个求解标准差，默认情况下是按照列。　　在MATLAB主窗口中输入std(A) 回车，结果如下：　　输出的是每一列的标准差。　　3、std(A，flag):　　这里flag代表的是用哪一个标准差函数，如果取0，则代表除以N-1，如果是1代表的是除以N，在MATLAB主窗口中输入std(A，1) 回车，std(A，0) 回车，可以看到如下结果　　4、std(A，flag，dim):　　第三个参数代表的是按照列求标准差还是按照行求标准差，std(A,1,1)代表的是按照列求标准差，std(A,1,2)代表的是按照行求标准差。　　在MATLAB主窗口中输入如下命令：std(A,1,1) 敲回车 std(A,1,2) 敲回车。

1、std函数是用来计算标准偏差的一个函数，由于其有不同的参数，我们就用下面的例子进行介绍：　　A=　　123　　111　　2、std(A):　　std（A）函数求解的是最常见的标准差，此时除以的是N-1。　　注意：此函数命令不能对矩阵求整体的标准差，只能按照行或者列进行逐个求解标准差，默认情况下是按照列。　　在MATLAB主窗口中输入std(A)回车，结果如下：　　输出的是每一列的标准差。　　3、std(A，flag):　　这里flag代表的是用哪一个标准差函数，如果取0，则代表除以N-1，如果是1代表的是除以N，在MATLAB主窗口中输入std(A，1)回车，std(A，0)回车，可以看到如下结果　　4、std(A，flag，dim):　　第三个参数代表的是按照列求标准差还是按照行求标准差，std(A,1,1)代表的是按照列求标准差，std(A,1,2)代表的是按照行求标准差。　　在MATLAB主窗口中输入如下命令：std(A,1,1)敲回车std(A,1,2)敲回车。

STD是C++语言命令。std其实就是standard标准的意思。例如std::cin就是标准输入std::cout就是标准输出的意思~:)一个标准库名字空间The std namespaceThe ANSI/ISO C++ standard requires you to explicitly declare the namespace in the standard library. For example,when using iostream.h,you do not have to specify the namespace of cout in one of the following ways:std::cout (explicitly)using std::cout (using declaration)using namespace std (using directive)Visual C++ continues to support the use of older header filenames with the .h extension. Such usage does not require invocation of the std namespace described above.扩展资料std（A）函数求解的是最常见的标准差，此时除以的是N-1。注意：此函数命令不能对矩阵求整体的标准差，只能按照行或者列进行逐个求解标准差，默认情况下是按照列。在MATLAB主窗口中输入std(A) 回车，结果如下：输出的是每一列的标准差。在MATLAB主窗口中输入如下命令：std(A,1,1) 敲回车 std(A,1,2) 敲回车。参考资料来源：百度百科-STD

what does Patty think of donkeys at last?Patty最后对驴是怎么想的

matlab 中 stem函数用于显示茎秆图，语法有：1、stem（Y） 将数据序列Y从x轴到数据值按照茎秆形式画出，以圆圈终止。如果Y是一个矩阵，则将其每一列按照分隔方式画出。2、stem（X，Y）在X的指定点处画出数据序列Y。3、stem（...，＇filled＇） 以实心的方式画出茎秆图。4、stem（...，＇LINESPEC＇） 按指定的线型画出茎秆图及其标记。扩展资料：stem() 函数主要用来画一张整图，画出来函数在数学层面属于离散函数。在均匀平铺的位置上生成轴对象，在Figure对象层次上起作用。在matlab的命令窗口中输入doc stem或者help stem就可以获得stem() 函数的帮助信息。缺省情况下，stem() 函数将图像中的数据存储为双精度类型(double)，64位浮点数，所需存储量很大；stem() 函数还支持另一种类型无符号整型(uint8)，即图像矩阵中每个数据占用1个字节。

对于离散序列，MATLAB用stem( )命令实现其绘制。 示例4：绘制从0 ~ 4 范围的正弦函数序列。 t = 0 : pi/10 : 4*pi;y = sin(t);stem(t, y);

单位阶跃函数。

matlab中step的作用是模块功能。STEP 可划分为两部分，STEP标准的数据模型和工具。数据模型包括通用集成资源、应用集成资源、应用协议；工具包括描述方法、实现方法、一致性测试方法和抽象测试套件。其中资源信息模型定义了开发应用协议基础的数据信息，包括通用的模型和支持特定应用的模型。它提供了一种不依赖具体系统的中性机制，只在实现产品数据的交换和共享。这种描述的性质使得它不仅适合于交换文件，也适合于作为执行和分享产品数据库和存档的基础。扩展资料step使用规则1、数据类型。counter 的数据类型通常是 Integer，但也可以是任何支持大于或等于 (>=)、小于或等于 (<=)、加法 (+) 和减法 (-) 运算符的类型。只要支持所有这些运算符，该数据类型甚至可以是用户定义的类型。2、start、end 和 step 表达式通常计算为 Integer 类型，但也可以计算为任何拓宽到 counter 的类型的数据类型。如果要将用户定义的类型用于 counter，这意味着您可能必须定义 CType 转换运算符，来将 start、end 或 step 的类型转换为 counter 的类型。参考资料来源：百度百科—STEP

如果要求出表达式，那么你应该先用符号运算。>> syms x>> strain1=x.^2;>> cycle1= strain1+diff(strain1)cycle1 =x^2 + 2*x然后再用x=x0(比如x=[1,2])替代原来的x来得到cyclex0=[1,2];>> strain=subs(strain1,x,x0)strain = 1 4>> cycle=subs(cycle1,x,x0)cycle = 3 8

struct在matlab中是用来建立结构体数组的。通常有两种用法：s = struct("field1",{},"field2",{},...) 这是建立一个空的结构体，field1，field2是域名称s = struct("field1",values1,"field2",values2,...) 这是在建立结构体的时候就赋值，values1,values2就是针对field1,field2所赋的值。例如：s = struct("type",{"big","little"},"color",{"red"},"x",{3 4})这个结构体有三个域"type","color","x"，根据英文意义可以看出分别是类型，颜色和值的意思。类型可以取的值有"big","little"，也就是类型有两种，分别是大和小。颜色只有"red"红色，而值x可以是3或者4matlab帮助里解释的很详细，你自己可以去check一下

b=find(a),a是一个矩阵,查询非零元素的位置，如果X是一个行向量，则返回一个行向量，否则，返回一个列向量。如果X全是零元素或者是空数组，则返回一个空数组，例子如下所示,也可以用b=find(a>2)，这句的意思是在a中找到比较2大的元素。 扩展资料 　　find在英语中的用法归纳如下： 　　一、作动词： 　　1.find的基本意思是“发现，找到”，作“发现”解时指发现已存在而以前不知道的事物，作“找到”解时指找到了所寻找的东西，知道了其下落，强调寻找的结果。引申可作“查明”“偶然遇见”“供给”“到达”“努力获得”“知道…是有的”等解。用于法律术语还可表示“裁决”。 　　2.find主要用作及物动词，可接名词、代词、带疑问词的动词不定式或从句作宾语，也可接双宾语，其间接宾语可以转换为介词for的宾语，可用于被动结构。 　　3.find还可接以动词不定式或“(to be+) adj./v -ed/ v -ing/ prep. -phrase”充当补足语的复合宾语。接现在分词时表示动作正在进行，接过去分词表示动作已经完成，接动词不定式则表示看到动作的全部过程，不定式符号to可有可无，to则不可省略，但可省略to be。 　　4.find如接以动词不定式、动名词作宾语或以形容词、介词短语充当宾语补足语的复合结构时，一般须采用形式宾语结构。 　　5.find后接that从句，可表示发现一个事实、道理; 后接out，可表示故意去找而找到某件难于找到的东西或经过打听、询问之后搞清楚某事，强调付出了一定的"努力。 　　6.find是瞬间动词，不可用于进行体，也不能与表示一段时间的状语连用。 　　二、作名词： 　　1.find用作名词意思是“发现”，指抽象的行为，是不可数名词，也可作“发现物”解，尤指有价值的东西，是可数名词。 　　2.find还可表示“被发现有惊人能力的人”，是可数名词。

从你的需要查找的数据来看，语句u=find（a>=max(a)）所得到的应该只是得到a向量中最大值所在原始向量中的序号。对于行向量或列向量，u=find（a>=max(a)）得到的都只是一个值，a(u)都代表的是这个值，我刚试了下没出现楼主出现的问题。

题主最好给出一个稍微具体点的应用场景，可能有更加优化的方法。我自己构造一个简单的例子如下：matlab:a=find(b>0)python:a=[iforiinrange(len(b))ifb[i]>0]另外注意，matlab里的数组索引从1开始，和python不同。

find函数用于返回所需要元素的所在位置 (位置的判定：在矩阵中，第一列开始，自上而下，依次为1，2，3...,然后再从第二列，第三列依次往后数)find（A）返回矩阵A中非零元素所在位置>> A = [1 0 4 -3 0 0 0 8 6];>> X = find(A)X =1 3 4 8 9find（A>5）返回矩阵A中大于5的元素所在位置>> find(A>5)ans =8 9[i,j,v]=find(A) 返回矩阵A中非零元素所在的行i,列j,和元素的值v(按所在位置先后顺序输出)>> A=[3 2 0; -5 0 7; 0 0 1];>> [i,j,v]=find(A)i =12123j =11233v =3-5271find（A>m,4）返回矩阵A中前四个数值大于m的元素所在位置也是自己理解，如有不正之处，还望指正！希望我可以帮上点小忙！

Grasshopper是数据化设计方向的主流软件之一软件，MATLAB主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。Grasshopper（简称GH）是一款可视化编程语言，它基于Rhino平台运行，是数据化设计方向的主流软件之一，同时与交互设计也有重叠的区域。与传统设计方法相比，GH的最大的特点有两个：一是可以通过输入指令，使计算机根据拟定的算法自动生成结果，算法结果不限于模型，视频流媒体以及可视化方案。二是通过编写算法程序，机械性的重复操作及大量具有逻辑的演化过程可被计算机的循环运算取代，方案调整也可通过参数的修改直接得到修改结果，这些方式可以有效的提升设计人员的工作效率。MATLAB是matrix和laboratory两个词的组合，意为矩阵工厂（矩阵实验室），软件主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式。

find 函数如何使用在士气球求助啊啊五环岛东路辅路

试试用lingo吧，或者下MATLAB规划的工具用箱试试。MATLAB不好做0-1规划。

把人家机器的gamma.m函数拷到你的文件目录下

程序一：syms tdsolve("Dy+20*y=0","y(0)=1")结果：exp(-20*t)程序二：t=linspace(0,1,1000);plot(t,exp(-20*t),"r.-")结果：

function [c, s] = wavefast(x, n, varargin) %WAVEFAST Perform multi-level 2-dimensional fast wavelet transform. % [C, L] = WAVEFAST(X, N, LP, HP) performs a 2D N-level FWT of % image (or matrix) X with respect to decomposition filters LP and % HP. % % [C, L] = WAVEFAST(X, N, WNAME) performs the same operation but % fetches filters LP and HP for wavelet WNAME using WAVEFILTER. % % Scale parameter N must be less than or equal to log2 of the % maximum image dimension. Filters LP and HP must be even. To % reduce border distortion, X is symmetrically extended. That is, % if X = [c1 c2 c3 ... cn] (in 1D), then its symmetric extension % would be [... c3 c2 c1 c1 c2 c3 ... cn cn cn-1 cn-2 ...]. % % OUTPUTS: % Matrix C is a coefficient decomposition vector: % % C = [ a(n) h(n) v(n) d(n) h(n-1) ... v(1) d(1) ] % % where a, h, v, and d are columnwise vectors containing % approximation, horizontal, vertical, and diagonal coefficient % matrices, respectively. C has 3n + 1 sections where n is the % number of wavelet decompositions. % % Matrix S is an (n+2) x 2 bookkeeping matrix: % % S = [ sa(n, :); sd(n, :); sd(n-1, :); ... ; sd(1, :); sx ] % % where sa and sd are approximation and detail size entries. % % See also WAVEBACK and WAVEFILTER. % Copyright 2002-2004 R. C. Gonzalez, R. E. Woods, & S. L. Eddins % Digital Image Processing Using MATLAB, Prentice-Hall, 2004 % $Revision: 1.5 $ $Date: 2003/10/13 01:14:17 $ % Check the input arguments for reasonableness. error(nargchk(3, 4, nargin)); if nargin == 3 if ischar(varargin{1}) [lp, hp] = wavefilter(varargin{1}, "d"); else error("Missing wavelet name."); end else lp = varargin{1}; hp = varargin{2}; end fl = length(lp); sx = size(x); if (ndims(x) ~= 2) | (min(sx) < 2) | ~isreal(x) | ~isnumeric(x) error("X must be a real, numeric matrix."); end if (ndims(lp) ~= 2) | ~isreal(lp) | ~isnumeric(lp) ... | (ndims(hp) ~= 2) | ~isreal(hp) | ~isnumeric(hp) ... | (fl ~= length(hp)) | rem(fl, 2) ~= 0 error(["LP and HP must be even and equal length real, " ... "numeric filter vectors."]); end if ~isreal(n) | ~isnumeric(n) | (n < 1) | (n > log2(max(sx))) error(["N must be a real scalar between 1 and " ... "log2(max(size((X)))."]); end % Init the starting output data structures and initial approximation. c = []; s = sx; app = double(x); % For each decomposition ... for i = 1:n % Extend the approximation symmetrically. [app, keep] = symextend(app, fl); % Convolve rows with HP and downsample. Then convolve columns % with HP and LP to get the diagonal and vertical coefficients. rows = symconv(app, hp, "row", fl, keep); coefs = symconv(rows, hp, "col", fl, keep); c = [coefs(:)" c]; s = [size(coefs); s]; coefs = symconv(rows, lp, "col", fl, keep); c = [coefs(:)" c]; % Convolve rows with LP and downsample. Then convolve columns % with HP and LP to get the horizontal and next approximation % coeffcients. rows = symconv(app, lp, "row", fl, keep); coefs = symconv(rows, hp, "col", fl, keep); c = [coefs(:)" c]; app = symconv(rows, lp, "col", fl, keep); end % Append final approximation structures. c = [app(:)" c]; s = [size(app); s]; %-------------------------------------------------------------------% function [y, keep] = symextend(x, fl) % Compute the number of coefficients to keep after convolution % and downsampling. Then extend x in both dimensions. keep = floor((fl + size(x) - 1) / 2); y = padarray(x, [(fl - 1) (fl - 1)], "symmetric", "both"); %-------------------------------------------------------------------% function y = symconv(x, h, type, fl, keep) % Convolve the rows or columns of x with h, downsample, % and extract the center section since symmetrically extended. if strcmp(type, "row") y = conv2(x, h); y = y(:, 1:2:end); y = y(:, fl / 2 + 1:fl / 2 + keep(2)); else y = conv2(x, h"); y = y(1:2:end, :); y = y(fl / 2 + 1:fl / 2 + keep(1), :); end

将路径和文件名称作为处理文件函数的参数输入即可function MMM(dir)比如写一个文件名称为MMM，里面有个函数名称为MMM，带一个dir的文件路径参数。在matlib窗口中输入MMM（"C:XXX.txt"），应当就可以吧

用resample函数

不能直接修改fs。需要重采样进行转换，resample函数。

采样点增益是指在模拟信号采样过程中，由于采样时序与信号本身的时序不完全匹配，导致采样结果中出现了噪声。采样点增益可以通过在数字信号处理中应用FIR滤波器来减小。通过滤波器设计，可以得到一个滤波器系数向量，利用这个向量可以对信号进行滤波处理，从而消除掉采样点增益所带来的影响。在Matlab中，可以使用fir1函数来设计FIR滤波器。该函数可以根据用户指定的滤波器类型、滤波器阶数和截止频率等参数，自动生成滤波器系数。然后，可以使用filter函数将信号与滤波器系数进行卷积，得到滤波后的信号。需要注意的是，采样点增益是由于采样时序与信号本身的时序不匹配而产生的噪声，因此在采样过程中应尽可能保证采样时序与信号时序的匹配度，以减小采样点增益对信号的影响。此外，在信号处理过程中也可以采用其他方法来消除噪声，例如小波变换、自适应滤波等。

调用函数resample（）就可以了resample为信号降采样处理，理解如下：B=resample(x,90,250); %采样从250Hz降到90Hz，如果250在前,就是插值从90到250,可以看B的长度,250Hz采样4000个数据等于90hz采样1440个数据,这就是降采样。

resample是抽取decimate和插值interp的两个结合具体完成如下操作，先插值9变成250*9Hz然后抽取25变成速率90Hz

resample函数，help一下，有详细介绍。

resample为信号降采样处理，理解如下：b=resample(x,90,250);%采样从250hz降到90hz，如果250在前,就是插值从90到250,可以看b的长度,250hz采样4000个数据等于90hz采样1440个数据,这就是降采样。resample是抽取decimate和插值interp的两个结合具体完成如下操作，先插值90变成250*9hz然后抽取250变成速率90hz。

resample为信号降采样处理，理解如下：B=resample(x,90,250); %采样从250Hz降到90Hz，如果250在前,就是插值从90到250,可以看B的长度,250Hz采样4000个数据等于90hz采样1440个数据,这就是降采样。resample是抽取decimate和插值interp的两个结合具体完成如下操作，先插值90变成 250*9Hz然后抽取250变成速率 90Hz。

这两个不是滤波器,而是重采样(改变采样频率), 虽然它们都使用了低通滤波器来防止重采样过程中的混叠现象. 至于这些低通滤波器(FIR)的具体参数,可以看HELP中的介绍. 比如 interp的:The length of the FIR lowpass interpolating filter is 2*l*r+1.The number of original sample values used for interpolation is 2*l.Ordinarily, l should be less than or equal to 10. The original signal is assumed to be band limitedwith normalized cutoff frequency 0≤alpha≤1,where 1 is half the original sampling frequency (the Nyquist frequency).The default value for l is 4 andthe default value for alpha is 0.5.

区间[0,1]上的均匀分布随机数的系统函数R=rand(n)R=rand(m,n)，产生m*n阶均匀随机矩阵

在matlab中实现函数的抽样用法：一、dyaddown功能：对时间序列进行二元采样，每隔一个元素提取一个元素，得到一个降采样时间序列。 格式：y = dyaddown(x, EVENODD)当EVENODD=0时，从x中第二个元素开始采样（偶采样）；当EVENODD=1时，从x中第一个元素开始采样（奇采样）。 2.y = dyaddown(x)EVENODD缺省，按EVENODD=0二、dyadup功能：对时间序列进行二元插值，每隔一个元素插入一个0元素，得到一个时间序列。 格式：y = dyadup(x, EVENODD)当EVENODD=0时，从x中第二个元素开始采样（偶采样）；当EVENODD=1时，从x中第一个元素开始采样（奇采样）。 2.y = dyadup(x)EVENODD缺省，按EVENODD=0三、interp功能：对时间序列进行整数倍插值，使得时间序列曲线更光滑。 格式： 1.y = interp(x, r)在x中插入一些数据，使得插值后的序列y的长度为x的r倍。 2.y = interp(x, r, l, alpha)插值后得到的序列y的长度为x的r倍。 3.[y, b] = interp(x, r, l, alpha)插值后同时得到一个低通插值滤波器的系数，长度为2rl+1. 说明： x--时间序列 r--插入点的倍数 l--插值滤波器长度alpha--滤波器的截止频率，0<alpha<=1，假设原序列的采样频率之半为1，缺省时l=4，alpha=0.5.y--插值后得到的时间学列b--低通插值滤波器的系数，长度为2rl+1四、downsample功能：对时间序列重采样，在原时间序列中等间隔地取出一些项，得到新序列。 格式：1.y = downsample(x, n)从第一项开始，等间隔n对x采样，得到的序列为y。 2.y = downsample(x, n, phase)从第phase+1项开始，等间隔n对x采样，得到的序列为y，而0<=phase<n.五、decimate功能：对时间序列进行整数倍采样处理，使得时间序列的长度降低。 格式：1.y = decimate(x, r)将时间序列x的采样频率降低为原来的1/r，即length(y)=length(x)/r。在抽取之前，默认地采用了8阶chebyshevI型低通滤波器压缩频带。 2.y = decimate(x, r, n)采用n阶chebyshevI型低通滤波器。 3.y = decimate(x, r, ‘fir")采用30阶的FIR型低通滤波器来压缩频带，对时间序列进行整数倍抽取。 4.y = decimate(x, r， n, ‘fir")指定当对时间序列进行整数倍抽取的时候，采用n点FIR型低通滤波器来压缩频带，对时间序列进行整数倍抽取。 说明： x--时间序列 r--采样要降低的倍数n--指定所采用的chebyshevI型低通滤波器的阶数 ‘fir"--FIR滤波器六、resample功能：对时间序列进行重采样。 格式：1.y = resample(x, p, q)采用多相滤波器对时间序列进行重采样，得到的序列y的长度为原来的序列x的长度的p/q倍，p和q都为正整数。此时，默认地采用使用FIR方法设计的抗混叠的低通滤波器。 2.y = resample(x, p, q, n)采用chebyshevIIR型低通滤波器对时间序列进行重采样，滤波器的长度与n成比例，n缺省值为10.3.y = resample(x, p, q, n, beta)beta为设置低通滤波器时使用Kaiser窗的参数，缺省值为5. 4.y = resample(x, p, q, b)b为重采样过程中滤波器的系数向量。 5.[y, b] = resample(x, p, q)输出参数b为所使用的滤波器的系数向量。 说明： x--时间序列p、q--正整数，指定重采样的长度的倍数。n--指定所采用的chebyshevIIR型低通滤波器的阶数，滤波器的长度与n成比列。 beta--设计低通滤波器时使用Kaiser窗的参数，缺省值为5。

你是不是在看KCF的代码啊。。也做目标跟踪嘛，我也在看。。哈哈哈suffix是后缀的意思，因为数据集中的Jogging视频序列中有两个跟踪目标，所以拆分为Jogging.1和Jogging.2来看，video（end-1:end）指的是video的后两个字符，也就是.1和.2，是要记一下video的后缀

用resample()，重新采样

imresize，该函数用于对图像做缩放处理。在matlab的命令窗口中输入doc imresize或者helpimresize即可获得该函数的帮助信息。函数功能该函数用于对图像做缩放

三次样条插值直接用spline函数做。边界条件加在y的首尾，第一个表示y"(x0)，最后一个表示y"(xt)。如果不加边界条件，默认是not-a-knot边界条件（注意不是自然边界条件）自然边界条件的插值要用csape函数才能得到。如果用interp1，则只能使用spline函数的默认边界条件，即not-a-knot条件。下面是例子x=0:3:9;y=x.*cos(x);xx=linspace(0,9);plot(x,y,"o");%样本点hold on;plot(xx,interp1(x,y,xx,"spline"),"r");%interp1只能使用默认边界条件plot(xx,spline(x,[0 y 0],xx),"r:");%spline可以使用第一类边界条件，这里y"(0)=y"(9)=0pp=csape(x,y,"second");plot(xx,fnval(pp,xx))%第二类边界条件要用csape做，这里自然边界条件legend("样本点","默认边界条件","一阶导为0","自然边界条件","location","south")

使用resample函数。举个例子：B=resample(x，90，250)%，采样从250降到90，如果250在前，就是插值从90到250，可以看B的长度，250采样4000个数据等于90采样1440个数据。

resample为信号降采样处理，理解如下：B=resample(x,90,250);%采样从250Hz降到90Hz，如果250在前,就是插值从90到250,可以看B的长度,250Hz采样4000个数据等于90hz采样1440个数据,这就是降采样。resample是抽取decimate和插值interp的两个结合具体完成如下操作，先插值90变成250*9Hz然后抽取250变成速率90Hz。

我曾经写过的一个程序,包含了画点,拟合,误差评估等你相应的删减一些,就能得到你需要的东西了.希望能看的懂：关键就一个函数：a=polyfit(x,y,n);请去matlab:helppolyfit保证10分钟弄回自己的程序.求人不如求己%tofindtheleast_squaresfitforaninputdataset.%tomakeastraighline,andprintitout.%todetemingtheconstantm,bofy=mx+b.%%clc;%clearthecommandwindowclear;%clearallthevariablehasbeendefineddisp("----thisprogramfindtheleast_sqaresfitofaninputdataset");disp("----astheformofy=a1*xn+a2*x(n-1)+.a(n)*x1+a(n+1)");disp("------------starting!----------------");%inputdatanum%inputtimesofpolyx=input("enterthearrayofxlike[1,2,3,4,.]:");y=input("enterthearrayofylike[1,2,3,4,.]:");n=input("entertimesofpolynomiallike2:")a=zeros(n+1);%tocounthowmanydatainxrow_cloum=size(x);m=row_cloum(2);fori=1:1:mplot(x(i),y(i),"bo");holdon;end%calculateaa=polyfit(x,y,n);y_best=0;%leastsqurausofybest%%plotthedatainputdataasbluecircles%plotthelineinred"-"%createthefittedlineoffigturexmin=min(x);xmax=max(x);t=xmin:0.01:xmax;fori=1:(n+1)y_a=(a(i)*(t.^((n+1)-i)));%elementofy_poly_arrayy_best=y_best+y_a;%sumthemend%rtheerrorarrayfork=1:mr_sum=0;y_kbest=0;fori=1:(n+1)y_k(i)=(a(i)*(x(k).^((n+1)-i)));%elementofy_k=ploy(xi)y_kbest=y_kbest+y_k(i);%sumthemploy(xi)endr_sum=(y(k)-y_kbest)^2+r_sum;end%thepolyis:disp("thepolyisy=a1*xn+a2*x(n-1)+.a(n)*x1+a(n+1):")disp("andthearrayAisA=:");disp(a);disp("the(y_best-y)^2isR=:")disp(r_sum);%figtureymin=min(y_best);ymax=max(y_best);axis=([xmin,xmax,ymin,ymax]);%setplotsettings%addatitleandlegendplot(t,y_best,"r-","LineWidth",2);%holdoff;title("fleast-squaresfity=a1*xn+a2*x(n-1)+.a(n)*x1+a(n+1)");xlabel("fitx");ylabel("fity");legend("inputdata","fittedline",4);gridon;print-dtifflsq_fit.tif;

a是uint8格式的，把a转换成double格式就可以了。surf(x,y,double(a))

‘的含义为转置，*为乘法；n"*k，就代表矩阵n转置后与矩阵k相乘参考一：计算离散时间傅里叶变换，并绘制图形。已知有限长序列x(n)={1,2,3,4,5}。n=-1:3;x=1:5;k=0:500;w=(pi/500)*k;X=x*(exp(-j*2*pi/500)).^(n"*k);subplot(211)plot(w/pi,imag(X)) subplot(212)plot(w/pi,abs(X))grid on

没错啊！很好的

mesh和surf都是画三维 图的函数，不适合画二维图---圆形。

你数据太少了。

plot3是最简单的三维图mesh是网格图surf是曲面图

surf(x,y,z)时，颜色越深（蓝）表示z值越小surf(x,y,z,c)时，颜色越深（蓝）表示c值越小

举个例子[X,Y,Z]=peaks;surf(X,Y,Z);hold onsurf(X,Y,0*Z-20,Z)

%syms alphb betab X Y Z; %这一行没有必要rE=4226;e=870;alphb=linspace(-pi/6,pi/6,11);betab=linspace(-pi/8,pi/4,10);[Ab,Bb]=meshgrid(alphb,betab);x=-(rE*cos(Bb)+e).*sin(Ab);y=(rE*cos(Bb)+e).*cos(Ab);z=rE*sin(Bb);surf(x,y,z) %区分大小写的

你的网格太密了。就是说，你的x,y取太多了。适当减少一些，就可以看得比较清楚。或者放大看效果。

surf创建带有填充的网格图，网格线是不变色的，靠填充的颜色表示值的大小变化。meSh创建网格图，不填充颜色，可以通过网格线的颜色表示数值大小的变化。大概的区别就在这里，还有一些它们的高级应用，它们的差别就很明显了，你可以查看帮助文档。通常如果数据点够密集，网格线就不很明显了，这两个函数作出的图就很相似了。可以

x=1:1:12;y=1:1:5;z=[0.20.240.250.260.250.250.250.260.260.290.250.29;0.270.310.30.30.260.280.290.260.260.260.260.29;0.410.410.370.370.380.350.340.350.350.340.350.35;0.410.420.420.410.40.390.390.380.360.360.360.36;0.30.360.40.430.450.450.510.420.40.370.370.37];surf(x,y,z)

%syms alphb betab X Y Z; %这一行没有必要 rE=4226; e=870; alphb=linspace(-pi/6,pi/6,11); betab=linspace(-pi/8,pi/4,10); [Ab,Bb]=meshgrid(alphb,betab); x=-(rE*cos(Bb)+e).*sin(Ab); y=(rE*cos(Bb)+e).*cos(Ab); z=rE*sin(Bb); surf(x,y,z) %区分大小写的

主要是载入网格数据，比如：clc;clearall;closeall;[X,Y,Z]=peaks(30);figure;surf(X,Y,Z)

1、首先，应该说明一下，EZPLOT在多数情况下画出的图还是比较准确的。不能因为某些情况下画出的图不准确就对它失去信心，毕竟EZ*这类函数用起来确实很方便。如果对于绘图结果有疑问，可考虑换用其它等效的表达方式试一试，看结果是否会有改进。2、对于你说的这个函数，可以使用下述写法：ezplot("x^2^(1/3)+y^2^(1/3)-1")画出来的图应该是你想要的了吧？3、说到原因，我也有点迷茫了。好像这涉及到一个非常基本的问题：-1开三次方结果应该是什么？按照常理，(-1)^(1/3)是不是应该等于-1？但是你在MATLAB里面试一试

可这样：t = 0 :0.01: tan(25/180*pi);[X,Y,Z] = cylinder(t);cdata=cat(3,ones(size(X)),zeros(size(X)),zeros(size(X)));%红色%cdata=cat(3,zeros(size(X)),ones(size(X)),zeros(size(X)));%绿色%cdata=cat(3,zeros(size(X)),zeros(size(X)),ones(size(X)));%兰色surf(X,Y,2*Z,cdata);%surf(X,Y,2*Z);%shading interpset(gca,"DataAspectRatio", [1 1 1]);

u=[-pi/2:pi/60:pi/2];u(:,[1,61])=[];v=[0:pi/60:2*pi];[U,V]=meshgrid(u,v);X=3*tan(U).*cos(V);Y=4*tan(U).*sin(V);Z=5*sec(U);surf(X,Y,Z);hold on;surf(X,Y,-Z);axis equal;

matlab中mesh（）和surf（）函数的主要区别是：mesh（）用于绘制不是特别精细的三维曲面网格图。同一层面的线条用相同的颜色表示。surf（）用于绘制比较光滑的三维曲面网格图。各线条之间的补面用颜色填充。

是matlab里面的专用函数名，翻译不出来啦，就是画三维曲面图surf(x,y,z)

1、使用view切换视角；2、不知道你所说的“插入云图”是什么意思？如果想手动修改颜色标识范围可能有些困难（整体上换用另一组调色板可以使用colormap）；3、不显示网格线包括两个方面，一是grid off，另一是把surface对象的边设为无色（需要获得对象的句柄）；4、使用File菜单的Export或Save as（不同版本的菜单不一样）保存。当然，如果需要，也可以写代码自动保存。示例代码（用peaks画曲面作为演示）：peaksview(0,90)colorbars=findall(gcf,"type","surface");set(s,"EdgeColor","none")grid offshading interp您可以逐句运行，看看每一句起到什么样的作用。希望对你有帮助，有问题可继续追问。

三维曲面（色）图 surf 和 surfc 是通过矩形区域来观测数学函数的函数。surf和surfc能够产生由X、Y、Z指定的有色参数化曲面，即三维有色图。 当x = 1:n、y = 1:m，并且[m,n]=size(Z)时，surf(Z)会产生一个矩阵Z的z部分的三维遮罩层，这里Z是一个定义在几何矩形网格内的单值函数。 surf(X,Y,Z)同样产生矩阵Z的有色遮罩层，XY可以是有xy定义的向量或矩阵 surf(X,Y,Z,C)是产生一个由C定义颜色的矩阵Z的有色遮罩层 例子k = 5;n = 2^k-1;[x,y,z] = sphere(n);c = hadamard(2^k);figure(1)surf(z);figure(2)surf(x,y,z);axis equalfigure(3)surf(x,y,z,c);colormap([1 1 0; 0 1 1])axis equal

xi=-10:0.5:10;yi=-10:0.5:10;[x,y]=meshgrid(xi,yi);z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2);surf(x,y,z)

画三维的，彩色的 面

x=1:1:12; y=1:1:5;z=[0.2 0.24 0.25 0.26 0.25 0.25 0.25 0.26 0.26 0.29 0.25 0.29; 0.27 0.31 0.3 0.3 0.26 0.28 0.29 0.26 0.26 0.26 0.26 0.29; 0.41 0.41 0.37 0.37 0.38 0.35 0.34 0.35 0.35 0.34 0.35 0.35; 0.41 0.42 0.42 0.41 0.4 0.39 0.39 0.38 0.36 0.36 0.36 0.36; 0.3 0.36 0.4 0.43 0.45 0.45 0.51 0.42 0.4 0.37 0.37 0.37];surf(x,y,z)

【MATLAB】用surf函数画出来：三维曲面（色）图。clear all[x,y]=meshgrid(linspace(-15,15));%设定xy范围z=sin((x.^2+y.^2).^0.5)./((x.^2+y.^2).^0.5);figure(1)surf(x,y,z)xlabel("X"),ylabel("Y"),zlabel("Z");figure(2)surf(x,y,z,"LineStyle","none")alpha(0.7)%设定透明度xlabel("X"),ylabel("Y"),zlabel("Z");扩展资料：曲面图具有下列图表子类型：三维曲面图 三维曲面图在连续曲面上跨两维显示数值的趋势。曲面图中的颜色并不代表数据系列，而是代表数值间的差别。三维曲面图（框架图） 不带颜色的三维曲面图称为三维曲面图（框架图）。注释：不带颜色的三维曲面图（框架图）不易于理解。您可能需要改用三维曲面图。曲面图和曲面图（俯视框架图） 曲面图和曲面图（俯视框架图）是曲面图的上视图。在曲面图中，颜色代表特定的数值范围。曲面图（俯视框架图）不显示颜色。注释：曲面图和曲面图（俯视框架图）不易于理解。参考资料来源：百度百科-曲面图

mesh 三维网格图；surf 三维表面图。

函数surf的作用：画三维曲面（色）图，起作用与mesh相似。surf(X,Y,Z)　X、Y、Z中Z通常是X,Y的函数，即Z(X,Y)。X、Y通常是通过调用meshgrid函数生成的数据网格。surf(Z)　surf(...,C)　surf(...,"PropertyName",PropertyValue,...)　surf(axes_handles,...)　相关函数：mesh，meshc,meshzxgrid=-10:0.1:10;ygrid=-10:0.1:10;[x,y]=meshgrid(xgrid,ygrid);z=0.5.*(x-exp(0.0154.*y-2));surf(x,y,z);

surf 三维着色表面图 三维表面图 表面图mesh 三维网线图 网线图

三维曲面（色）图 surf 和 surfc 是通过矩形区域来观测数学函数的函数。surf和surfc能够产生由X、Y、Z指定的有色参数化曲面，即三维有色图。 当x = 1:n、y = 1:m，并且[m,n]=size(Z)时，surf(Z)会产生一个矩阵Z的z部分的三维遮罩层，这里Z是一个定义在几何矩形网格内的单值函数。 surf(X,Y,Z)同样产生矩阵Z的有色遮罩层，XY可以是有xy定义的向量或矩阵 surf(X,Y,Z,C)是产生一个由C定义颜色的矩阵Z的有色遮罩层 例子k = 5;n = 2^k-1;[x,y,z] = sphere(n);c = hadamard(2^k);figure(1)surf(z);figure(2)surf(x,y,z);axis equalfigure(3)surf(x,y,z,c);colormap([1 1 0; 0 1 1])axis equal

1、解压tensor_toolbox_2.5.zip2、将解压出来的文件夹tensor_toolbox_2.5拷贝到malab安装文件中的toolbox文件夹下，并将tensor_toolbox_2.5改名为tensor_toolbox3、启动malab，在matlab窗口中中找到Current Folder，从根目录中找到tensor_toolbox文件夹4、在matlab窗口中运行以下命令：addpath(pwd) %<-- Add the tensor toolbox to the MATLAB pathcd met; addpath(pwd) %<-- Also add the met directorysavepath %<-- Save for future MATLAB sessions参考网址及下载地址：http://www.sandia.gov/~tgkolda/TensorToolbox/thankyou-2.5.html

1、解压tensor_toolbox_2.5.zip 2、将解压出来的文件夹tensor_toolbox_2.5拷贝到malab安装文件中的toolbox文件夹下，并将tensor_toolbox_2.5改名为tensor_toolbox 3、启动malab，在matlab窗口中中找到Current Folder，从根目录中找到tensor_tool.

原因之一是因为没有把tensor toolbox 的路径加到matlab工具包下将tensor包拷贝到matlab的安装文件下的toolbox文件目录下，然后在matlab中运行addpath D:/MATLAB/toolbox/tensor toolbox（中间的目录依个人而定） ，然后运行 help ‘tensor"，若成功则安装成功

原因之一是因为没有把tensortoolbox的路径加到matlab工具包下将tensor包拷贝到matlab的安装文件下的toolbox文件目录下，然后在matlab中运行addpathD:/MATLAB/toolbox/tensortoolbox（中间的目录依个人而定），然后运行help‘tensor"，若成功则安装成功

clock是一个获取系统当前时间的函数，是一个行向量。可以直接使用，也可以赋给某个变量名，比如date=clock;clock本身不带参数。元素依次是 年 月 日 时 分秒。至于clock元素的选用，可以像对数组操作一样。比如date(:,4)就是获得当前时间的小时数。

法国龙船干白小龙船干红葡萄酒。酒，不管是“纯粮酿造”还是“食用酒精勾兑”，其中的酒精都是通过粮食发酵而来的。dragonboatlegend是法国龙船干白小龙船干红葡萄酒，产自法国波尔多产区，甄选赤霞珠为原料酿制。

matlab里的support(212)是什么意思？support（）函数的作用就是分区绘图。support(212)是在第2行第1列的位置放置第2张图形。（或第2区第1列的位置放置第2张图形。）

matlab中text 函数1、作用：可用该函数在图形中指定的位置上显示字符串。2、使用方法： 如 text(x,y,"string")：在二维图形中指定的位置(x,y)上显示字符串stringtext(x,y,z,"您好！") ：在三维图形空间中的指定位置(x,y,z)上显示字符串"您好！"text(x,y,z,"string"."PropertyName",PropertyValue…) ：对引号中的文字string定位于用坐标轴指定的位置，且对指定的属性进行设置。

程序没有问题

在布局问题求解中,为了好表达约束条件,需要引入一个决策变量Vik(i表示设备序号,i=1,2,3,....15!K表示第几行,k=1,2,3)因为一个设备只能在一行,而一行中最多布置设备数量不能超过设备总数15当设备i在第k行的时候Vik=1,else Vik=0 注意(i与k是两个不同的表示量)if ( ) Vik=1else Vik=0end%约束1^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^fVik=0;for k=1:3 %表示从第1到第3行循环 fVik=fVik+Vik;endfV1ik=fVik-1; % fV1ik=0就满足约束1%约束2^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^fVik=0;for i=1:15 fVik=fVik+Vik;endfV2ik=15-fVik; %fV2ik>=0就满足约束2现在问题是if 后面括号的程序应该如何写?

因为那个默认的ans是int类型的，四舍五入。所以得到1

图片不存在啊 检查图片是否在当前文件夹下否则 写完整路径 例如 imread("E:lena.bmp" )