ATL

以有理数形式输出结果。。。

MATLAB里的formula的意思是，返回系统内嵌的CFIT、SFIT或FITTYPE对象的公式或自定义函数公式。其使用格式：formula(fun)例如：f = fittype("weibull");formula(f)ans =a*b*x^(b-1)*exp(-a*x^b)

如果lower(val_method)这个变量对应的值为"manual"，并且segment不是元胞数组，那么就输出那段提示，并且使Model这个变量为空数组，然后函数返回并退出。另外，Matlab的逻辑与运算应该是&&，而不是&。

hrollfcoef这个函数不是matlab自带的function [xh] = hrollfcoef(irfn,ipoint,sr,alfs,ncc)%****************** variables *************************% irfn : Number of symbols to use filtering% ipoint : Number of samples in one symbol% sr : symbol rate% alfs : rolloff coeficiense% ncc : 1 -- transmitting filter 0 -- receiving filter% *****************************************************xi=zeros(1,irfn*ipoint+1);xq=zeros(1,irfn*ipoint+1);point = ipoint;tr = sr ; tstp = 1.0 ./ tr ./ ipoint;n = ipoint .* irfn;mid = ( n ./ 2 ) + 1;sub1 = 4.0 .* alfs .* tr; % 4*alpha*R_sfor i = 1 : n icon = i - mid; ym = icon; if icon == 0.0 xt = (1.0-alfs+4.0.*alfs./pi).* tr; % h(0) else sub2 =16.0.*alfs.*alfs.*ym.*ym./ipoint./ipoint; if sub2 ~= 1.0 x1=sin(pi*(1.0-alfs)/ipoint*ym)./pi./(1.0-sub2)./ym./tstp; x2=cos(pi*(1.0+alfs)/ipoint*ym)./pi.*sub1./(1.0-sub2); xt = x1 + x2; % h(t) plot((1:length(xh)),xh) else % (4alphaRst)^2 = 1plot((1:length(xh)),xh) xt = alfs.*tr.*((1.0-2.0/pi).*cos(pi/4.0/alfs)+(1.0+2.0./pi).*sin(pi/4.0/alfs))./sqrt(2.0); end % if sub2 ~= 1.0 end % if icon == 0.0 if ncc == 0 % in the case of receiver xh( i ) = xt ./ ipoint ./ tr; % normalization elseif ncc == 1 % in the case of transmitter xh( i ) = xt ./ tr; % normalization else error("ncc error"); end % if ncc == 0 end % for i = 1 : n %******************** end of file ***************************网上找的，你看看能不能拼到你那个程序里去

在for语句中使用plot画图，画出来是只有点。因为plot(SNR,BER,"-ro")中的BER变量为单一数值。所以如想绘画出直线段的话，则应考虑数值变量，即BER变量为一数组变量。因此，题主的代码应作如下改动：1、for循环语句头改动，即为for i=0:1:152、BER=nerr/(cols-rows)代码改动，即为BER（i）=nerr/(cols-rows)3、plot函数位置改动，放在for循环语句外，即。。。end%LDPC性能图plot(SNR,BER,"-ro")set(gca,"Yscale","log");hold on

听不懂 .

能解释一下为什么要用模拟退火法吗？这个函数的最小值直接观察就能看出来在x=y=0处，或者用fminunc也可以求解：x=fminunc(@(x)sin(x(1)*x(2))+x(1)^2+x(2)^2,[0 0])当然，用模拟退火算法也很简单，换一下函数即可：x=simulannealbnd(@(x)sin(x(1)*x(2))+x(1)^2+x(2)^2,[0 0])

可以通过外接矩形来得到直径信息，参考regionprops函数。

需建立整车模型。 准备好Matlab/Sinilmk 模型，本教程以 Cause 软件中自带的一个模型，将其命名为 Lxm.mdlo。设置端门，注意：DLL 文件中的输入与输出使用的是 Simuluik 中的印T与"Out”模块，此处与 API 文件中的输入输出接1不同，之后就可完成了。

在实际的工作和生活过程中，优化问题无处不在，比如资源如何分配效益最高，拟合问题，最小最大值问题等等。优化问题一般分为局部最优和全局最优，局部最优，就是在函数值空间的一个有限区域内寻找最小值；而全局最优，是在函数值空间整个区域寻找最小值问题。 matlab中的提供的传统优化工具箱（Optimization Tool），能实现局部最优，但要得全局最优，则要用全局最优化算法（Global Optimization Tool），主要包括： GlobalSearch 全局搜索和 MultiStart 多起点方法产生若干起始点，然后它们用局部求解器去找到起始点吸引盆处的最优点。 ga 遗传算法用一组起始点（称为种群），通过迭代从种群中产生更好的点，只要初始种群覆盖几个盆，GA就能检查几个盆。 simulannealbnd 模拟退火完成一个随机搜索，通常，模拟退火算法接受一个点，只要这个点比前面那个好，它也偶而接受一个比较糟的点，目的是转向不同的盆。 patternsearch 模式搜索算法在接受一个点之前要看看其附近的一组点。假如附近的某些点属于不同的盆，模式搜索算法本质上时同时搜索若干个盆。 下面我就一些具体例子，来说明各种优化方法： 可以看出，初值x0不同，得到的结果截然不同，这说明这种求解器，能寻找局部最优，但不一定是全局最优，在起点为8时，取得全局最优。 我们换一种求解器：fminbound,这种求解器不需要给点初值。 因此全局最优的方法能够获取全局最优。 结果：最小二乘拟合结果误差较大 可以看出全局优化结果较好，误差较小。 这种算法的运行时间：Elapsed time is 6.139324 seconds. 使用并行计算的方式解决 结果：14 out of 100 local solver runs converged with a positive local solver exit flag. Elapsed time is 4.358762 seconds.Sending a stop signal to all the labs ... stopped.可以看出，运行时间减少，提高了效率。 这种方法只能寻找局部最优。 现在用全局优化算法：

ANSYS就可以啊，不过ANSYS是一种专家软件，比较熟悉才能最大程度的发挥它的作用，ANSOFT相对而言就比较人性化，也比较傻瓜化点。不过如果想锻炼自己，还是用C吧，fortran也行，有限元里面fortran的例子确实非常多。

把电脑系统时间改到2016年

算了你改时间吧什么版本啊

另存为license.lic覆盖到安装目录下/licenses/目录中，亲测2010b和2012b都好用，建议文件名与之前的一致。%%%%%%%%%%%%%%%%%%%%%%%%%%%%INCREMENT Aerospace_Blockset MLM 99 permanent uncounted A05070F00D1EB1F92326 VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=216 SN=888888 TS_OKINCREMENT Aerospace_Toolbox MLM 99 permanent uncounted 6090F0C08395D4289512 VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=237 SN=888888 TS_OKINCREMENT Bioinformatics_Toolbox MLM 99 permanent uncounted 40E0B0406DE56D23A426 VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=190 SN=888888 TS_OKINCREMENT Cert_Kit_IEC MLM 99 permanent uncounted 9080309041D2DCCB2B10 VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=179 SN=888888 TS_OKINCREMENT Communication_Blocks MLM 99 permanent uncounted 80E010304ACCEEB5E0AA VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=77 SN=888888 TS_OKINCREMENT Communication_Toolbox MLM 99 permanent uncounted C0003000770A1A086530 VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=162 SN=888888 TS_OKINCREMENT Compiler MLM 99 permanent uncounted 60D0C0F0DC03C2F72FF7 VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=232 SN=888888 TS_OKINCREMENT Control_Toolbox MLM 99 permanent uncounted 6020E0B00B7ECE0893BB VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=7 SN=888888 TS_OKINCREMENT Curve_Fitting_Toolbox MLM 99 permanent uncounted 10703080061FF9DA5A81 VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=229 SN=888888 TS_OKINCREMENT Data_Acq_Toolbox MLM 99 permanent uncounted E0B0E0E05124CF4A1A8D VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=21 SN=888888 TS_OKINCREMENT Database_Toolbox MLM 99 permanent uncounted 0010D0B02CA7353F8314 VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=230 SN=888888 TS_OKINCREMENT Datafeed_Toolbox MLM 99 permanent uncounted 40905060D1ECD939538F VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=185 SN=888888 TS_OKINCREMENT Dial_and_Gauge_Blocks MLM 99 permanent uncounted A0F05070AA9AEBF40588 VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=226 SN=888888 TS_OKINCREMENT Distrib_Computing_Toolbox MLM 99 permanent uncounted 5000D00031D2C6B89F9C VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=229 SN=888888 TS_OKINCREMENT EDA_Simulator_Link MLM 99 permanent uncounted E070C0C02B4B90B5A8D0 VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=205 SN=888888 TS_OKINCREMENT Econometrics_Toolbox MLM 99 permanent uncounted 505010F00EB1289A73BA VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=149 SN=888888 TS_OKINCREMENT Embedded_IDE_Link MLM 99 permanent uncounted 90E01090A9BDA4FF67D7 VENDOR_STRING=QQ=47399897 HOSTID=ANY ck=254 SN=888888 TS_OK

你看看这个 http://www.piaodown.com/down/soft/22826.htm

