欢迎光临散文网 会员登陆 & 注册

数学建模比赛MATLAB从入门到精通教程

2021-12-12 21:55 作者:公孙灼月  | 我要投稿



Empty state-space model.

是因为你的文件命名和matlab库里的重名了。改一下你的名字就好了



2.1.1 MATLAB系统环境 P2 - 10:44





结构与单元矩阵 一个变量套变量,一个直接变量套数据
















远赴人间惊鸿宴,一睹人间盛世颜









左乘乘在行,右乘乘在列









特征值:eigenvalue






脚本是程序









黄金分割比


长宽比=宽比长减宽













&&是与符号,相当于串联电路,要两个灯泡都没坏才能接通

&&就相当于是数学中的且,&和&&稍有区别在这里都可以,但是编程中一般都是用&&

因为题目说了是一个字符,所以不能是两位数,自然只能是0到9了












这个用处可大了 到了主成分分析 和降噪的时候会用到协方差矩阵和其特征值以及其特征向量,推荐去看看3blue1brown的线性代数的本质




fill是填充颜色的

fill是绘图然后填充

subplot指定了绘图子区域







sp转稀疏






自动排列了而已

类似字典把

就是存储嘛,前面是坐标



带状稀疏矩阵不就是对角阵嘛



B的每一列是非零的对角线,d的每一行非零对角线的位置序号

k=-3,k=0,k=3时的三条对角线对应d

前面有一节讲过,对角线是0,往下是负的,往上是正的


真变态对角线是0从上面数是正下面是负。



spdiags提取并创建稀疏带状和对角矩阵


kf1和前面(3)的第一列写法不一样


工程上阶数比较高的三对角矩阵肯定用这个方法最好


swich选择


单元数据 大括号


switch只执行一个语句,第二个语句不论怎么样都不会执行



s的意思是把输入的东西当成字符串存起来

‘s’是字符串形式

's'表示 返回输入的文本,而不会将输入作为表达式来计算。






C = num2cell(A) % 将数值阵列 A 转成异质阵列 C

num2cell 将数值矩阵转为单元矩阵

第一步产生的是数值向量,先将数值向量转为单元数据

{1,2,3}写入MATLAB为cell类型不转化的话是double类型














是arctanx的泰勒展开,当x=1的情况










x,y 所投的坐标




第一个 行向量 循环 4次

第二个列向量 一次

因为第二句是一个列向量,而第一句可以看作是四个列向量

for语句是执行所有元素,当变量k赋值为矩阵的时候,列就是元素,列数就是元素个数,当k执行到最后一个元素,有值就赋到k



msum存放最后的累加和


n 输入值的个数

~= 不等于





continue此时挨if最近,所以只是跳出if而不是for







输入输出参数都没用分配存储空间 形式参数



建议二者统一




形式参数与实际参数


顺序个数一致


多个输入参数用,分割






matlab中log指的就是数学上的ln










前8项平方和等于第8*9项


sum函数只有是矩阵的时候才是对每一列的元素求和













横坐标是x的下标,纵坐标是x的值





列等行为纵坐标,行数为条数

经鄙人验证,方阵时以每一列为纵坐标


有一个x,就有一个sin

这里的y其实由三行组成,每行100个,其列数与x向量长度相等,所以以x向量为横坐标,y的每行为纵坐标画图

linspace(a,b,c) a,b是区间首尾 c是中间取的拟合点的个数






















可以使用代码直接在图形窗口里面的高级选项修改


大括号








第二个指出特定点






先横后纵







grid 用于切换

默认无网格




默认有边框




保持

刷新











背景板纵向分成四份,横向分成四份,选第八个区

利用坐标轴对象操作对图形窗口做任意分割

subpolt中前两个数表示把窗口分成几乘几份,然后第三个数表示选中其中的第几个



第一个 横坐标为logx,纵坐标不变



旋转图形 极角+旋转角度 角度为正 顺时针















hist 直角坐标系

rose 极坐标系








将y转为弧度制




x通常为向量

省略explode时图为一个整体


secore,ex 均为向量

逆时针

ex中1指分离

location 指定图例位置

eastoutside 图例于右边外侧

默认于图内,可能会遮挡图




如果觉得这些画图比较难的,matlab左上角绘图里有现成在工作区

选择要绘图的变量后就可以在绘图中选择需要的样式,和excl相似的模板


filled数据点填充标记 实心

省略为空心




大小方向轨迹

省略想想x,y 则在xy平面均匀取若干点






axis 确定x,y,z的范围,若不写则为0




1 xyz条数=矩阵列数

2 向量长度与矩阵相同 行向量的长度=矩阵列数 列向量=矩阵行数




t=t‘ 将行向量装置为列向量’
















前面的,这个三行二列是从下到上、从左到右数的第几条线


那个y是列向量,'是转置的意思!





mesh函数画网格图,surf函数画曲面图

结构体内容引用自非结构体数组对象




 






s 球面 圆心于原点,半径为一

n默认20 员的光滑程度

c r为标量时


圆柱面





改为-1这样柱面高度为原来两倍,方便接下来相交

因为默认值都是1,改了就是为了方便看图的,并且只有改了才相切

