matlab-bb

matlab字符串的⽐较,MATLAB字符串⽐较总结% 字符串处理a=' a';b='b b';c='cccc';m=''% 获取字符串长度length(a)% 连接两个字符串,每个字符串最右边的空格被裁切d=strcat(a,c)length(d)% 连接多⾏字符串,每⾏长度可不等，⾃动把⾮最长字符串最右边补空格% 使与最长字符串相等，会忽略空字符串e=strvcat(a,b,m)size(e)% char连接，空字符串会被空格填满f=char(a,b,m)size(f)% strcmp ⽐较两个字符串是否完全相等，是，返回真，否则，返回假% strncmp ⽐较两个字符串前n个字符是否相等，是，返回真，否则，返回假% strcmpi ⽐较两个字符串是否完全相等，忽略字母⼤⼩写% strncmpi ⽐较两个字符串前n个字符是否相等，忽略字母⼤⼩写% isletter 检测字符串中每个字符时否属于英⽂字母% isspace 检测字符串中每个字符是否属于格式字符(空格，回车，制表，换⾏符等)% isstrprop 检测字符每⼀个字符是否属于指定的范围a='d sdsdsd 15#';b=isletter(a)c=isspace(a)% 字符串替换和查找% strrep 进⾏字符串替换，区分⼤⼩写% strrep(str1,str2,str3)% 它把str1中所有的str2字串⽤str3来替换% strfind(str,patten) 查找str中是否有pattern，返回出现位置，没有出现返回空数组% findstr(str1,str2) 查找str1和str2中，较短字符串在较长字符串中出现的位置，没有出现返回空数组% strmatch(patten,str) 检查patten是否和str最左侧部分⼀致% strtok(str,char) 返回str中由char指定的字符串前的部分和之后的部分，mm='youqwelcome';[mm1,mm2]=strtok(mm,'q')% blanks(n) 创建有n个空格组成的字符串% deblank(str) 裁切字符串的尾部空格% strtrim(str) 裁切字符串的开头和尾部的空格，制表，回车符% lower(str) 将字符串中的字母转换成⼩写% upper(str) 将字符串中的字母转换成⼤写% sort(str) 按照字符的ASCII值对字符串排序% num2str 将数字转换为数字字符串% str2num 将数字字符串转换为数字% mat2str 将数组转换成字符串% int2str 把数值数组转换为整数数字组成的字符数组------------------------------CELL数组中的字符串⽐较：c=cell(2,1);c(1,1)=cellstr('xxx');c(2,1)=cellstr('yyyyyyy');strcmp(c{1,1},c{2,1});------------------------------isequal Test arrays for equality, 可⽤来⽐较两个字符数组是否相同。

MATLAB操作命令大全1.基本操作：- clear: 清除工作区中的所有变量。

- clc: 清除命令窗口的内容。

- close all: 关闭所有图形窗口。

- help function-name: 显示与函数相关的帮助文档。

- who: 显示当前工作区中的所有变量。

- save file-name: 保存当前工作区中的所有变量到指定的文件。

- load file-name: 从文件中加载变量到当前工作区。

2.变量操作：-=:赋值操作符，将右边的值赋给左边的变量。

-+:加法操作符。

--:减法操作符。

-*:乘法操作符。

-/:除法操作符。

-^:幂运算操作符。

- sqrt(x): 计算 x 的平方根。

- abs(x): 计算 x 的绝对值。

- max(x): 返回 x 中的最大值。

- min(x): 返回 x 中的最小值。

- sum(x): 计算 x 中所有元素的和。

3.数组操作：- zeros(m, n): 创建一个 m 行 n 列的全零数组。

- ones(m, n): 创建一个 m 行 n 列的全一数组。

- eye(n): 创建一个 n 行 n 列的单位矩阵。

- size(x): 返回 x 的维度。

- length(x): 返回 x 的长度。

- reshape(x, m, n): 将 x 重新排列为一个 m 行 n 列矩阵。

- transpose(x): 将 x 的行和列互换。

4.控制流程：- if-else: 条件语句，根据条件执行不同的代码块。