MATLAB软件，基本是电子类和计算机类专业必备的编程软件。其强大的矩阵功能不用说出来，学习这些专业的的同学知道这个软件的重要性，所以，在这里给说明一下如何下载并且安装MATLAB软件。工具/原料MATLAB7.1方法/步骤1首先，可以在作者我的百度云网盘里找到MATLAB7.1文件夹里有下载链接，用下载工具打开并下载即可（至于怎么找到我的百度云网盘，网上搜一搜教程就有，这里就不说了）。2下载完后如图界面显示，首先先解压PLP文档，解压至当前文件夹（放着就行，暂时用不到，但是后面会用到）。3打开MATLAB7.1 CD1（选择winRAR方式打开），如下图界面显示，然后看见最下面的程序setup.exe，双击，选择安装，安装界面会出现一堆英文，不用管，点击next4之后会出现一个让你输入姓名和公司名的页面的（随便输入）将步骤2解压的PLP文件夹打开选择PLP笔记本文档，打开将里面的一串数字编号，复制粘贴到安装界面的PLP框即可。完成后选择NEXT就行。5然后会出现使用协议之类的界面，选择YES，在选择NEXT即可。点击完NEXT后可能会弹出一个网页，不用管直接关掉，继续安装界面，选择Typical自定义安装，然后NEXT。6然后选择你想要安装的位置，然后点击NEXT进入安装页面，接下来是比较麻烦的步骤了。7安装了一段时间后会弹出让你插入CD2的提示，选择你刚开始下载MATLAB的所放置的文件夹，这时你在下载好的文件中将CD2文件解压（就解压在你下载的地方），然后点击OK，千万不要点击SKIP CD。然后就会继续安装。8继续安装后再过一段时间会提示安装CD3，这时的步骤和步骤7操作是一样的，但是先把CD2解压的文件移动到一个新文件加在解压CD3在当前文件夹。CD1、CD2、CD3都搞定后则会继续安装，如下图二。步骤阅读9等待一段时间后，会出现两个界面，第一个直接选择NEXT，第二个是安装完成界面选择FINISH即可。10安装完成后，双击MATLAB发现出现如图界面，不要着急，因为我的是WIN7系统，而软件默认是Windows XP系统，所以右键选择属性，在属性界面选择兼容性，在兼容性选择Windows Vista,然后在打开MATLAB，你会发现可以打开了。步骤阅读END注意事项安装完要是打不开不要着急，仔细看最后一步。根据每个人系统不同，选择不同的兼容版本，WIN7选择Windows Vista就可以了

有关更新 MATLAB 的许可证(license)的步骤，请参阅下面的信息。注意：如果您的许可证中添加了新产品，则您将需要运行安装过程来将这些新产品添加到 MATLAB 中。有关此过程的更多帮助，请参阅相关解决方案 1-1CBD3 -“如何将其他工具箱安装到现有 MATLAB 安装中？”。对于 R2008a 及更高版本：==========如果您要更新计算机上的同一许可证并且具有 Internet 连接，请转到 MATLAB 的“帮助”(Help)菜单并选择“许可”(Licensing)>“更新当前许可证”(Update Current License)来自动下载最新许可证文件。否则，您可以通过重新激活来更新MATLAB。要激活MATLAB，您将需要以下各项之一：- Internet 连接和拥有激活权限的 MathWorks 帐户- Internet 连接和您的许可证的激活密钥- 您的许可证的更新的 license.dat 文件如果您确实满足这些选项之一，请与您的许可证管理员联系以获得许可证文件或激活密钥。要激活 MATLAB，请执行以下操作：1.启动 MATLAB，转到“帮助”(Help) 菜单并选择“许可”(Licensing) >“激活软件”(Activate Software)。如果您无法启动 MATLAB，可以手动启动激活客户端。有关说明，请参阅相关解决方案 1-614OPR -“如何手动启动 MATLAB 激活客户端？”。2.如果您具有 Internet 连接，请选择“使用 Internet 自动激活”(Activate automatically using the Internet)，然后点击“下一步”(Next)。按照剩余页面上的说明输入适当信息。如果您不具有 Internet 连接，请选择“不使用 Internet 手动激活”(Activate manually without using the Internet) 选项，然后点击“下一步”(Next)。在提示时输入许可证文件的路径或浏览到该文件。点击“下一步”(Next) 完成激活。对于 R2007b 及更早版本：==========对于 Windows 安装：对于 MATLAB R14 至 R2007b，您将需要 PLP 来更新 MATLAB。如果您没有 PLP，请与您的许可证管理员联系来获取它。在您得到 PLP 之后：1.从 MATLAB CD 或下载的安装文件（确保除安装程序本身外还下载了您的许可证上的每个产品文件）启动 MATLAB 安装程序。2.选择“输入新的个人许可证密码 (PLP) 来更新许可证而不安装任何内容”(Enter a new Personal License Password (PLP) to update your license without installing anything) 而非“安装”(Install) 选项，然后点击“下一步”(Next)。3.在出现提示时，输入 PLP 并按照安装程序的指示操作，直到完成。对于 UNIX/Linux/Mac 安装：对于 R2008a 之前的版本，您将需要您的 MATLAB 版本的许可证文件来更新 MATLAB。如果您没有此文件，请与您的许可证管理员联系。在您拥有许可证文件之后，按照下面的说明更新您的 MATLAB：1.在 $MATLABROOT/etc（其中 $MATLABROOT 是 MATLAB 安装目录）中找到现有的 license.dat 文件并将该文件重命名为 license.old2.使用文本编辑器打开 license.old。3.删除文件中除了以 SERVER 和 DAEMON 开头的行以外的所有其他行。4.粘贴您从许可证管理员处获得的许可证信息。5.将新文件另存为 license.dat6.在保存文件之后，您将需要在 Unix 和 Linux 计算机上重新启动许可证管理器。资料来源：http://cn.mathworks.com/matlabcentral/answers/115947-matlab

try again~

空着就可以,没必要输入的,不输入也可以安装的如果你需要那我就提供几组吧14-35392-57842-08484-51918-65127-20615-58251-25863-41734-36749-62468-18588-17416-34304-09016-01360-60742-14636-39593-62798-29971-50921-64950-39730-63578-40265-12488-31210-33908-42704-31275-23962-29047 14-13531-19296-24560-20147-44308-03958-00844-31365-30221-49822-53101-49371-40464-62160-05002-40310-37359-48633-64466-15884-57760-64649-27202-63216-21440-58183-16381-31228-38995-55241-00649-11266-22286 14-13079-17698-23560-45609-07352-05352-05350-56924-06021-29849-49518-65460-21329-16030-08711-42918-36765-10752-22284-40230-60254-23812-15890-10926-28557-33310-64529-52031-44719-51491-31733-52511-04580 14-33171-17783-64592-35456-50316-07400-44955-36962-37980-13959-26362-17710-64056-38199-63996-17554-20508-45692-11950-26271-04118-15594-56475-57425-15435-47845-62908-17443-47971-65086-33513-16377-01657 14-24872-36179-32308-22113-61986-33825-00564-45479-60426-10395-51330-19488-62201-37785-19497-45389-18974-51073-03706-04875-59691-49786-28969-00719-61582-14402-53787-33213-56814-33775-57022-14254-56498

Matlab 7 (R14) 注册码1：14-58204-39252-07634-11570-16849-09455-22809-05445-13616-29058-08276-06885-12215-41987-21894-60423-57622-18647-58411-24238-20443-59027-07209-27706-28292-14609-15393-48293-13036-12293-43713-57876-43362 Matlab 7 (R14) 注册码2：14-44889-04614-04275-46147-23559-43066-41714-23083-65272-04997-17469-27919-17226-59862-27901-53983-56217-20094-53460-62647-58166-24499-35558-19511-44882-53016-25658-61109-03776-34505-00776-15813-07183 Matlab 7 (R14) 注册码3：14-02863-32167-49274-14620-55383-23033-26960-31585-34411-64505-59377-01535-25859-02729-42340-44002-31180-19826-51572-37426-25833-53451-02530-20898-18863-41455-29228-13667-31335-59199-04825-64974-59539 Matlab 7 (R14) 注册码4：14-31062-57999-64507-28421-43456-35967-55178-20933-12777-34026-44684-07146-17266-64175-62985-50264-38373-35045-48372-03550-51628-06609-24618-64094-55458-49747-04848-24494-63995-46820-01807-20764-37086 Matlab 7 (R14) 注册码5：14-17107-22787-29968-14354-40195-52542-37833-61505-03296-57542-10587-61927-24639-23185-47049-43791-51460-27360-26190-65454-27576-01212-32058-00132-19787-28479-08728-21126-59415-29872-44199-41000-65357Mathworks Matlab V7.0 R14 PLP 14-35392-57842-08484-51918-65127-20615-58251-25863-41734-36749-62468-18588-174 16-34304-09016-01360-60742-14636-39593-62798-29971-50921-64950-39730-63578-402 65-12488-31210-33908-42704-31275-23962-29047 14-13531-19296-24560-20147-44308-03958-00844-31365-30221-49822-53101-49371-404 64-62160-05002-40310-37359-48633-64466-15884-57760-64649-27202-63216-21440-581 83-16381-31228-38995-55241-00649-11266-22286 14-13079-17698-23560-45609-07352-05352-05350-56924-06021-29849-49518-65460-213 29-16030-08711-42918-36765-10752-22284-40230-60254-23812-15890-10926-28557-333 10-64529-52031-44719-51491-31733-52511-04580 14-33171-17783-64592-35456-50316-07400-44955-36962-37980-13959-26362-17710-640 56-38199-63996-17554-20508-45692-11950-26271-04118-15594-56475-57425-15435-478 45-62908-17443-47971-65086-33513-16377-01657 14-24872-36179-32308-22113-61986-33825-00564-45479-60426-10395-51330-19488-622 01-37785-19497-45389-18974-51073-03706-04875-59691-49786-28969-00719-61582-144 02-53787-33213-56814-33775-57022-14254-56498

guou

PLP是注册序列号的意思！我这里有7.0的序列号 14-35392-57842-08484-51918-65127-20615-58251-25863-41734-36749-62468-18588-17416-34304-09016-01360-60742-14636-39593-62798-29971-50921-64950-39730-63578-40265-12488-31210-33908-42704-31275-23962-29047 14-13531-19296-24560-20147-44308-03958-00844-31365-30221-49822-53101-49371-40464-62160-05002-40310-37359-48633-64466-15884-57760-64649-27202-63216-21440-58183-16381-31228-38995-55241-00649-11266-22286 14-13079-17698-23560-45609-07352-05352-05350-56924-06021-29849-49518-65460-21329-16030-08711-42918-36765-10752-22284-40230-60254-23812-15890-10926-28557-33310-64529-52031-44719-51491-31733-52511-04580 14-33171-17783-64592-35456-50316-07400-44955-36962-37980-13959-26362-17710-64056-38199-63996-17554-20508-45692-11950-26271-04118-15594-56475-57425-15435-47845-62908-17443-47971-65086-33513-16377-01657 14-24872-36179-32308-22113-61986-33825-00564-45479-60426-10395-51330-19488-62201-37785-19497-45389-18974-51073-03706-04875-59691-49786-28969-00719-61582-14402-53787-33213-56814-33775-57022-14254-56498 你挑一个看能不能用

