赞
踩
本文适合哪些对象?
最后一句:基础内容其实并不难,不要被外部的一些机构,书籍等误导,请对本文的所有内容完全掌握,内容很少,花两天学会,后续请花时间学习算法。
公众号/CSDN/知乎:川川菜鸟
基本加减乘除:
a=5+5 %加法,同理减法
b=2^3 %立方,
c=5*2 %乘法
x = 1; %分号用于不显示
y = x + 5
常见变量:
在上面小节,我们看到的左边就是字母就叫做变量。很简单,这里随便举几个例子就学会了。
(1)计算b=a*2,a=5
clear all %清空
clc % 情况工作区
a=5;
b=a*2;
b
(2)通过用户输入来获取变量x的值,计算y=x*3
x=input('请输入数字') %input用于接收数字
y=x*2
(3)使用预定变量,常见的比如pi,inf,eps,NaN
r=1
L=2* pi *r
S=pi*r*r
(4)使用常见数学运算,比如开根、求sin(x)
sqrt(78); %百分号用于不显示结果
sin(ans)
这里需要注意,不给表达式赋值时,会自动赋予ans,比如上面写全了就是"ans = sqrt(78)"。
(5)一行多个分配。前面的步骤加封号,只显示最后一个语句的运行结果。
a = 2; b = 7; c = a * b
(6)使用who和whos命令。who查看已使用的所有变量名称;whos查看变量类型、内存、是否为复杂变量等详细信息,:
who
whos
(7)清楚指定变量,或全部变量。使用clear命令。
clear a %清楚a变量
clear %清楚所有变量
(8)可以使用省略号 (…) 将长赋值扩展到另一行。
initial_velocity = 0;
acceleration = 9.8;
time = 20;
final_velocity = initial_velocity...
+ acceleration * time
默认情况下,MATLAB 四个小数位值显示数字。这就是所谓的 short format。如果想更精确,需要使用 format 命令。长(long ) 命令格式显示小数点后16位。
注:一般来说,我们不会太在乎精度,能保留四位数就可以了,所以这个不怎么用得到。
(1)通过以下两种方式对比即可理解,注意逗号分开:
format long ,pi
format short ,pi
再来个例子:
format short %直接声明为短型
a=3.1415926
(2)format bank 命令将数字四舍五入到小数点后两位。
format bank
daily_wage = 177.45;
weekly_wage = daily_wage * 6
(3)格式short e 命令允许以指数形式显示四个小数位加上指数。
format short e
4.678 * 4.9
(4) 格式 long e 命令允许以指数形式显示,小数点后四位加上指数。
format long e
x = pi
(5) format rat 命令给出了计算结果最接近的有理表达式。
format rat
4.678 * 4.9
常见大于等于小于,如果为真返回1,否则返回0:
a=5
b=2
a>b %大于
a==b %等于
a<b %小于
就是数组,矩阵的形式,举两个例子就知道了。
(1)手动输入数组。元素时间用空格/逗号,分号分割。其中逗号与空格等价。分号用于换行,制作多行多列矩阵。
r = [7 8 9 10 11]
r2 = [7 8 9 10 11]
c = [7; 8; 9; 10; 11]
d=[1 2 ;3 4 ;5 6]
两个向量相加等价于对应元素相加,比如
(2)冒号创建。格式:x=开始:间隔:结尾
x=2:2:10
(3)linspace函数创建线性间隔向量。基本格式:linspace(开始,结束,元素个数),默认个数为100。
y=linspace(1,20,10)
矩阵相乘,自然是点积,不是直接的乘法。以几个例子作为说明。
(1)直接.*
a=[1 2 3 4 5];
b=[6 7 8 9 10];
c=a.*b
(2)使用dot函数点积,对应元素相乘再相加,相当于上面的sum(c )
d=dot(a,b)
使用cross函数,基本格式为:cross(a,b)。前提:a和b必须为三维向量。例如:
a=[1 2 3 ];
b=[6 7 8 ];
c=cross(a,b)
其实在前面已经提到过了,这里再强调一次。
x=input("请输入数字:");
x
基本格式:disp(x),显示变量x的值。例如:
a=5*3;
b=a*5;
disp(b)
pause用于暂停程序。例如下面只执行pause上面语句:
a=5*2
pause
b=a*2
语法格式为:
if 表达式
语句组
end
例如:比较a和b的大小
clear all
clc
a=5;
b=2;
if a>b
fprintf("a大于b\n")
end
例如:
clear all
clc
a=5;
b=2;
if a>b
fprintf("a大于b\n")
else
fprintf("a不大于b\n")
end
例如:输入一个数字,判断奇数和偶数
clear all
clc
x=input('请输入数字:');
c=rem(x,2) %求余数的函数
if c==1
fprintf('x为奇数\n')
elseif c==0
fprintf('x为偶数\n')
else
fprintf('非数字')
end
基本格式:
if 表达式
语句组
if 表达式
语句组
end
end
例如:判断学号是否正确,且为偶数
clear all
clc
y=input('请输入学号:')
if y==124
if rem(y,2)==0
fprintf('学号正确,且为偶数\n')
end
else
fprintf('学号错误\n')
end
witch 模块执行多个选项中一组语句。每个选项都包含表达式中。有一情况为真时,MATLAB 执行相应的语句,然后退出开关模块。
switch <开关表达式>
case <条件表达式1>
<语句>
case <条件表达式2>
<语句>
...
...
otherwise
<语句>
end
给出一个评价成绩等级的例子:
grade = 'B';
switch(grade)
case 'A'
fprintf('Excellent!\n' );
case 'B'
fprintf('Well done\n' );
case 'C'
fprintf('Well done\n' );
case 'D'
fprintf('You passed\n' );
case 'F'
fprintf('Better try again\n' );
otherwise
fprintf('Invalid grade\n' );
end
基本格式:
switch(ch1)
case 'A'
fprintf('This A is part of outer switch');
switch(ch2)
case 'A'
fprintf('This A is part of inner switch' );
case 'B'
fprintf('This B is part of inner switch' );
end
case 'B'
fprintf('This B is part of outer switch' );
end
一个嵌套switch语句的例子:
a = 100;
b = 200;
switch(a)
case 100
fprintf('This is part of outer switch %d\n', a );
switch(b)
case 200
fprintf('This is part of inner switch %d\n', a );
end
end
fprintf('Exact value of a is : %d\n', a );
fprintf('Exact value of b is : %d\n', b );
基本格式:
for 变量=表达式
执行语句
end
表达式的通常格式为f:s:l,s默认为1,m表示开始,l表示结束。举个例子说明:遍历数字1到10
for a = 1:10
fprintf('值a: %d \n', a);
end
同理可以设置间隔:
for a = 1:2:10
fprintf('值a: %d \n', a);
end
间隔也可以为负数:
for a = 1.0: -0.1: 0.0
disp(a)
end
右侧也可以是向量:
for a = [2 3 4 5 6]
disp(a)
end
基本格式:
while 表达式
语句组
end
举个简单的例子:遍历小于10的数字,并累加求和
i=0;
sum=0;
while(i<10)
sum=sum+i;
i=i+1;
end
fprintf('累加和为: %d \n', sum);
(1)for嵌套循环基本格式:
for m = 1:j
for n = 1:k
语句;
end
end
一个例子:使用嵌套的 for 循环来显示从 1 到 100 的所有质数。
for i = 2:100
for j = 2:100
if(~mod(i,j))
break; % 可以被除尽,就不是质数
end
end
if(j > (i/j))
fprintf('%d 是质数\n', i);
end
end
编写函数一般在.m文件中编写。函数文件由function语句引导,基本结构为:
function 输出形参表 = 函数名(输入形参表)
% 写一点关于函数的摘要
% 写一点参数说明
函数体语句
end
手动新建函数:
执行函数的形式:
[输出实参表]=函数名[输入实参表]
例如编写一个函数:如果x>0,则y=2x;如果x<0,则y=xx;如果x=0,则y=0
function y= f(x)
% x为一个整数
if x>0
y=2*x;
elseif x==0
y=0;
else
y=x*x;
end
disp(y)
end
执行:
这就是一个常见的普通函数。文件名与函数名保存是一致的。
匿名函数基本格式:
f = @(变量列表)表达式
举两个例子如下(niming.m):
close all;
clc;
b=@(x)(sin(x)+cos(x)); %创建关于x的匿名函数
c=b(pi) %带入pi
d=@(x,y)(sin(x)+cos(y)) %创建x和y的匿名函数
e=d(pi,pi) %带入pi值测试
简单的说,就是函数内部再写了函数,一般写两层。建议没事不要写嵌套函数,函数尽量分开写。
嵌套函数基本格式:
function x = A(p1, p2)
...
function y = B(p3)
...
end
...
end
举个例子:计算一元二次方程(qiantao2.m)
function [x1,x2] = qiantao2(a,b,c)
%定义函数qiantao2计算一元二次函数
d=0; %初始化
function disc
d = sqrt(b^2 - 4*a*c); %计算d
end
disc; %输出d
x1 = (-b + d) / (2*a);
x2 = (-b - d) / (2*a);
end
测试案例:qiantao2(1,-2,1)
再来个测试案例:qiantao2(2,4,-4)
function avg = average(nums)
global TOTAL
avg = sum(nums)/TOTAL;
end
工作区输入以下代码,依据全局变量计算其平均值:
global TOTAL;
TOTAL = 10;
n = [34, 45, 25, 45, 33, 19, 40, 34, 38, 42];
av = average(n)
(1) plot()命令打开基本图窗
导入自变量和因变量,绘制一次函数:
x = [0:5:100];
y = x;
plot(x, y)
(2)xlabel、ylabel、title、grid on命令
x = [0:0.01:10];
y = sin(x);
plot(x, y), xlabel('x'), ylabel('Sin(x)'), title('Sin(x) Graph'),
grid on, axis equal
(3) plot(x,y,x1,y1)
同一个图形绘制多个函数,只需在plot()命令中传入多个x,y即可。
x = [0 : 0.01: 10];
y = sin(x);
g = cos(x);
plot(x, y, x, g, '.-'), legend('Sin(x)', 'Cos(x)')
(4) 图形上设置颜色
下面是常见的颜色代码表:
一个运用不同颜色绘图的示例:
x = [-10 : 0.01: 10];
y = 3*x.^4 + 2 * x.^3 + 7 * x.^2 + 2 * x + 9;
g = 5 * x.^3 + 9 * x + 2;
plot(x, y, 'r', x, g, 'g')
(5)axis ()设置轴刻度
设置轴刻度基本格式:
axis ( [xmin xmax ymin ymax] )
一个手动设置轴刻度的例子
x = [0 : 0.01: 10];
y = exp(-x).* sin(2*x + 3);
plot(x, y), axis([0 10 -1 1])
(6) Sub-Plots生成子图
subplot()基本格式
subplot(x,y,z) %共x行y列现在第z个
一个创建子图的例子:
x = [0:0.01:5];
y = exp(-1.5*x).*sin(10*x);
subplot(1,2,1)
plot(x,y), xlabel('x'),ylabel('exp(–1.5x)*sin(10x)'),axis([0 5 -1 1])
y = exp(-2*x).*sin(10*x);
subplot(1,2,2)
plot(x,y),xlabel('x'),ylabel('exp(–2x)*sin(10x)'),axis([0 5 -1 1])
举个例子,有如下的数据,分别表示为年龄,人体血压高压:
x=[75;78;51;82;77;88;41;78;78;61;71;74;62;81;75;64;80;72;51;80;56;73]
y1=[208;146;168;149;208;102;130;155;163;154;145;147;143;161;145;120;153;158;123;163;177;148]
最简单的方法,比如绘制年龄与高血压关系,使用scatter函数:
scatter(x,y1) %直接绘制
完整为:
x=[75;78;51;82;77;88;41;78;78;61;71;74;62;81;75;64;80;72;51;80;56;73];
y1=[208;146;168;149;208;102;130;155;163;154;145;147;143;161;145;120;153;158;123;163;177;148];
scatter(x,y1)
绘制如下:
想要它好看点呢?绘制为红色,添加个r参数即可,它为red缩写:
scatter(x,y1,'r') %填充颜色绘制为红色
绘制如下:
如果你还想绘制更加好看呢?多彩一点!
c = linspace(1,10,length(x)); %创建向量
scatter(x,y1,[],c)
绘制如下:
再来,填充并标记的散点色彩图:
sz = 25;
c = linspace(1,10,length(x));
scatter(x,y1,sz,c,'filled')
绘制如下:
本节完整代码:
x=[75;78;51;82;77;88;41;78;78;61;71;74;62;81;75;64;80;72;51;80;56;73];
y1=[208;146;168;149;208;102;130;155;163;154;145;147;143;161;145;120;153;158;123;163;177;148];
scatter(x,y1)
scatter(x,y1,'r') %填充颜色绘制为红色
c = linspace(1,10,length(x)); %创建向量
scatter(x,y1,[],c)
sz = 25;
c = linspace(1,10,length(x));
scatter(x,y1,sz,c,'filled')
例如张三的各门课成绩分别为:94,54,65,87,使用这个数据,我们将绘制条形图。
x = [1:4]; %x个数一定要与y个数对应
y = [94,54,65,87];
bar(x,y) %使用bar函数
xlabel('科目')
ylabel('成绩')
title('条形图')
绘制如下:
请参考:数学建模学习(5):最全概率分析之图表绘制详解
本片文章有着详细的教程,因此我便不再重述。
请参考:数学建模学习(46):读取excel数据plot画图
参考matlab官方文档:
如果你还想学习更多关于Matlab的一些函数使用方法,请持续学习文档。
MATLAB提供各种方法求解微积分和积分问题。这里我们讨论计算函数的极限和验证极限的性质。
(1)基本格式:limit 函数将表达式作为参数,并在自变量变为零时找到表达式的极限。
举个例子,计算当x趋于0时,函数 f(x) = (x3 + 5)/(x4 + 7) 的极限
syms x
limit((x^3 + 5)/(x^4 + 7))
(2)计算x趋于0以外的极限,使用带参数的limit命令。limit(f(x),a),这里a是要趋向的值。
limit((x - 3)/(x-1),1)
再举个例子:
limit(x^2 + 5, 3)
使用符号包计算极限值:
x = sym("x"); %让x趋向于0
subs((x^3+5)/(x^4+7),x,0)
首先给出极限的基本性质:
考虑【f(x)=(3x+5)/(x - 3)】、【g(x) =x^2+1】这两个函数,分别计算在5处的极限,以及和、差、积、商四种情况下的极限。给出的代码如下所示:
syms x
f = (3*x + 5)/(x-3);
g = x^2 + 1;
l1 = limit(f, 4)
l2 = limit (g, 4)
lAdd = limit(f + g, 4)
lSub = limit(f - g, 4)
lMult = limit(f*g, 4)
lDiv = limit (f/g, 4)
运行结果较长,这里直接给出答案:17,17,34,0,289,1。
使用符号包时,运用的代码如下:
x = sym("x");
f = (3*x + 5)/(x-3);
g = x^2 + 1;
l1 = subs(f, x, 4)
l2 = subs (g, x, 4)
lAdd = subs (f+g, x, 4)
lSub = subs (f-g, x, 4)
lMult = subs (f*g, x, 4)
lDiv = subs (f/g, x, 4)
结果为17.0,17.0,34.0,0.0,289.0,1.0。
当函数 f(x)的极限在 x = a 处不连续,即x 从左侧接近时的极限值不等于 x 从右侧接近时的极限值。
绘制f(x) = (x - 3)/|x - 3|函数的图片,并求x=3的极值如下所示:
f = (x - 3)/abs(x-3);
ezplot(f,[-1,5])
l = limit(f,x,3,'left')
r = limit(f,x,3,'right')
微分是计算表达式的导数,并在图上寻找局部最大值、最小值,并且求解微分方程。
diff(f)命令计算导数,基本格式:
diff(f) %f为要求导的函数
举例如下:
syms t
f = 3*t^2 + 2*t^(-2);
diff(f)
运行结果如下:
常用求导规则列举如下:
syms x
syms t
f = (x + 2)*(x^2 + 3)
der1 = diff(f)
f = (t^2 + 3)*(sqrt(t) + t^3)
der2 = diff(f)
f = (x^2 - 2*x + 1)*(3*x^3 - 5*x^2 + 2)
der3 = diff(f)
f = (2*x^2 + 3*x)/(x^3 + 1)
der4 = diff(f)
f = (x^2 + 1)^17
der5 = diff(f)
f = (t^3 + 3* t^2 + 5*t -9)^(-6)
der6 = diff(f)
运行结果较长,这里不进行解释,读者只需要知道求解一个函数的导数,①申明求导变量:syms = x ②定义要求导函数:f ③调用diff(f)命令对函数f进行求导
这里给出1个小例子,其他形式一模一样,读者自己尝试就好。
syms x
y = exp(x)
diff(y)
diff(f, n)计算高阶导数,f为求导函数,n为求导的相应阶数。
基本格式:
syms x
f = ...;
diff(f, n)
对y = f(x) = x*e^(-3x)求二阶导如下:
syms x;
f = x*exp(-3*x);
diff(f, 2)
已知原函数,对y = f(x) = 3 sin(x) + 7 cos(5x).求导,满足f" + f = -5cos(2x)。思路是首先计算f",之后计算f" + f,接着计算-5cos(2x),最后检查两遍是否相等即可。代码实现如下:
syms x
y = 3*sin(x)+7*cos(5*x); % defining the function
lhs = diff(y,2)+y; %evaluting the lhs of the equation
rhs = -5*cos(2*x); %rhs of the equation
if(isequal(lhs,rhs))
disp('Yes, the equation holds true');
else
disp('No, the equation does not hold true');
end
disp('Value of LHS is: '), disp(lhs);
若要寻找 f(x) = 2x^ 3 + 3x^ 2 − 12x + 17函数的最大值与最小值,我们可以①用ezplot(y,[ ])绘制出函数最大值与最小值,y是函数,[ ]是自变量取值区间。②调用diff(y)命令求导。③solve(g)命令求解导函数,即令导函数为0,即可解出对应的自变量x。
syms x
y = 2*x^3 + 3*x^2 - 12*x + 17; %定义函数
ezplot(y, [-2, 2])
g = diff(y);
s = solve(g)
之后调用subs()命令评估函数f在1,和-2处的值。
subs(y, 1), subs(y,-2)
基本格式:
dsolve('eqn','cond1', 'cond2',…)
eqn中导数用 D 表示,计算y" - y = 0, y(0) = -1, y’(0) = 2,我们可以输入如下代码:
dsolve('D2y - y = 0','y(0) = -1','Dy(0) = 2')
MATLAB 将多项式表示为包含按降幂排序的系数的行向量。例如,方程 P(x) = x4 + 7x3 - 5x + 9 可以表示为 p = [1 7 0 -5 9];
polyval()函数计算指定值处的多项式,基本格式:
p = [1 7 0 -5 9];
polyval(p,4)
输入自变量为一个矩阵式,可调用polyvalm(y,X),其中y是函数,X是自变量矩阵。
p = [1 7 0 -5 9];
X = [1 2 -3 4; 2 -5 6 3; 3 1 0 2; 5 -7 3 8];
polyvalm(p, X)
root(f)求多项式的根,其中f为要求解的函数,基本格式:
p = [1 7 0 -5 9];
r = roots(p)
ploy®返回多项式系数,其中r是我们已知的根,比如在上例中,r就是我们已知根,要想求对应系数,计算对应系数的代码如下:
p2 = poly(r)
已知x和y为两个向量,是要拟合到 n 次多项式的 x、y 数据,可使用polyfit()函数,基本格式如下:
x = [1 2 3 4 5 6]; y = [5.5 43.1 128 290.7 498.4 978.67]; %data
p = polyfit(x,y,4) %get the polynomial
借助前面绘图的知识,我们可将拟合出的曲线打印在图上,代码如下:
x = [1 2 3 4 5 6]; y = [5.5 43.1 128 290.7 498.4 978.67]; %data
p = polyfit(x,y,4) %get the polynomial
x2 = 1:.1:6;
y2 = polyval(p,x2);
plot(x,y,'o',x2,y2)
grid on
请参考本专栏:数学建模丛小白到精通之路,本专栏一直更新,从未断过。
本专栏使用的编程语言:
我对你的期望:
我的祝愿:希望大家拿奖到手软!——川川菜鸟
关注公众号/CSDN/知乎:川川菜鸟
获取更多数模资料
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。