z的改变将柱形高度改为和底面圆的直径一样,使得相交图形的效果好











方位角 正值 逆时针


默认方位角-37.5

仰角30





2 2维 方位角 0 仰角 90

3 3维

4 默认



















最后一个最光滑


cylinder - 生成圆柱

圆柱半径从pi到0变化,密度是10











绘图时先先选中为x的值,y z ……











A(:)即为将A的所有元素排成列向量





















dim是按列或按行计算








列与列之间的相关性,只有4*4种结果

因为要找最接近销量的方案,第一列是销量

















g1 是为了保证f,g的幂次相等

相乘后x最高变成了7次,所以项数变成了7+1=8项

用g1也是同样的效果,只不过结果前面会多出几个零,而实际对方程没有任何影响,因此也可用g1

加减运算是向量运算,要补足;乘除是调用函数,软件已经帮你写好了过程




diff也可以求导









先求商可能产生余式





一次并行求多个多项式的值,自变量有多少就求几个,各个元素互不干扰





ones产生全一矩阵,eye产生对角线为1的矩阵



输入单位矩阵就相等












6.3.1 数据插值 P31 - 01:23


6.3.1 数据插值 P31 - 01:23


6.3.1 数据插值 P31 - 01:23










感觉加一个 axis equal 更能表示机翼线条的形状










对于(2)来说 如果……的值










这个计算量应该是说计算的复杂程度,因为只计算相邻的两个数值。










6.4.1 数据插值应用举例 P32 - 00:27


6.4.1 数据插值应用举例 P32 - 00:27


6.4.1 数据插值应用举例 P32 - 00:27












倒数第4,绝对值

倒数第3,升序排列




































6.6.1 曲线拟合应用举例 P34 - 01:22


6.6.1 曲线拟合应用举例 P34 - 01:22


6.6.1 曲线拟合应用举例 P34 - 01:22
















7.1.1 数值微分与数值积分 P35 - 03:29


7.1.1 数值微分与数值积分 P35 - 03:29


7.1.1 数值微分与数值积分 P35 - 03:29










x =


    1    2    3    4

    2    3    4    5

    3    4    5    6


>> dx=diff(x,1,1)


dx =


    1    1    1    1

    1    1    1    1












1e-8=1E-8=10^-8

















1:end-1的个数=2:end





由里到外
















7.4.1 非线性方程求解与函数极值计算 P37 - 02:13


7.4.1 非线性方程求解与函数极值计算 P37 - 02:13


7.4.1 非线性方程求解与函数极值计算 P37 - 02:13






fplot(x,x-1./x+5)

>> hold

已锁定最新绘图

>> grid

>>







x(1)=x; x(2)=y x(3)=z










x0是x1和x2的迭代初值向量

A=,b=, 定义线性不等约束, lb 解的下限

x0是初值,A b线性不等式约束,Aeq beq线性等式约束,lb下边界,ub上边界,nonlcon非线性约束条件


每两个逗号对应前面文字部分的约束类别 且用系数矩阵和向量表示









7.5.1 常微分方程数值求解 P38 - 01:16


7.5.1 常微分方程数值求解 P38 - 01:16


7.5.1 常微分方程数值求解 P38 - 01:16
















ode后面的(函数,x轴上取值范围,初值条件),大致可以总结成这样

第一行是矩阵相乘变成了-2*x(2);x(1)对应前一页PPT公式

左边那个是在时间域下系统频率随时间变化的曲线,右边那个是在幅频域下的相频特性曲线

这里x=[x1,x2]=[x',x]







tic记录语句开始时间

toc记录语句结束时间



‘计算的点数’num……

ode15s 's'专门求解刚性问题





















































































8.1.1 符号对象 P40 - 01:12

























8.2.1 符号微积分 P41 - 01:30
















8.3.1 级数 P42 - 00:08
















8.4.1 符号方程求解 P43 - 01:54




































9.2.1 曲线与曲面对象 P45 - 02:06





圆的参数方程




color 默认白色

gouraud 先对顶点颜色插补,再对顶点勾画的面上的颜色插补,用于表现曲面

phong 对顶点处的法线插值,再计算各个像素的反光,光照效果好



axis 设置坐标轴范围和纵横比

axis( )设置方位角

hs1 曲面无网格线,差值填充

light 放置相同光源 位置为【0,0,8】

lighting 设置光照模式

specularstrength 光泽


9.3.1 图形用户界面设计方法 P46 - 01:34







回调函数 定义对象怎样处理信息并响应某事件

source源控件对象句柄

ev~存储事件数据






上为对mline ——type的定义




先运行例子1

后运行例子2


9.4.1 用户界面设计工具 P47 - 01:15


guide

图形用户界面开发环境




black 空模板

UIcontrols 带控件对象的图形用户界面模板

axes and menu 带坐标轴与菜单的图形用户界面~

modal 带模式对话框的图形~
















注意大小写,Value





9.5.1 APP设计工具 P48 - 13:08


我用的是2018a版本,上面那个同学说的对,是把APP.ButtonGroup下面第一个改成RadioButton就可以了。









数学建模比赛MATLAB从入门到精通教程的评论 (共 条)

分享到微博请遵守国家法律