求复数虚部的函数，imag(2+3i)ans = 3

这是由于两矩阵相乘时内维不相造成的，检查d1,d2的维数是否有问题

clc,clearfun=inline("c(1)./(1+((c(1)-37779.1)/3779.1)*exp(-c(2)*c(1)*x))","c","x");x=[1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011];y=[3779.10 3824.69 3860.12 3905.01 3953.22 3990.75 4016.72 4043.06 4083.64 4124.97 4163.07 4199.90 4238.00 4269.50 4304.40 4334.80 4363.70 4389.00 4413.00 4434.80 4456.20 4475.40 4679.91 4728.80 4776.40 4856.80 4925.20 4990.90 5071.80 5154.90 5212.40 5275.50 5446.51 5463.00];c=lsqcurvefit(fun,[0,0],x,y)结果：Maximum number of function evaluations exceeded; increase options.MaxFunEvalsc = -10.0000 -0.0001

quad("1.1687*10^51*1./(sqrt(d)).*exp(-7.24*10^14*d-4.216*10^(-7)./(sqrt(d)))",0.001,10000)ans = 0

貌似ezplot不支持调用变量，R(2),R(3)......这些东西导致了错误。

函数

matlab function函数的用法如下：为function [输出变量] = 函数名称(输入变量）。function用来定义函数，一般一个函数放在一个.m文件里。举个简单的例子建立一个myfunction.m，然后在文件中写。function y=myfunction(a,b)其中a,b是输入函数的参数，y是函数返回的值。当需要返回多个值时，可以将y看作一个数组。或者直接将函数的开头写成如function [x,y]=myfunction(x,y)的形式。然后就是定义函数的内容，通俗的说就是怎样由输入参数a,b得到返回值y。比如最简单的y=a+b。所有函数的开头必须加上“function()”给函数命名：最好的命名方法是，函数拥有的名称和它所表现的功能相一致。名字可以包含字母或下划线(不可以包含数字)。加上一个“{”：函数的代码必须写在“{”符号之后。插入一段函数代码，加上一个“}”：函数书写完毕后，必须加上“}”符号。扩展资料：MATLAB用function定义一个函数后，怎么让这个函数生效是把这个M文件存到哪？你可以通过在命令栏中输入function，点击右键或者按F1来查询function函数的帮助文档。在帮助文档中对function函数进行了详细的功能描述，使用方法的介绍。function函数基本形式为：function[]=fun();即为：[输出]=fun(输入)。Function函数是包含在Function和End Function语句之间的一组VBScript语句. Function函数与Sub过程类似。但Function函数可以有返回值.可以使用参数.如果Function函数没有任何参数,则Function语句必须要包含空括号。

function用来定义函数，一般一个函数放在一个.m文件里。建立一个myfunction.m，然后在文件中写。function y=myfunction(a,b)其中a,b是输入函数的参数，y是函数返回的值。当需要返回多个值时，可以将y看作一个数组，或者直接将函数的开头写成如function [x,y]=myfunction(x,y)的形式。然后就是定义函数的内容，通俗的说就是怎样由输入参数a,b得到返回值y。比如最简单的y=a+b。function y=myfunction(a,b)y=a+b;保存为一个m文件myfunction.m，调用方式如下。num=myfunction(3,4);这样就可以由函数中的y=a+b得到num的值是3+4=7。在matlab中调用函数只能在.m文件中定义，并且文件名和函数名形同。

1、先通过在命令栏中输入function，点击右键或者按F1来查询function函数的帮助文档。然后在帮助文档中对function函数进行了详细的功能描述，使用方法的介绍。2、M函数除了直接用函数名调用之外，也可以进行参数传递，使得Matlab应用更加方便。M函数文件以function开头，格式为function输出变量=函数名称（输入变量）语句。3、function输出变量=函数名称（输入变量）函数内容end例如，利用matlab的程序编辑器，来编写函数文件。4、但Function函数可以有返回值。可以使用参数。如果Function函数没有任何参数，则Function语句必须要包含空括号。Function函数通过函数名返回一个值。

1.新建一个m文件，把function函数复制进去；然后保存（名称应该是自动生成的，即为sinh.m，不过最好不要用和系统中的文件一样的名称，可以先help查看一下matlab是否有和sinh重复的名称，若有的话，得把function函数改个名称），保存的地址为当前工作目录。2.在命令窗口调用function函数：sinh(5)

fi — Defines a fixed-point numeric object in the MATLAB workspace. Each fi object is composed of value data and a numerictype object.fimath — Governs how overloaded arithmetic operators work with fi objectsfipref — Defines the display, logging, and data type override preferences of fi objectsnumerictype — Defines the data type and scaling attributes of fi objectsquantizer — Quantizes data sets可以创建以上的对象，fi 定义一个定点对象 fimath 重载操作符 fipref 定义显示，记录数据类型重叠特性 numerictype 定义数据类型和尺度属性 quantizer 量化数据集

几乎没有一个男人看到这部电影_but_被英雄的精神。

validation checks是用来防止网络训练过适的，即是训练时网络的验证误差大于训练误差的次数，默认值为6就是大于次数达到6次训练就会停止。网络训练过适时，即使performance还在减小网络性能也不会提高，反而会降低，因此要设定过适次数来提前训...

dhku

请教一下专业人士吧

原因：embedded matlab function模块要进行代码生成，所以语法要求比一般的MATALB程序严格。你现在遇到的问题是因为if...elseif...end语句结构中缺少一个else分支，这样就存在隐患——如果前面的if和elseif条件都不满足，y就会没有赋值。解决方法（选择其一即可）：1、在if语句最后增加一个else分支，指定前面条件都不满足的时候y取什么值；2、在if语句之前增加一个赋值语句，指定y的默认值。

他真是一个不为所知的人，住在不为所知的地方，不为谁而作他不为所知的计划。还没有一点想法吗？不知道他将去哪。他难道不是有一点箱我们吗？不为所知的人，世界为你掌控。他盲目地只看他想看的。不为所知的人，你能看见我吗？还没有一点想法吗？不知道他将去哪。他难道不是有一点像我们吗？不为所知的人，不要担心。慢慢来，不要着急。顺其自然，直到有人向你伸手帮助，他真是不为所知的人，住在不为所知的地方，不为谁而作他不为所知的计划

Nowhere man 无处男子 Beatles 披头士 He"s a real nowhere Man,他是一个真正的无处男 Sitting in his Nowhere Land,坐在他无处土地 Making all his nowhere plans使他无处计划 for nobody.对谁都没有好处. Doesn"t have a point of view,没有一个角度来看, Knows not where he"s going to,也不知他的去, Isn"t he a bit like you and me?不是他有点像你和我? Nowhere Man, please listen,无处男子,请听, You don"t know what you"re missing,你不知道你错过, Nowhere Man, the world is at your command.nowhereman,世界就在你的指挥. He"s as blind as he can be,他是盲人,因为他可以, Just sees what he wants to see,只看到他所不愿看到的, Nowhere Man can you see me at all?无处男子你看我吗? Nowhere Man, don"t worry,nowhereman,别急, Take your time, don"t hurry,慢慢来,不急于一时, Leave it all till somebody else离开这一切直到别人 lends you a hand.lends你一臂之力. Doesn"t have a point of view,没有一个角度来看, Knows not where he"s going to,也不知他的去, Isn"t he a bit like you and me?不是他有点像你和我? Nowhere Man, please listen,无处男子,请听, You don"t know what you"re missing,你不知道你错过, Nowhere Man, the world is at your command.nowhereman,世界就在你的指挥. He"s a real Nowhere Man,他是一个真正的无处男 Sitting in his Nowhere Land,坐在他无处土地 Making all his nowhere plans使他无处计划 for nobody.对谁都没有好处. Making all his nowhere plans使他无处计划 for nobody.对谁都没有好处. Making all his nowhere plans使他无处计划 for nobody.对谁都没有好处.

matlab at frequencyMATLAB在频率matlab at frequencyMATLAB在频率

22

我在十二点吃午饭，现在是十二点了。

根据名字猜是开根吧

吃辣子中文吃吃吃吃吃，这个挺好的，去看看吧，特别好，特别有这个

"ideal_bs" 这个函数你没有写?还是写了没有放在同一个目录下? 顺便BS一下楼主下载人家的程序,还不下载完全...

这段 MATLAB 代码用来检测人脸。具体来说，它会执行以下操作：使用 webcam 函数打开摄像头，并获取一张图片，保存在变量 pic 中。使用 vision.CascadeObjectDetector 函数创建一个对象检测器，用于检测人脸。使用 imshow 函数显示图片。进入循环，每次都会获取一张新的图片，并将其转换为灰度图，保存在变量 pic2 中。使用 step 函数检测图片中的人脸，并将结果保存在变量 bbox 中。使用 insertObjectAnnotation 函数在图片中插入标注，表示检测到的人脸的位置。使用 imshow 函数显示图片。该代码将不断重复这些步骤，直到用户手动停止程序。

给你一个正确的程序，你自己参考一下吧！function [x,k]=Jacobimethod(A,b,x0,N,emg)% A:线性方程组左端矩阵% b:线性方程组右端向量% x0:迭代初值% N:迭代次数上界，若迭代次数大于n，则迭代失败% emg:精度指标% k:迭代次数% x:用迭代法求得的线性方程组的近似解n=length(A);x1=zeros(n,1);x2=zeros(n,1);x1=x0; k=0;r=max(abs(b-A*x1));while r>emgfor i=1:nsum=0;for j=1:nif i~=jsum=sum+A(i,j)*x1(j);endendx2(i)=(b(i)-sum)/A(i,i);endr=max(abs(x2-x1));x1=x2;k=k+1;if k>Ndisp("迭代失败，返回")；return;endend祝朋友好运！

WRF模式为完全可压缩以及非静力模式，为了兼顾计算效率和稳定性，水平方向采用Arakawa C(荒川C)网格点（重点考虑1-10km），垂直方向则采用地形跟随质量坐标。Arakawa-C网格点为交错格点，以下为图示说明： 因此，在以下U、V变量定义中

meshgrid用来生成x-y平面上的小矩形顶点坐标值的矩阵，也称为格点矩阵.meshgrid也适用于三元函数u=f(x,y,z)[X,Y]=meshgrid(x,y):绘制二维图形时生成小矩形的格点[X,Y]=meshgrid(x):等价于[X,Y]=meshgrid(x,x)[X,Y,Z]=meshgrid(x,y,z):绘制三维图形时生成空间曲面的格点[X,Y]=meshgrid(x):等价于[X,Y,Z]=meshgrid(x,x,x)

和楼主类似的问题。。。mark一下我也是想避免在给structure赋值时使用循环，但我的错误提示和楼主的不太一样。。我尝试的时候也出现过你那个错误提示。。楼主提问这么长时间了应该明白了吧= =想请教一下。。我的问题http://zhidao.baidu.com/question/1510263486424463300.html?sort=6#answer-1785853173

目测你这个函数定义有误Function Eliminate%逐点剔除法建议这样定义：Function []=Eliminate()"Error: Missing variable or function"说明你函数定义是缺自变量，即便不需要自变量也要用空括号表示即Eliminate()

桌面上鼠标右击，选择属性，看你的图片的存放位置，然后写入imread里。例如：I=imread("E:斑块(中文)截图白斑3.tif");存放位置+图片名字+.图片格式

相关图片 1969年，在美国，新停战动员委员会(New Mobilization Committee to End the War)组织了一场全国性的三天反战示威，11月5日在华盛顿特区的盛大游行是整个活动的最高潮。 此主题相关图片 2004年10月的，美国的联邦法院法官，日前接纳历史学家Jonathan Wiener申请，判FBI将十页记载了约翰连侬死前政治活动的机密文件解禁。 自约翰在80年被枪後，Jonathan分别在84年及2000年推出两本有关约翰的书，而第二本就是根据他自83年起以「资讯自由」为由控告美国政府，其後在97年解禁的二百四十八页所参考写成的。 被联联法官下令解禁的十页文件，记载了70年代约翰参与反越示威的经过，当中的其他内容，有些更是牵涉英国政府的。美国政府对法官的判决可在60日内提出上诉。 1998年蓝侬之子：美国政府，杀我爸 已故披头四合唱团主唱约翰蓝侬的儿子席恩蓝侬，向纽约客杂志表示，他相信美国政府是父亲遭枪杀的幕後元凶，二十二岁的席恩指出，由於父亲是一位反文化的革命份子，对当时保守的美国政府而言，是极其危险的，如果他下令明日炸白宫的话，很可能有一万歌迷会真的付诸行动。 2000.02.19 路透社 未公开的蓝侬档案，似与美国国家安全有关 轰动一时的摇滚乐团披头四成员之一约翰．蓝侬（John Lennon），在遭疯狂歌迷暗杀後20年，1名历史系教授调阅相关资料，但目前尚有10份机密文件深锁於档案柜，因美国司法单位认为一旦公布，恐将危及国家安全。 据路透社报导，加州大学历史系教授韦纳申请调阅蓝侬的3封私人信件，为的是作为和美国政府打官司的证物，他真正的目的是要求美国联邦调查局（FBI）让他调阅10份有关蓝侬的机密文件。但美国司法部检察官卡巴雷罗18日警告美国政府，若公布蓝侬支持左派政治活动的资料，可能危及国家安全。 其实早在多年前，韦纳就已积极申请调阅蓝侬相关资料，在几番诉讼波折後，美国联邦调查局（FBI）终在1997年与韦纳和解，公布所有蓝侬相关文件，只除了目前正在打官司的10份极机密文件。因 FBI 表示，有1外国政府要求文件不得公开於世。报导指出，这个所谓的「外国政府」据信是英国政府。 韦纳相信，美国政府千方百计想保密的机密文件，内容与1960年代蓝侬资助左派组织的活动有关。 判决此司法诉讼案的美国地方法院法官鲁宾斯，倾向作出对韦纳有利的判决。但卡巴雷罗在听证会後向韦纳的律师表示，美国政府仍会针对此案上诉，因公布10份蓝侬文件，可能危害国家安全。 韦纳闻讯後义愤填膺地说：「这不是国家全安的资料，只是1名已故摇滚巨星的政治活动纪录。」摇滚巨星蓝侬1980年12月8日在归返曼哈顿上西区住宅途中，遭一名疯狂歌迷暗杀身亡。 2000.02.21 明日报 「X档案」曝光 蓝侬死因再被怀疑 1980年蓝侬被刺案本已尘封多时，2000年却因凶手服刑届满、FBI、英国情报局秘密档案曝光而再掀风暴。 1980年约翰蓝侬被枪击死於自家门口的案件，本来已经成为乐迷尘封的伤心往事，不过2000年刺杀他的凶手将服刑届满出狱，再加上作家维纳(Wiener)的惊爆新书《给我真相：蓝侬FBI秘密档案》的发行问世，让英美媒体再度挖出两国政府情报局对蓝侬进行的众多违法监控事实，并再度让世人怀疑蓝侬的真正死因 综合路透社、美联社报导，作家维纳的新书《给我真相：蓝侬FBI秘密档案》中披露了100页FBI隐藏多年的违法监控蓝侬报告，其中内容包括美国政府如何监控、阻挠蓝侬捐款或帮助左派共党组织与少数族裔组织。 无独有偶地，21日英国政府违法秘密监控、拦截蓝侬捐给「爱尔兰共和军」(IRA)的捐款与实质协助，「披头四」传记作者杭特(Hunter Davies)获悉後表示：「我一点都不意外，我所认识的他就是这样一个慷慨而大无畏支持自由主义的人。」。 「You can say I"m a dreamer. But I"m not the only one....」，约翰蓝侬是个梦想家，而全球认同他的乐迷如今有增无减，英美秘密档案曝光让蓝侬的死因又成为2000年最热门的话题之一，英美两国的人权组织都认为这个案子有必要公布更多的真相。 2000.02.22 明日报 洋子否认蓝侬援助爱尔兰共和军 蓝侬之死再度引起注意，蓝侬遗孀洋子出面否认蓝侬援助爱尔兰共和军。 综合路透社、美联社报导，作家维纳的新书《给我真相：蓝侬FBI秘密档案》中披露了100页FBI隐藏多年的违法监控蓝侬报告，21日英国媒体报导英国政府也违法秘密监控、拦截蓝侬援助「爱尔兰共和军」(IRA)，不过约翰蓝侬的遗孀小野洋子(Yoko Ono)台北时间22日终於忍不住出面澄清：「我丈夫从未援助爱尔兰共和军，他只帮助需要帮助的人。」。 这位曾为了抗议战争、主张和平而与蓝侬公开在床上裸体拥抱面对媒体的「前英伦摇滚第一夫人」说：「即使在六○年代，(IRA)暴力的手段也不为人支持。」，她并表示蓝侬的金钱援助向来是以孤儿、饥民或身处困境的老弱妇孺为主要对象。 不过《给我真相：蓝侬FBI秘密档案》的作者加州大学柏克莱分校教授维纳仍然坚持：「有证据显示蓝侬的确曾对爱尔兰共和军表示同情。」。 2000.02.22 外电报导 小野洋子：蓝侬从未资助北爱共和军 英国媒体指蓝侬曾金援北爱共和军，他的遗孀小野洋子出面否认。 据路透社报导，前披头四合唱团成员约翰蓝侬的遗孀小野洋子，21日出面反驳英国媒体有关蓝侬曾资助爱尔兰共和军(IRA)的报导。 据英国观察家周报20日报导，前英国间谍雪勒表示，曾经看过跟蓝侬有关的机密文件，里面提到蓝侬於70年代初期捐钱给IRA。 小野洋子告诉西班牙媒体，她和蓝侬曾经捐钱给许多国家，包括爱尔兰在内，但那些钱是为了协助「需要帮助的儿童、孤儿及妇女」。她说：「早在60年代，我就不相信暴力是做事情的方法。」小野洋子的艺术作品目前正在西班牙北部城市萨拉哥萨展览。 蓝侬生前曾在英国一场集会中高举「IRA对抗英国帝国主义胜利」的标语，他支持IRA终结英国对北爱尔兰统治的目标，众所皆知。但他的遗孀小野坚持：「我先生没有拿钱给IRA。」 美国联邦调查局的档案中，有一些与蓝侬有关的机密文件。研究蓝侬的加州大学历史学教授怀纳则正在努力争取其中10件档案，最近法院便判决3封机密信函应交给怀纳。怀纳向英国广播公司BBC表示，尽管蓝侬对爱尔兰的民权运动表示同情，但就蓝侬与IRA有关联一事，英国情报单位恐怕有误。

歌名是Hello,Goodbye吧you say yes, i say no.you say stop and i say go go go, oh no.you say goodbye and i say hellohello helloi don"t know why you say goodbye, isay hellohello helloi don"t know why you say goodbye, i say hello.i say high, you say low.you say why and i say i don"t know, oh no.you say goodbye and i say hellohello goodbye hello goodbye hello hellohello goodbye i don"t know why you say goodbye, i say hellohello goodbye hello goodbye hello hellohello goodbye i don"t know why you say goodbyehello goodbye i say goodbye.why why why why why why do you say goodbye goodbye, oh no?you say goodbye and i say hellohello helloi don"t know why you say goodbye, i say hellohello helloi don"t know why you say goodbye, i say hello.you say yes i say yes i say no but i may mean no.you say stop i can stay and i say go go go till it"s time to go oh, oh no.you say goodbye and i say hellohello helloi don"t know why you say goodbye, i say hellohello helloi don"t know why you say goodbye, i say hellohello helloi don"t know why you say goodbye, i say hello hello.hela heba helloa cha cha.................

在一个subsystem中放置一个goto，goto与一个经过运算后的输入信号连接；对goto中的tag命名，在tag visibility选择local/scope/global类型。在另外一个subsystem或与goto同一个subsystem中放置from，在from中的goto tag处输入goto相同的名称，点击update tags。

( fName, "reader", [cache] )是创建用于读取seq文件的接口sr当我在matlab中以sr = seqIo( "V001.seq", "r", 2048 )

如何利用MATLAB提取datasheet中曲线的数据1、首先设置编程得到的曲线，并赋值给变量h。1h=findobj(gca,"type","line");2、得到变量h的横坐标。1x=get(h,"xdata");3、得到变量h的纵坐标。1y=get(h,"ydata");4、变量h的横纵坐标，就是编程得到曲线对应的数据。

首先把变量名改成 x1 x2 ，并且确保之前 syms 过，即在之前写上：syms x1 x2;也就是 f 的表达式变成：f = (x1 + x2) ^ 4 + x2 ^ 2然后，梯度和海森矩阵分别为：g = jacobian(f);h = jacobian(g);

>>a=imread("D:Program Files oolboximagesimdemos ice.png");>> b=gradient(double(a));转化图像为double型，即可，gradient中h(:,ones(p,1))估计是double，我转化后可以运行，得到的是水平方向的梯度。>>[ax,ay]=gradient(double(a));可以得到水平和垂直方向的梯度。

./和/有什莫区别？不是对矩阵里每个元素都除吧

I和X Y的灰度值不在一个范围内，分开显示看看吧，就是显示在2个图里

求助matlab中gradient的用法[x,y]=meshgrid([-5:0.5:5])z=1./(x.^2-2*x+4)+1./(y.^3-2*y+4)g=gradient(z)subplot(1,2,1),surf(x,y,z)subplot(1,2,2),surf(x,y,z,g)

Gradient是梯度的意思，BP神经网络训练的时候涉及到梯度下降法，表示为梯度下降的程度与训练过程迭代次数（步长）的关系。Performance是神经网络传递误差大小的意思，表示为均方差与训练过程迭代次数（步长）的关系。

Maple是目前世界上最为通用的数学和工程计算软件之一，在数学和科学领域享有盛誉，有“数学家的软件”之称。Maple 在全球拥有数百万用户，被广泛地应用于科学、工程和教育等领域，用户渗透超过96％的世界主要高校和研究所，超过81％的世界财富五百强企业。Maple系统内置高级技术解决建模和仿真中的数学问题，包括世界上最强大的的符号计算、无限精度数值计算、创新的互联网连接、强大的4GL语言等，内置超过5000个计算命令，数学和分析功能覆盖几乎所有的数学分支，如微积分、微分方程、特殊函数、线性代数、图像声音处理、统计、动力系统、等。 　　Maple不仅仅提供编程工具，更重要的是提供数学知识。Maple是教授、研究员、科学家、工程师、学生们必备的科学计算工具，从简单的数字计算到高度复杂的非线性问题，Maple都可以帮助您快速、高效地解决问题。用户通过Maple产品可以在单一的环境中完成多领域物理系统建模和仿真、符号计算、数值计算、程序设计、技术文件、报告演示、算法开发、外部程序连接等功能，满足各个层次用户的需要，从高中学生到高级研究人员。Maple 主要技术特征1. 强大的求解器：数学和分析软件的领导者2. 技术文件环境：重新定义数学的使用性3. 知识捕捉：不仅是工具，更是知识4. 外部程序连接 MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。 基本功能MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连 matlab开发工作界面接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 　　MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且mathwork也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++ ，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。优势（1）友好的工作平台和编程环境（2）简单易用的程序语言（3）强大的科学计算机数据处理能力 （5）应用广泛的模块集合工具箱（6）实用的程序接口和发布平台（7）应用软件开发（包括用户界面）

价格都贵到你不能承受，即使是一些大企业，也不能承受。。。不过中国有盗版。。。而且为了学术研究中国说了，那不算盗版。。。据我所知matlab是按年卖使用权的，怎么也要几万到几十万美元/每年

第一种情况：选择 Maple 或者 MuPAD Engine 引擎你可以使用Maple engine 代替默认的MuPAD engine。首先你需要一个匹配版本的Maple，例如Maple 12.02版本支持Matlab 2008b及以前版本。1、在Matlab命令行输入 symengine2、显示图形窗口3、选择Maple engine，点击 Maple 按钮，然后点击 Maple Location 按钮，找到你机器上Maple的安装位置。第二种情况：机器(Windows)上安装了Maple Toolbox for Matlab如果你希望Maple engine称为默认的符号计算引擎，你需要设置环境变量。本帖隐藏的内容需要回复才可以浏览在Windows中环境变量设置如下：鼠标右击“我的电脑”——>选择“属性(R)”——>“高级”标签栏——>点击底部的“环境变量(N)”。将MATLAB_SYMBOLIC变量值设为“Maple”，如果没有则创建一个新的系统变量。如果希望使用MuPAD，则将“MATLAB_SYMBOLIC”的值改为“Matlab”。还有第一步：并安装Maple；第二步：安装Maple Toobox for Matlab 13，安装到你的maple安装目录下；第三步：拷贝MapleToolboxForMatlab.dat 到 ～Maple 13 oolboxMapleToolboxForMatlablicense ("~"指maple安装根目录)。第四步：运行matlab2010a，cmd窗口执行`toolbox_version`。会返回结果toolbox_version=maple;上面说明maple工具箱安装成功，但是不经设置的话matlab仍然使用内置的MuPAD, 即便是toolbox_version=maple，并不是说你现在的符号运算都是使用maple内核，其实现在的使用还是matlab默认的MuPad。现在我们开始设置最关键的一个东西，也是让maple内核工作的关键！第五步：在MATLAB命令提示行里面输入pathtool，这是会弹出一个对话框，下拉到最后，选中最后两个，然后点击左边的”Move to Top“（移动到顶端）按钮，然后选择保存。exit退出matlab，现在设置好了。重新启动matlab，现在可以用maple内核了点击“确认”按钮，关闭对话框。重启Matlab，通过命令确认你选择了何种符号计算引擎。1） MuPAD 的情况： >> evalin(symengine,"polylib::discrim(a*x^2+b*x+c,x)")2） Maple 的情况：1. >> maple("discrim(a*x^2+b*x+c,x)")

在命令栏输入 help bias就知道了

小题1:A小题2:C小题3:B 试题分析：文章主要描写了大西洋惊人的力量。与内陆海相比，大西洋似乎有着永不枯竭的力量。大西洋就像一个有呼吸、有心跳、有心理、有性格的人一样。小题1:根据“Unlike real inland seas, which seem strangely still, the Atlantic Ocean is rich in oceanic liveliness all the time. ...the Atlantic Ocean is never dead and dull ...It is an ocean that moves, impressively and endlessly.”可知，真正的内陆海都是出奇的平静，而大西洋却一直都是充满活力，故选A。小题2:根据“As the poem suggests, the Atlantic Ocean is never dead and dull. It is an ocean that moves, impressively and endlessly. It makes all kinds of noise---it is forever thundering, boiling, crashing and whistling.”可知，本段是对大西洋无穷力量的描述，故选C。小题3:根据“It is easy to imagine the Atlantic Ocean trying to draw breath ...giving to the waters a feeling of heartbeat ...And it has a psychology. It has personalities”可知，作者把大西洋比作了一个有血有肉的人，故选B。

方法/步骤打开Simulink仿真环境。搭建Simulink仿真模型：根据原理搭建仿真模型，并检查仿真模型的错误纠正之。模块封装：选择需要封装的部分，单击鼠标右键选择“create subsystem”选项即可封装模块。模块封装外部接口：模块封装好之后，可以检查外部接口是否合适，主要检查外部接口顺序、模块大小等。尽量使之美观、简洁。外部接口顺序调整：如果需要调整外部接口顺序，可以打开封装好的模块，调整图中所示的接口顺序即可实现外部顺序的调整。外部美化：一个好的Simulink仿真图形会令人有兴趣有心情看下去，外部的美化也需要注意。步骤阅读

如果是系统自带的，你可以直接用，如果是外部的或者是自编的你需要先把文件夹拷贝到tools文件夹下，再设置路径。x0dx0aMatlab常用工具箱介绍(英汉对照)x0dx0aMatlab Main Toolbox——matlab主工具箱x0dx0aControl System Toolbox——控制系统工具箱x0dx0aCommunication Toolbox——通讯工具箱x0dx0aFinancial Toolbox——财政金融工具箱x0dx0aSystem Identification Toolbox——系统辨识工具箱x0dx0aFuzzy Logic Toolbox——模糊逻辑工具箱x0dx0aHigher-Order Spectral Analysis Toolbox——高阶谱分析工具箱x0dx0aImage Processing Toolbox——图象处理工具箱x0dx0aLMI Control Toolbox——线性矩阵不等式工具箱x0dx0aModel predictive Control Toolbox——模型预测控制工具箱x0dx0aμ-Analysis and Synthesis Toolbox——μ分析工具箱x0dx0aNeural Network Toolbox——神经网络工具箱x0dx0aOptimization Toolbox——优化工具箱x0dx0aPartial Differential Toolbox——偏微分方程工具箱x0dx0aRobust Control Toolbox——鲁棒控制工具箱x0dx0aSignal Processing Toolbox——信号处理工具箱x0dx0aSpline Toolbox——样条工具箱x0dx0aStatistics Toolbox——统计工具箱x0dx0aSymbolic Math Toolbox——符号数学工具箱x0dx0aSimulink Toolbox——动态仿真工具箱x0dx0aSystem Identification Toolbox——系统辨识工具箱x0dx0aWavele Toolbox——小波工具箱x0dx0ax0dx0a例如：控制系统工具箱包含如下功能：x0dx0a连续系统设计和离散系统设计x0dx0a状态空间和传递函数以及模型转换x0dx0a时域响应（脉冲响应、阶跃响应、斜坡响应）x0dx0a频域响应（Bode图、Nyquist图）x0dx0a根轨迹、极点配置x0dx0ax0dx0a较为常见的matlab控制箱有：x0dx0ax0dx0a控制类：x0dx0ax0dx0a控制系统工具箱(control systems toolbox)x0dx0a系统识别工具箱（system identification toolbox）x0dx0a鲁棒控制工具箱（robust control toolbox）x0dx0a神经网络工具箱（neural network toolbox）x0dx0a频域系统识别工具箱（frequency domain system identification toolbox）x0dx0a模型预测控制工具箱（model predictive control toolbox）x0dx0a多变量频率设计工具箱（multivariable frequency design toolbox）x0dx0ax0dx0a信号处理类：x0dx0a信号处理工具箱（signal processing toolbox）x0dx0a滤波器设计工具箱（filter design toolbox）x0dx0a通信工具箱（communication toolbox）x0dx0a小波分析工具箱（wavelet toolbox）x0dx0a高阶谱分析工具箱（higher order spectral analysis toolbox）x0dx0ax0dx0a其它工具箱：x0dx0a统计工具箱(statistics toolbox)x0dx0a数学符号工具箱（symbolic math toolbox）x0dx0a定点工具箱（fixed-point toolbox）x0dx0a射频工具箱（RF toolbox）x0dx0ax0dx0a1990年，MathWorks软件公司为Matlab提供了新的控制系统模型化图形输入与仿真工具，并命名为Simulab，使得仿真软件进入了模型化图形组态阶段，1992年正式命名为Simulink，即simu（仿真）和link（连接）。matlab7.0里的simulink为6.0版本，matlab6.5里的simulink为5.0版本。x0dx0ax0dx0aMATLAB的SIMULINK子库是一个建模、分析各种物理和数学系统的软件，它用框图表示系统的各个环节，用带方向的连线表示各环节的输入输出关系。x0dx0a　　启动SIMULINK十分容易，只需在MATLAB的命令窗口键入“SIMULINK”命令，此时出现一个SIMULINK窗口，包含七个模型库，分别是信号源库、输出库、离散系统库、线性系统库、非线性系统库及扩展系统库。x0dx0a　　1.信号源库x0dx0a　　包括阶跃信号、正弦波、白噪声、时钟、常值、文件、信号发生器等各种信号源，其中信号发生器可产生正弦波、方波、锯齿波、随机信号等波形。x0dx0a　　2.输出库x0dx0a　　包括示波器仿真窗口、MATLAB工作区、文件等形式的输出。x0dx0a　　3.离散系统库x0dx0a　　包括五种标准模式：延迟，零-极点，滤波器，离散传递函数，离散状态空间。x0dx0a　　4.线性系统库x0dx0a　　提供七种标准模式：加法器、比例环节、积分环节、微分环节、传递函数、零-极点、状态空间。x0dx0a　　5.非线性系统库x0dx0a　　提供十三种常用标准模式：绝对值、乘法、函数、回环特性、死区特性、斜率、继电器特性、饱和特性、开关特性等。x0dx0a　　6.系统连接库　　包括输入、输出、多路转换等模块，用于连接其他模块。x0dx0a　　7.系统扩展库x0dx0a　　考虑到系统的复杂性，SIMULINK另提供十二种类型的扩展系统库，每一种又有多种模型供选择。x0dx0a　　使用时只要从各子库中取出模型，定义好模型参数，将各模型连接起来，然后设置系统参数，如仿真时间、仿真步长、计算方法等。SIMULINK提供了Euler、RungeKutta、Gear、Adams及专用于线性系统的LinSim算法，用户根据仿真要求选择适当的算法。x0dx0ax0dx0a当然，不同版本的Matlab/Simulink内容有所不同。x0dx0ax0dx0a另外，Simulink还提供了诸如航空航天、CDMA、DSP、机械、电力系统等专业模块库，给快速建模提供了很大的便利。

下载一个工具箱，添加到matlab中，直接运用：在matlab的file下面的set path把它加上，把路径加进去后在file→Preferences→General的Toolbox Path Caching里点击update Toolbox Path Cache更新一下，就OK了1.1 如果是Matlab安装光盘上的工具箱，重新执行安装程序，选中即可；1.2 如果是单独下载的工具箱，一般情况下仅需要把新的工具箱解压到某个目录。2 在matlab的file下面的set path把它加上。3 把路径加进去后在file→Preferences→General的Toolbox Path Caching里点击update Toolbox Path Cache更新一下。4 用which newtoolbox_command.m来检验是否可以访问。如果能够显示新设置的路径，则表明该工具箱可以使用了。或者将工具箱复制到toolbox目录然后在set path里面添加这个目录就可以用了。在matlab的file下面的set path把它加上，把路径加进去后在file→Preferences→General的Toolbox Path Caching里点击update Toolbox Path Cache更新一下，就OK了以下是我在别的地方看到的，转过来你参考一下吧。首先说说添加到matlab搜索路径好处：1 对 n——你只需要存储一个副本，就可以在其他地方使用。具体来说，假设你在数据盘D上新建了两个目录abc和def，这两个工程（每个目录下的所有程序相应地称为一个工程）都需要调用同一个（些）函数（简称工具箱），这时候，如果你没有把该工具箱添加到matlab的搜索路径下，则需要分别把工具箱中所有用到的文件都复制到目录abc和def下才能正确运行。这显然浪费空间，所以，matlab提供了一个搜索路径（默认在matlab安装目录下的toolbox中），只要把工具箱对应的整个文件夹复制到搜索路径对应的目录下，并且通知matlab一声（把该路径正确添加到搜索路径中），就可以在abc和def中使用这个工具箱了（即无论你的工程文件在哪个目录（有效的目录）下都可以访问这个工具箱中的函数）。下面就以matlab安装目录下的toolbox目录作为默认的添加路径进行详细说明。1. 如何添加工具箱？以下是添加工具箱的方法：（论坛很多人转载过，这里就不作区分了，下面的举例也一样）如果是Matlab安装光盘上的工具箱，重新执行安装程序，选中即可。如果是单独下载的工具箱，则需要把新的工具箱（以下假设工具箱名字为svm）解压到toolbox目录下，然后用addpath或者pathtool把该工具箱的路径添加到matlab的搜索路径中，最后用whichnewtoolbox_command.m来检验是否可以访问。如果能够显示新设置的路径，则表明该工具箱可以使用了。具体请看工具箱自己代的README文件。1.1 举例：要添加的工具箱为svm，则解压后，里边有一个目录svm，假设matlab安装在D:MATLAB6p5，将svm目录拷贝至D:MATLAB6p5 oolbox，然后运行matlab，在命令窗口输入addpathD:MATLAB6p5 oolboxsvm回车，来添加路径。然后在svm目录下，任意找一个m文件，以svcinfo.m为例，在命令窗口中输入which svcinfo.m。如果显示出该文件路径，如 D:MATLAB6p5 oolboxsvmsvcinfo.m，则安装成功，当然也可以在命令窗口输入path来查看。上面的说明和例子基本上介绍了在matlab中如何添加工具箱，下面是其他补充：1.2 添加方式总结：事实上，有两种添加工具箱到matlab搜索路径的方法：其一是用代码，其二是用界面。其实无论用哪种方法，都是修改pathdef.m这个文件，阁下如果是高手，可以直接打开该文件修改，呵呵，对此这里不作讨论。1.2.1 代码方式：1.2.1.1 适用于添加下载的工具箱（别人的）：在命令窗口输入addpath D:MATLAB6p5 oolboxsvm 或者addpath(‘D:MATLAB6p5 oolboxsvm");但是这种方法只能添加svm目录，如果该目录下有其他子文件夹，并且运行时候“隐式”调用到这些子文件夹（例如假设svm目录下存在子文件夹matdata，该子文件夹下有logo.mat这个文件，且在m文件代码中使用了诸如 load logo 这样的句子，即没有显式给出logo.mat的具体路径，则称为“隐式”），则不能正确访问。因此，有必要在添加时使用以下语句把svm目录下所有文件夹都添加到搜索路径中：500){this.resized=true;this.style.width=500;}"> 程序代码addpath(genpath(D:MATLAB6p5 oolboxsvm));另外，如果只使用以上代码，则退出matlab后，新添加的路径不会被保存下来，下次重新启动matlab后又需要重新添加。可以用savepath来解决这个问题，即在命令窗口中使用 savepath 便可。不过保存路径操作只对matlab 2006a、2006b有效（应该对版本7.0以上的都有效，因没有测试，所以不敢肯定），对6.5版本的matlab无效（该版本没有savepath这个命令）。1.2.1.2 适用于添加自己的工具箱（即工具箱自己编写，然后希望别人下载后当运行主文件时自动把路径添加到matlab搜索路径中）：在主文件中加入如下代码；500){this.resized=true;this.style.width=500;}"> 程序代码sCurrPath = fileparts(mfilename("fullpath"));addpath(genpath(sCurrPath)); %如果该工具箱没有其他子目录，则可以不需要用genpathsavepath; %这句可根据各人需要自行选择，相信上一节的解释已经很清楚1.2.2 界面方式：首先把主界面调出来，有两种方式，一种是在命令窗口输入pathtool，另一种是在matlab菜单调用（File -> Set Path…）。打开主界面后，相信不用太多解释了。Add Folder… 或者 Add with subfolders… 就可以了，添加完毕后，如果需要保存该路径，则在close前save一下就可以了。这种方式的好处是如果路径有问题，matlab会马上报错，告诉用户添加失败。2 工具箱添加失败：2.1 work目录和toolbox目录问题：你当然可以单独把一个或多个文件（不含文件夹）放在work目录下来实现“1对n”，因为work目录是其中一个matlab默认的搜索路径，但显然不能包含文件夹，如果包含文件夹，则同样需要把该文件夹添加到搜索路径中。相反，你不能单独把一个或多个文件（不含文件夹）放在toolbox目录下来实现“1对n”，因为toolbox这个目录并非matlab的一个默认搜索路径，除非你把toolbox文件夹添加到搜索路径中。2.2 由于路径名称而导致添加失败的总结：2.2.1 路径存在空格：错误：addpath C:Program FilesMATLABR2006b oolboxfinity正确：addpath(‘C:Program FilesMATLABR2006b oolboxfinity") 或者使用界面方式添加说明：在6.5版本的matlab不推荐使用带空格的路径，因为matlab6.5的安装路径是不允许有空格的，因此不晓得带空格的工具箱路径会出现什么问题。7.0以上应该没有问题。界面方式则估计没有问题，有问题matlab会自动报错的。2.2.2 路径存在中文：用addpath和界面方式均可以成功添加，但不推荐，最好使用英文路径。2.2.3 路径存在“@”字符：均不成功，不要使用“@”字符小结：出现其他与添加相关的错误时，用标准的路径和文件名试试吧。3 正确添加了工具箱，但运行调用时出错：3.1 版本问题：对于工具箱检测到matlab版本不兼容的问题，先引用一下happy教授的回复：=====================happy================================这个是有用的，估计是作者考虑到随着matlab的升级会出现不兼容的问题，因此加了类似的设置error("unsupported matlab version. please send an email.");以提示大家。从这一设置上看4.x和5.x版本的是完全没有问题的，之后的版本可能会出现问题，是否真的会出现需要进行测试=========================================================我也遇到过这类问题（如3.3节提到的spm5，下载得到的文件如果遇到2006版本的matlab则无法成功运行），但不是matlab的直接回复说“版本不支持”，而是普通的语法出错。我的经验是根据错误提示调试一下，看看问题出在哪里，然后再对程序作相应修改。例如在我的问题中，由于不支持最新的版本，所以2006a和2006b版本分别用下面的语句进行检查：strcmp(version("-release"),"2006a") 和 strcmp(version("-release"),"2006b")。3.2 程序中其他语法错误：估计部分语法不兼容，例如7.0版本以上的matlab可以使用&&、||、@(x)等等之类的符号，但是在6.5版本下无效，这种错误也需要自己手动修改一下代码。3.3 工具箱中的函数重名问题：Happy教授也提到过，除了改名，没有太好的办法。但是如果改的地方不止一个，例如遇到该文件进行了多次的自我调用（一个典型例子是以switch和case语句进行区分不同的操作），或者其他文件存在对该文件的调用时，这是个很头痛的事情。根据我的经验，也没有找到太好的解决办法。以我所用的两个工具箱为例，它们的文件夹名字不同（一个是spm2，一个是spm5），但是里面的主m文件名字都一样，为spm.m，如果两个工具箱同处于matlab的搜索路径中，会导致其中一个工具箱失效，更别说两者之间通过切换来进行调用了。我一般的解决方法是安装两个matlab，例如一个是6.5版，一个是2006版，把两个工具箱分别添加到不同的matlab中，这样就可以启动不同版本的matlab使用不同版本的工具箱了，呵呵，算是一个补救的方法吧。3.4 找不到m文件的问题：如果你已经把工具箱正确添加到搜索路径下，这种情况一般不会发生。也就是说，通过“文件夹copy -> 运行matlab ->添加路径”步骤后，该文件夹下所有文件应该都可以访问到。但是，当你对该文件夹下的m文件更新（包括修改和新增）了以后，此时如果你不重启matlab，则可能会出现找不到m文件的问题（特别是该工具箱中的文件没有依赖关系，它们只是被放在一起方便调用，当你新增一个m文件到该文件夹下而不重启matlab的时候，会造成这一新增m文件访问失败）。解决这个问题的方法之一显然是重启matlab，如果你不想重启，也可用第二个方法：在命令窗口输入: rehash toolbox ，“强制”matlab刷新toolbox目录下的所有文件，这样就可以正确访问了。

