常用基础命令
- 通过linspace(start, end, num)可以生成一个包含start和end总共num个元素的等差向量。
>>result = linspace(1, 5, 5)
result = 1 2 3 4 5
>>result = linspace(0, 10, 10)
result = 0 1.1111 2.2222 3.3333 4.4444 5.5556 6.6667 7.7778 8.8889 10.0000
- nargin、nargout可以在函数体内表示输入参数、接收参数的个数
function [a,b] = nargTest(c,d)
if nargin == 1 %可以根据调用此函数时输入参数的个数进行不同的操作
a=c;
b=c;
else
a = c+d;
b = c-d;
if nargout==2 %可以根据调用此函数时接收参数的个数进行不同的操作
a = 0;
end
end
- 一些简单函数的求值不必编写function函数,采用inline即可实现
f = inline('x^2+2*x+1', 'x'); %后面的参数x也可不写
f(0) %求解函数在0处的函数值
f(1) %求解函数在1处的函数值
plot(x,f(x)) %画出函数图像
inline函数在新版本中将要被取缔了,建议使用匿名函数定义方式
@(参数列表)(函数表达式)
f = @(x)(x^2+2*x+1);
以上两种方式中,如果想要同时计算多个自变量对应的函数值(输入参数为向量),则函数表达式中的* \ / ^等符号需要加点
。
- 可以通过dot(A, B)来计算两向量的点积
- rand(m, n) 可以被用于产生一个在[0, 1]之间均匀分布的m行n列随机矩阵;而 randn(m, n) 则用于产生一个均值为0方差为1的m行n列高斯随机矩阵(
白噪声
) - zores(m, n)、ones(m, n) 根据上面可见名知意
一元函数绘图
- 直角坐标下作函数图像的方法有plot、ezplot、fplot等,同时可通过legend()命令对图像进行注释说明(图例),运行如下代码:
%% 通过plot来绘图
x = 0:0.1:4*pi;
y1 = exp(-0.5*x);
y2 = y1.*sin(5*x); %由于x是向量,所以y1和后面的部分均为向量,需要用 .*
plot(x,y2, x,y1,'--r', x,-y1,'--r') %画图并设置线型
legend('exp(-0.5x)sin(5x)', 'exp(-0.5x)') %设置图例
%% 通过ezplot来绘图
ezplot('exp(-0.5*x)*sin(5*x)', [0, 14, -1, 1]) %画图并设置显示区间
此外还可用title()为图形设置图名,用xlabel()、ylabel()设置坐标轴。
- 在极坐标下作函数图象采用polar,运行如下代码:
%% 极坐标下绘图
theta = 0:pi/20:4*pi;
r = 1+2*theta;
subplot(1,2,1);
polar(theta, r, 'r');
theta = 0:0.001:2*pi;
r = sin(5*theta);
subplot(1,2,2);
polar(theta, r, 'b');
- 通过bar、bar3命令画条形统计图,运行如下代码:
%% 条形图绘制
x = -2.9:0.2:2.9;
y = exp(-x.*x);
subplot(1,2,1);
bar(x,y);
title('Figure:2D Bar Chart');
subplot(1,2,2);
bar3(x,y);
title('Figure:3D Bar Chart');
- 通过xlim、ylim等函数来指定坐标轴范围,代码如下:
x = linspace(-10,10,200);
y = sin(4*x)./exp(x);
plot(x,y)
xlim([0 10])
ylim([-0.4 0.8])
二元函数绘图
- 通过plot3命令进行空间曲线的绘制,运行如下代码:
%% 空间曲线的绘制
t = linspace(0,10*pi);
x = sin(t);
y = cos(t);
z = t;
plot3(x, y, z);
- 空间曲面的绘制一般分为三步:生成平面网格、计算网格点上函数值、绘制网面
执行如下代码:
%% 空间曲面的绘制
[x,y] = meshgrid(-2 : 0.1 : 2);
z = x.*exp(-x.^2-y.^2);
subplot(1,2,1);
mesh(x,y,z);
%通过meshc绘制具有等势线的曲面图
subplot(1,2,2);
meshc(x,y,z);
- 通过contour、contour3命令来绘制等高线图,运行如下代码:
%% 等值线图绘制
[x,y] = meshgrid(-2 : 0.1 : 2);
z = x.*exp(-x.^2-y.^2);
subplot(1,2,1);
contour(x,y,z); %后面还可以加一个数字参数,指定画多少根等值线
title('Figure:2D Counter Plot');
%通过meshc绘制具有等势线的曲面图
subplot(1,2,2);
contour3(x,y,z); %后面还可以加一个数字参数,指定画多少根等值线
title('Figure:3D Counter Plot');