- for loop: 循环语句，执行指定次数的代码块。

- while loop: 循环语句，根据条件反复执行代码块。

- break: 在循环中使用，用来跳出当前循环。

- continue: 在循环中使用，用来跳过当前循环的剩余部分。

5.统计分析：- mean(x): 计算 x 的平均值。

- median(x): 计算 x 的中位数。

- std(x): 计算 x 的标准差。

在MATLAB 中，你可以执行各种数学运算来操作变量。

以下是一些常见的例子：1. 加法：你可以使用"+" 符号来执行加法运算。

例如，如果你有两个变量a 和b，你可以通过"a + b" 来计算它们的和。

2. 减法：你可以使用"-" 符号来执行减法运算。

例如，"a - b" 将返回a 和b 的差。

3. 乘法：你可以使用"*" 符号来执行乘法运算。

例如，"a * b" 将返回a 和b 的乘积。

4. 除法：你可以使用"/" 符号来执行除法运算。

例如，"a / b" 将返回a 除以b 的结果。

5. 指数运算：你可以使用"^" 符号来执行指数运算。

例如，"a ^ 2" 将返回a 的平方。

6. 开方运算：你可以使用"sqrt" 函数来执行开方运算。

例如，"sqrt(a)" 将返回a 的平方根。

7. 对数运算：你可以使用"log" 或"ln" 函数来执行对数运算。

例如，"log(a)" 将返回以e 为底的对数，而"ln(a)" 将返回以自然对数e 为底的对数。

8. 正弦、余弦和正切运算：你可以使用"sin"、"cos" 和"tan" 函数来执行这些运算。

例如，"sin(a)" 将返回 a 的正弦值，其中 a 是以弧度计的角度。

9. 其他三角函数：你可以使用"asin"、"acos" 和"atan" 函数来计算反正弦、反余弦和反正切值。

10. 绝对值：你可以使用"|" 符号来获取一个数的绝对值，或者使用"abs" 函数。