要将用Matlab语言编写的函数文件编译成可独立执行的*.exe文件(即可脱离Matalab环境的执行程序)，首先要安装和配置好Matlab Compiler，一般来说，在安装Matlab时就已经安装了相应版本的Matlab Compiler。只是不同版本的Matlab，其编译器的使用方法有一定的差异，这一点要引起一定的注意。在确定安装好Matlab Compiler后，还需要对Compiler进行适当的配置，方法是在Matlab命令窗口输入：Mbuild –setup然后根据提示执行相应的操作，使用者可根据自己计算机中现有编译器的情况选择合适的编译器，如VC++ 6.0、VC++7.0、Bland C的编译器等，目前Matlab好象还不支持VC++8.0(我计算机安装的就是VC++2005，Matlab就无法识别)。当然，如果你的计算机里根本就没有安装其他任何语言的编译器，也可选择Matlab自带的Lcc编译器，其实这个编译器对大多数用户已经够用了(我就是选择的Matlab自带的Lcc编译器)。配置好编译器后，自然就是对自己编写的M文件进行编译了。将M文件编译为独立可执行文件的语法是：>>mcc –m fun1.m fun2.m…..其中fun1就是最后的可执行文件的名称。另外，也可通过采用命令开关-o指定编译最终目标文件的名称，如 mcc –m main.m –o mrank_main，就是将编译后的文件指定为mrank_main.exe。编译后的生成文件根据编译器的版本不同而不同。具体的可参阅相关资料。如果要在没有安装matlab的计算机上执行编译后的程序，首先要将MATLAB701 oolboxcompilerdeploywin32中的MCRinstaller.exe安装到该计算机上(7.0以前的版本是mglinstaller.exe)。其次是要将“MCRinstaller.exe安装目录 untimewin32”这个路径添加到该计算机的环境变量中，添加的方法是：右击“我的电脑”“属性”“高级”“环境变量”“添加”指定一个变量名，然后将上述路径复制到里面就可以了。第三步是将编译生成的相相关文件拷贝到同一目录下(当然其他目录亦可)。第四步是打开MS-dos操作窗口，进入到编译后的*.exe程序所在的目录，执行编译生成的*.exe文件即可。需要说明的一个问题：如果你的程序还附带有图片(如Version上的图标)、Web页面(如help文档)，在编译的时候可能无法与M文件一起编译(因为我没有试过，也没有看到有关这方面的介绍)，这没有关系，先在计算机上执行一次你编译后生成的独立可执行文件，这时在同一目录下会生成一个以你编译后的程序名+_mcr结尾的文件夹，这时，你只需要把与程序相关的图标和Web页面拷贝到该文件夹中的Matlabprogram下面的一个子目录(一般与用于编译的文件夹名相同)中即可。其他的方法还有：(以下是网上搜索到的方法，没有试验过，有兴趣的朋友可以试一试)方法一：在matlab下的workspace里打comtool，点file-open project将我们先前建好的comtest.cbl工程文件打开，再点component--package component就实现了打包，此时到comtestdistrib文件夹里看，生成的comtest.exe就是打包后的解压程序，双击它会解压出一些文件，再点击解压出来的_install.bat就可以实现安装 方法二：如何将gui生成exe ？已有gui.m文件和gui.fig文件1 在matlab的command窗口中输入 mcc -B sgl GUI.m 2.将上步生成的文件包括*.m 文件和*.fig文件一起考到待运行的机器 此时仍需matlab所必需的动态连接库。3. 将 <matlab path >/extern/lib/win32/mglinstallar.exel拷贝到到待运行机器上4.在机器上先运行mglinstallar.exe, 然后选择解压目录，将在指定目录下解压缩出bin和toolbox两子目录，其中在binwin32目录下就是数学库和图形库脱离MATLAB运行所需的所有动态连接库，共有37个。可以将这些.dll考入system32, 也可以直接放在应用程序目录下。而toolbox目录则必须与应用程序同一目录。5.大功告成。 方法三：运行mcc -B sgl GUI.m 生成exe文件，把GUI.exe文件和GUI.ctf文件拷贝到目标电脑上，然后在目标电脑上安装 MCRInstaller.exe ,安装完以后，我点击GUI.exe,它显示：this application has requested the Runtime to terminate it in an unusual way。这个可能是哪里出错了呢？Solution:The segmentation violation can be caused due to a corruption in either the profile of the current user or by corrupted Windows libraries. In order to determine which is the cause, try logging in as a different username and run MATLAB. If you are able to successfully launch MATLAB under a different username:1. Log onto the machine, as the user that is receiving the error.2. On your desktop, click on Start > Run3. Type "regedit" in the Open text field and click OK, this will bring up the Registry Editor4. In the Registry Editor expand the following folders: HKEY_CURRENT_USER > Software > MathWorks > Matlab5. Under the MATLAB key, find the folder having the same name as the version number of the MATLAB that is crashing; and delete it. This Registry key will be regenerated the next time MATLAB is started. Click on the link below to view an image of the Registry Editor, which shows how to delete a registry key for MATLAB 6.5 (R13).After following these steps, you should be able to launch MATLAB successfully. If MATLAB crashes regardless of the user who is logged in, the crash may be caused by corrupted system libraries. For Windows NT, 2000, or XP, some customers have found that reinstalling Windows Service Packs have resolved the issue这个是启动6.5时出类似错误的解决方法，你试试7.0可否成功吧

