机械优化设计matlab编程指导201_...

（一）matlab基本用法介绍
5. 循环、结束命令
¾break是直接跳出该层循环
¾continue是直接进入该层循环的下一次迭代¾return是直接退出程序或函数返回
大概的等级关系如下
¾return>break>continue
（二）matlab语法应用实例
1. 循环的用法（用于优化迭代）
%%%%%%%%%%%%%%%%%% while 循环
i=0;
while(i<10)
i=i+1;
end
%%%%%%%%%%%%%%%%%% for 循环
i=0;
for j=1:10
i=i+1;
end
%%%%%%%%%%%%%%%%%% break终止本层循环
for i=1:2
i
k=0;
for j=1:10
k=k+1
if k==5;break;end
end
end
%%%%%%%%%%%%%%%%%% continue直接进入下一次循环变量的迭代
k=0;
for j=1:10
k=k+1;
if k==5;continue;end
k
end %注意，k=5时，后面的循环体未被执行，即k=5不显示在输出窗口中
2 函数调用（利用外部m文件定义函数）（用于定义目标函数、约束函数、优化子程序等）
建立wad2.m文件：
function F1=wad2(x)
F1=x(1)^2+x(2)^2;
在主程序m文件中或者matlab的commond window中输入
x0=[3 5]';
fx0= wad2(x0)
运行结果如下：fx0 = 34 ，得到了函数的值
※另外函数中的F1也可以定义成向量。

3 全局变量（可用于优化过程中的重要参数定义，在程序任意部位都可以访问）
建立wad3.m文件：
function F1=wad3(x)
global ac bc
F1=ac*x(1)^2+bc*x(2)^2;
在主程序m文件中定义
x0=[3 5]';
ac=0.2; bc=0.8;
fx0= wad3(x0)
结果：fx0 = 14.6000
4 判断语句（例如可以用于优化计算收敛的终止条件判断）
k=2
if k>3
'yes k>3'
else
'no ,k<3'
end
if k==2
'yes k=2'
end
5 数组（例如可以用于优化计算收敛的终止条件判断）
% 1. 数组：matlab 程序以矩阵为运算单位，任何变量都是矩阵（数组）% 数组不需要声明，直接定义、使用即可
a=[1 2 3]
b=[ 1 2 3
4 5 6]
c=b'; %求矩阵的转置矩阵
b23=b(2,3) %取出矩阵b中第2行第3列元素
x=[1
2
3]
Hk=eye(4);% 生成4维的单位阵
6 清屏
clc
clear
7 二维绘图
x=1:0.1:10;
y1=sin(x);
y2=cos(x);
figure;
plot(x,y1,'r');
hold on;
plot(x,y2,'b');
grid on;
xlabel('x axis'),ylabel('y axis');
8 三维绘图
x=1:0.1:10;
y1=sin(x);
y2=cos(x);
z=x+cos(y1);
figure;
plot3(x,y1,z,'b');
hold on;
grid on;
xlabel('x axis'),ylabel('y axis'),zlabel('z axis');
9 二维等值线图
%%对函数： F=(xk(1)-2)^4+(xk(1)-2*xk(2))^2; [x3,y3]=meshgrid(-5:0.2:10,-10:0.2:10);
z3=(x3-2).^4+(x3-2*y3).^2;
contour(x3,y3,z3,40);%目标函数的等值线
hold on;
axis equal;
xlabel('自变量x1','FontSize',12);
ylabel('自变量x2','FontSize',12);
title('×××××××方法', 'FontSize',15);。