关于MATLAB整数规划分支定界法一、编程利用Matlab的线性规划指令:[x,fval]=linprog(f,A,b,Aeq,beq,lb,ub)编写计算整数规划函数，输入与输出与上述指令相同分枝定界法(递归实现)function [x,fval,status] = intprog(f,A,B,I,Aeq,Beq,lb,ub,e)%整数规划求解函数 intprog()% 其中 f为目标函数向量% A和B为不等式约束 Aeq与Beq为等式约束 % I为整数约束% lb与ub分别为变量下界与上界% x为最优解，fval为最优值%例子:% maximize 20 x1 + 10 x2 % S.T.% 5 x1 + 4 x2 <=24 % 2 x1 + 5 x2 <=13 % x1, x2 >=0 % x1, x2是整数 % f=[-20, -10];% A=[ 5 4; 2 5];% B=[24; 13];% lb=[0 0];% ub=[inf inf];% I=[1,2];% e=0.000001;% [x v s]= IP(f,A,B,I,[],[],lb,ub,,e) % x = 4 1 v = -90.0000 s = 1% 控制输入参数if nargin < 9, e = 0.00001;if nargin < 8, ub = [];if nargin < 7, lb = [];if nargin < 6, Beq = [];if nargin < 5, Aeq = [];if nargin < 4, I = [1:length(f)];end, end, end, end, end, end%求解整数规划对应的线性规划,判断是否有解 options = optimset('display','off'); [x0,fval0,exitflag] = linprog(f,A,B,Aeq,Beq,lb,ub,[],options);if exitflag < 0disp('没有合适整数解');x = x0;fval = fval0;status = exitflag;return;else%采用分支定界法求解bound = inf;[x,fval,status] =branchbound(f,A,B,I,x0,fval0,bound,Aeq,Beq,lb,ub,e);end子函数function [newx,newfval,status,newbound] =branchbound(f,A,B,I,x,fval,bound,Aeq,Beq,lb,ub,e)% 分支定界法求解整数规划% f,A,B,Aeq,Beq,lb,ub与线性规划相同 % I为整数限制变量的向量% x为初始解，fval为初始值options = optimset('display','off');[x0,fval0,status0]=linprog(f,A,B,Aeq,Beq,lb,ub,[],options);%递归中的最终退出条件%无解或者解比现有上界大则返回原解 if status0 <= 0 || fval0 >= bound newx = x;newfval = fval;newbound = bound;status = status0;return;end%是否为整数解,如果是整数解则返回 intindex = find(abs(x0(I) -round(x0(I))) > e);if isempty(intindex)newx(I) = round(x0(I));newfval = fval0;newbound = fval0;status = 1;return;end%当有非整可行解时，则进行分支求解 %此时必定会有整数解或空解%找到第一个不满足整数要求的变量n = I(intindex(1));addA = zeros(1,length(f));addA(n) = 1;%构造第一个分支 x<=floor(x(n))A = [A;addA];B = [B,floor(x(n))];[x1,fval1,status1,bound1] =branchbound(f,A,B,I,x0,fval0,bound,Aeq,Beq,lb,ub,e);A(end,:) = [];B(:,end) = [];%解得第一个分支，若为更优解则替换，若不是则保持原状status = status1;if status1 > 0 && bound1 < boundnewx = x1;newfval = fval1;bound = fval1;newbound = bound1;elsenewx = x0;newfval = fval0;newbound = bound;end%构造第二分支A = [A;-addA];B = [B,-ceil(x(n))];[x2,fval2,status2,bound2] =branchbound(f,A,B,I,x0,fval0,bound,Aeq,Beq,lb,ub,e);A(end,:) = [];B(:,end) = [];%解得第二分支，并与第一分支做比较，如果更优则替换 if status2 > 0 && bound2 < boundstatus = status2;newx = x2;newfval = fval2;newbound = bound2;end二、求解利用上述指令求解下列问题:汽车厂生产三种类型的汽车，已知各类型每辆车对钢材、劳动时间的需求，利润及工厂每月的现有量小型中型大型现有量钢材(吨) 1 2 5 1000劳动时间(小时) 250 125 150 120000利润(万元) 3 5 121、若每月生产的汽车必须为整车，试制订月生产计划，使工厂的利润最大2、如果生产某一类型汽车，则至少要生产50辆，那么最优的生产计划应作何改变, 解:1f = [-3 -5 -12];A = [1 2 5;250 125 150];B = [1000 120000];I = [1:length(f)];lb = [0 0 0];[x,fval,status] = intprog(f,A,B,I,[],[],lb)答案为 x =307 344 1 fval = -2653 status =12lb = [50 50 50]答案为 x =350 200 50 fval =-2.6500e+003 status =1用分枝定界法求解整数线性规划问题%%问题的标准形式:%% min c'*x%% s.t. A*x<=b%% Aeq*x=beq%% 要求是整数%%function [y,fval]=BranchBound(c,A,b,Aeq,beq)%NL=length(c);UB=inf;LB=-inf;FN=[0];AA(1)={A};BB(1)={b};k=0;flag=0;while flag==0;[x,fval,exitFlag]=linprog(c,A,b,Aeq,beq); if (exitFlag == -2) | (fval >= UB)FN(1)=[];if isempty(FN)==1flag=1;elsek=FN(1);A=AA{k};b=BB{k};endelsefor i=1:NLif abs(x(i)-round(x(i)))>1e-7kk=FN(end);FN=[FN,kk+1,kk+2];temp_A=zeros(1,NL);temp_A(i)=1;temp_A1=[A;temp_A];AA(kk+1)={temp_A1};b1=[b;fix(x(i))];BB(kk+1)={b1};temp_A2=[A;-temp_A];AA(kk+2)={temp_A2};b2=[b;-(fix(x(i))+1)];BB(kk+2)={b2};FN(1)=[];k=FN(1);A=AA{k};b=BB{k};break;endendif (i==NL) & (abs(x(i)-round(x(i)))<=1e-7) UB=fval;y=x;FN(1)=[];if isempty(FN)==1flag=1;elsek=FN(1);A=AA{k};b=BB{k};endendendendy=round(y);fval=c*y;。