yuir ney

function p = TSPmodify(G)n = size(G,1);p = [1:n 1];flag = true;while flag for i = 1:n-1 for j = i+1:n if G(p(i),p(j)) + G(p(i+1),p(j+1)) < G(p(i),p(i+1)) + G(p(j),p(j+1)) flag = false; p(i+1:j) = fliplr(p(i+1:j)); end end endend

不好意思，没看到维纳滤波，程序就不删了 %谱减法语音增强 %输入参数s 语音数据，fs 采样频率，p 下面有说明，共11个，可不输入，有默认值 %“过度减法(oversubtraction)”作减法的时候,保留一小部分原来的背景噪音,用这部分背景噪音来掩盖住音乐噪音的谱峰,从而消除了令人不悦的音乐噪音。%通过给的参数p，估计噪音，做谱减法。从而消除噪音。 function [ss,po]=specsubm(s,fs,p) %利用频谱相减（spectral subtraction）增强 [SS,PO]=(S,FS,P) % % implementation of spectral subtraction algorithm by R Martin (rather slow) % algorithm parameters: t* in seconds, f* in Hz, k* dimensionless % 1: tg = smoothing time constant for signal power estimate (0.04): high=reverberant, low=musical % 2: ta = smoothing time constant for signal power estimate %used in noise estimation (0.1) % 3: tw = fft window length (will be rounded up to 2^nw samples) % 4: tm = length of minimum filter (1.5): high=slow response to noise increase, low=distortion % 5: to = time constant for oversubtraction factor (0.08) % 6: fo = oversubtraction corner frequency (800): high=distortion, low=musical % 7: km = number of minimisation buffers to use (4): high=waste memory, low=noise modulation % 8: ks = oversampling constant (4) % 9: kn = noise estimate compensation (1.5) % 10:kf = subtraction floor (0.02): high=noisy, low=musical % 11:ko = oversubtraction scale factor (4): high=distortion, low=musical%检查函数的输入参数，如果输入少于三个，po为默认值，po的参数上面有说明 if nargin<3 po=[0.04 0.1 0.032 1.5 0.08 400 4 4 1.5 0.02 4]."; else po=p; endns=length(s); ts=1/fs; ss=zeros(ns,1);ni=pow2(nextpow2(fs*po(3)/po(8))); ti=ni/fs; nw=ni*po(8); nf=1+floor((ns-nw)/ni); nm=ceil(fs*po(4)/(ni*po(7)));win=0.5*hamming(nw+1)/1.08;win(end)=[]; zg=exp(-ti/po(1)); za=exp(-ti/po(2)); zo=exp(-ti/po(5));px=zeros(1+nw/2,1); pxn=px; os=px; mb=ones(1+nw/2,po(7))*nw/2; im=0; osf=po(11)*(1+(0:nw/2)."*fs/(nw*po(6))).^(-1);imidx=[13 21]"; x2im=zeros(length(imidx),nf); osim=x2im; pnim=x2im; pxnim=x2im; qim=x2im;for is=1:nf idx=(1:nw)+(is-1)*ni; x=rfft(s(idx).*win); x2=x.*conj(x);pxn=za*pxn+(1-za)*x2; im=rem(im+1,nm); if im mb(:,1)=min(mb(:,1),pxn); else mb=[pxn,mb(:,1:po(7)-1)]; end pn=po(9)*min(mb,[],2); %os= oversubtraction factor os=zo*os+(1-zo)*(1+osf.*pn./(pn+pxn));px=zg*px+(1-zg)*x2; q=max(po(10)*sqrt(pn./x2),1-sqrt(os.*pn./px));ss(idx)=ss(idx)+irfft(x.*q); end if nargout==0 soundsc([s; ss],fs); end

原来想注释完，有事要出去，基本的思想已写完，后面的看程序吧。不外乎是，分帧，加窗，估计噪音，傅里叶变换，反傅里叶变换等。%输入参数s 语音数据，fs 采样频率，p 下面有说明，共11个，可不输入，有默认值%“过度减法(oversubtraction)”作减法的时候,保留一小部分原来的背景噪音,用这部分背景噪音来掩盖住音乐噪音的谱峰,从而消除了令人不悦的音乐噪音。%通过给的参数p，估计噪音，做谱减法。从而消除噪音。function [ss,po]=specsubm(s,fs,p)%利用频谱相减（spectral subtraction）增强 [SS,PO]=(S,FS,P)%% implementation of spectral subtraction algorithm by R Martin (rather slow)% algorithm parameters: t* in seconds, f* in Hz, k* dimensionless% 1: tg = smoothing time constant for signal power estimate (0.04): high=reverberant, low=musical% 2: ta = smoothing time constant for signal power estimate% used in noise estimation (0.1)% 3: tw = fft window length (will be rounded up to 2^nw samples)% 4: tm = length of minimum filter (1.5): high=slow response to noise increase, low=distortion% 5: to = time constant for oversubtraction factor (0.08)% 6: fo = oversubtraction corner frequency (800): high=distortion, low=musical% 7: km = number of minimisation buffers to use (4): high=waste memory, low=noise modulation% 8: ks = oversampling constant (4)% 9: kn = noise estimate compensation (1.5)% 10:kf = subtraction floor (0.02): high=noisy, low=musical% 11:ko = oversubtraction scale factor (4): high=distortion, low=musical%检查函数的输入参数，如果输入少于三个，po为默认值，po的参数上面有说明if nargin<3 po=[0.04 0.1 0.032 1.5 0.08 400 4 4 1.5 0.02 4]."; else po=p; endns=length(s);ts=1/fs;ss=zeros(ns,1);ni=pow2(nextpow2(fs*po(3)/po(8)));ti=ni/fs;nw=ni*po(8);nf=1+floor((ns-nw)/ni);nm=ceil(fs*po(4)/(ni*po(7)));win=0.5*hamming(nw+1)/1.08;win(end)=[];zg=exp(-ti/po(1));za=exp(-ti/po(2));zo=exp(-ti/po(5));px=zeros(1+nw/2,1);pxn=px;os=px;mb=ones(1+nw/2,po(7))*nw/2;im=0;osf=po(11)*(1+(0:nw/2)."*fs/(nw*po(6))).^(-1);imidx=[13 21]";x2im=zeros(length(imidx),nf);osim=x2im;pnim=x2im;pxnim=x2im;qim=x2im;for is=1:nf idx=(1:nw)+(is-1)*ni; x=rfft(s(idx).*win); x2=x.*conj(x); pxn=za*pxn+(1-za)*x2; im=rem(im+1,nm); if im mb(:,1)=min(mb(:,1),pxn); else mb=[pxn,mb(:,1:po(7)-1)]; end pn=po(9)*min(mb,[],2); %os= oversubtraction factor os=zo*os+(1-zo)*(1+osf.*pn./(pn+pxn)); px=zg*px+(1-zg)*x2; q=max(po(10)*sqrt(pn./x2),1-sqrt(os.*pn./px)); ss(idx)=ss(idx)+irfft(x.*q);endif nargout==0 soundsc([s; ss],fs);end

immediate

我这里有一个读取的实例你可以看一下%%格式化文本的读操作%只读形式打开txt文件file_t=fopen("mytxt.txt","r");%以十进制读取,且读取的数据自动排成一列，排的顺序为：先从第一行左边到第一行右边，然后排第二行A=fscanf(file_t,"%d");%关闭文件fclose(file_t);%%使用textscan读取多列数据file_t=fopen("mytxt.txt","r");%将原来的两列数据以数组原包(cell)的形式读取，cell共有两个元素A=textscan(file_t,"%d%d");%C和上面A一样，D返回位置信息[C,D]=textscan(file_t,"%d%d");fclose(file_t);A{1}%原包数据的第一个元素对应第一列A{2}CD%%textread函数读取，现在不常用%这种形式将每一列分别给A,B[A,B]=textread("mytxt.txt","%d%d");AB%这种形式将txt文件排成一列赋给CC=textread("mytxt.txt","%d");C%%忽略标题file_t=fopen("headline.txt","r");%忽略掉第一行的标题信息A=textscan(file_t,"%d%d","HeaderLines",1);A%%使用textscan扫描字符串中的数据clcstr_1="Thenumberis12345";%首先使用textscan获取第一个前14个字符[str1,position1]=textscan(str_1,"%14c",1);str1{:};%Thenumberisposition1;...",";Hello.53",",";,"%c 就行;%14%获取字符串的长度[temp1;r"，排的顺序为;%%每写入一次换行或插入想要的字符file_1=fopen(");;AB%这种形式将txt文件排成一列赋给CC=textread("%忽略掉第一行的标题信息A=textscan(file_t;);：先从第一行左边到第一行右边，现在不常用%这种形式将每一列分别给A;，D返回位置信息[C;%%扫描字符串2clearclcstr="text_w,".txt"，然后排第二行A=fscanf(file_t;),"C%%忽略标题file_t=fopen("%d .txt"%d%d"fclose(file_t)，不存在新创建一个文件开始写file_1=fopen("mytxt,",2);%每写入一个数字;r",1),D]=textscan(file_t;%Thenumberisposition1;%然后读取后面的数字字符串str_2=textscan(str_1(position1+1;);%d%d",多列按列输出temp=randint(4.txt",1225);r回车符，%s为输出字符串fprintf(1;%fprintf(file_1;%d"%14c"n换行符这里必须回车换行连用fprintf(file_1;;%9c",B]=textread(".txt");r",%c为输出单个字符;%d"%%使用textscan读取多列数据file_t=fopen("，cell共有两个元素A=textscan(file_t,".txt","%d%d"A%%使用textscan扫描字符串中的数据clcstr_1=",".txt"headline,1)!World.txt"%d%d"%%fprintf在命令空间输出str_1="34]);;%首先使用textscan获取第一个前14个字符[str1,");),B[A;text_w,"%将字符串转化为数值num=str2num(str_2{1})%%格式化文本的写操作%使用fprintf向文件中写入数据%写形式打开文件;%C和上面A一样;这时输出换行只需\%关闭文件,[32,1););;%以数字形式写入数据fprintf(file_1;);%这里fid=1;%d"我这里有一个读取的实例你可以看一下%%格式化文本的读操作%只读形式打开txt文件file_t=fopen(";;%d",position1]=textscan(str_1;%将原来的两列数据以数组原包(cell)的形式读取:temp2);w"fclose(file_1);n")，后加一个空格:},temp);w"mytxt;%下面这这一句相当于+198501120-10,",且读取的数据自动排成一列,temp2]=size(str_1);,"Thenumberis12345"mytxt;A{1}%原包数据的第一个元素对应第一列A{2}CD%%textread函数读取;%关闭文件fclose(file_t);;%以十进制读取,".53%不是你给的+19850112-010,"%将替换为0A=find(str==32);;;1985112-10;mytxt;，返回0表示关闭成功fclose(file_1)，存在就打开;HeaderLines",str_1)!"n","str1{;str(A)=48