MATLAB常用语句总结7

MATLAB总结7:常见错误归纳

本篇专门用于记录一些应试技巧


文章目录

  • MATLAB总结7:常见错误归纳
  • 前言
  • 一、一些小定义和小技巧
  • 二、蒙塔卡罗求解方法
    • 1.函数的定义
    • 2.函数引用
    • 3.代码量较少的蒙塔卡罗
  • 三、函数引用与多变量
  • 四、矩阵引用
  • 五、非线性函数:fmincon的使用
  • 六、线性规划函数
    • 1.linprog
    • 2.fminbnd、fminsearch
  • 七、积分和微分
    • 1.quad积分
    • 2.积分int 微分diff
  • 八、lsqcurvefit
  • 九、矩阵求解
  • 十、polyfit 函数进行多项式拟合
  • 十一、ode23和ode45
  • 十二、fzero
  • 十三、meshgrid和surf
  • 十四、fsolve


前言

记忆力比较差,把学习过程中常用的语法记录于此,方便回忆时索引


一、一些小定义和小技巧

  1. rand的使用

rand 是用于生成在(0,1)上随机数的函数,而后面加入 rand(1,2),则意为一个一行两列的(0,1)矩阵,不可乱加。需要使用不同范围的随机数则使用 A*rand() 即可

    x=5*rand();y=5*rand();
  1. 符号
取余:mod(A,b)%A是要进行取余的数字,b是除数
不等于:~=
  1. find()的使用
    用r=find(A<=0.6),直接找到A中小于0.6的下标

  2. 分段函数的使用

y=(x<=8).* (2*x)+(x>8&&x<12).*(3*x)+(x>12).* (4*x);
  1. upper, lower的大小写转换
for i=1:length(v)if v(i)=='a'r(i)=upper(r(i));elseif v(i)=='h'r(i)='2';elseif v(i)=='v'r(i)='3';end
end

二、蒙塔卡罗求解方法

在这里插入图片描述

1.函数的定义

代码如下(示例):

function c=myobj(x)
c=10*x(1)^3+x(3)*(x(1)^2+x(2)^2);
function r=mycon(x)
r=[x(1)^2+x(2)^2-x(3)-10 ; x(1)^2+x(2)^2+x(3)-3];
function m=myrand
while 1m=[10*rand-5 ;10*rand-5 ;10*rand-5 ];if mycon(m)<=0return end
end

2.函数引用

代码如下(示例):

function [val,x]=myfun2
rand('seed',11);betterx=myrand;
betterval=myobj(betterx);for i=1:10000betterx=myrand;if myobj(betterx)>betterval;betterval=myobj(betterx);end
end
x=double(betterx);
val=double(betterval);

3.代码量较少的蒙塔卡罗

function [val,x]=myfun2
rand('seed',11);
val=0;
for i=1:100000s=-10*rand(1,3)+5;if s(1)^2+s(2)^2-s(3)-10<=0 && s(1)^2+s(2)^2+s(3)-3<=0if 10*s(1)^3+s(3)*(s(1)^2+s(2)^2)>valval=10*s(1)^3+s(3)*(s(1)^2+s(2)^2);x=s;endend
end

三、函数引用与多变量

要分清多元函数的变量要用循环引用

代码如下(示例):

function v=testmain
x=-2:0.1:2;
v=fun(x);function r=fun(x)
for i=1:length(x)if x(i)<-1r(i)=x(i)^2+2*x(i)+1;elseif x(i)>-1 && x(i)<=2r(i)=3*x(i);elser(i)=2*sqrt(x(i));end
end

四、矩阵引用

引用矩阵可以用行列索引,也可以单行单列(注意单行单列是逗号),而且引用矩阵必须是括号

代码如下(示例):

A(1,:)%引用第一行
A(:,2)%引用第二列
A(1,:)=A(1,:).^2;%给矩阵赋值要明白给谁赋值

五、非线性函数:fmincon的使用

  1. 使用时要注意哪里是只输入系数,哪里是方程
  2. 求最小值的是方程,其他的都只输出系数
  3. 同种变量用逗号隔开,不同种用分号
  4. 注意noncon的位置是使用@,引用里面有()
  5. lb ub里面是逗号
  6. 函数系数的对应关系

在这里插入图片描述

代码如下(示例):

function [min,val]=fun3
Fun=@(x) 2*(x(1)-1)^2+3*(x(2)-1)^2+0.2*x(1)*x(2)+(2*x(3)-2)^2
A=[3 2 6;4 5 2;2 9 7];
b=[50;40;100];
Lb=[0,0,0];
ub=[15,9,25];
[min,val]=fmincon(Fun,[0,0,0],A,b,[],[],Lb,ub,@noncon);
function [c,ceq]=noncon(x)
c=[-2*x(1)^-x(2)^2-x(3)^2+5,x(1)^2+x(2)^2+2*x(3)^2-200];
ceq=[0];

六、线性规划函数

1.linprog

求解常规线性规划
注意!!!f 函数调用是系数

代码如下(示例):
在这里插入图片描述

function [x,fval]=myfun19
f=[3,-4,2];
A=[1 -1 1;3 2 4;3 2 0];
b=[10;22;10];
Aeq=[];
beq=[];
lb=[0,0,0];
ub=[inf,inf,inf];
[x,fval] = linprog(f,A,b,Aeq,beq,lb,ub);

2.fminbnd、fminsearch

用于找到极小值,fminbnd找到区间内,fminsearch找到从某一点开始的极小值
fminbnd

[x,val] = fminbnd(fun, x1, x2)

% 定义要最小化的函数
fun = @(x) x^2 + 2*x + 1;% 在区间[-2, 2]内寻找最小值
[x,val] = fminbnd(fun, -2, 2);

fminsearch

% 定义要最小化的函数
fun = @(x) x(1)^2 + x(2)^2;% 指定初始点
x0 = [1, 1];% 在多维空间中寻找最小值
x = fminsearch(fun, x0);

七、积分和微分

1.quad积分

y=@(x) (3.*x+2)./(x.^3+2.*x+5);
r=quad(y,0,3);

2.积分int 微分diff

在这里插入图片描述
在这里插入图片描述

F=@(x,y) exp(2*x.^2+3*y.^2);
I=integral2(F,0.5,1,0,1);%取值范围

八、lsqcurvefit

  1. xdata,ydata储存数据,x(i)存放未知变量
  2. 调用时要把 fun(x,xdata) 全部带上
  3. 注意调用时,x和x1的意义是已经计算的系数和要代入计算的函数值,他们的位置和使用fun时的位置相对应!!!

代码如下(示例):

xdata=[10.03,20.04,30.23,40.37,50.16,60.4,70.3,80.66,90.69,100.16,110.94,121.18,130.66,141.23,150.88,160.24];
ydata=[50.44,82.95,133.07,202.76,286.39,381.34,467.36,540.54,591.7,624.28,647.62,660.77,668.06,672.83,675.36,676.82];
fun = @(x,xdata) x(1)./(1+((x(1)/x(2))-1).*exp(-x(3)*xdata));
x0 = [1,1,1];
x = lsqcurvefit(fun,x0,xdata,ydata);
x1=[150:0.5:190];
y=fun(x,x1);

九、矩阵求解

  1. 可以使用linsolve
  2. 也可以求逆
  3. 或者是符号“\”
    在这里插入图片描述

代码如下(示例):

a=[1 1 -1 2;-1 -2 1 -3;1 1 1 -3;2 1 -3 1];
e=[-3;1;-7;4];
x=inv(a)*e;
a=[1 1 -1 2;-1 -2 1 -3;1 1 1 -3;2 1 -3 1];
e=[-3;1;-7;4];
x=linsolve(a,e);
a=[1 1 -1 2;-1 -2 1 -3;1 1 1 -3;2 1 -3 1];
e=[-3;1;-7;4];
x=a\e;

十、polyfit 函数进行多项式拟合

  1. p = polyfit(x, y, n);
  2. x 是数据点的 x 值(自变量)向量,y 是数据点的 y 值(因变量)向量,n 是要拟合的多项式的次数。
  3. polyval(p,50); 用于计算特定位置的值

代码如下(示例):

% 示例数据点
t=[1,6,11,16,21,26,31,36];
x=[5.35,10.51,15.664,20.76,25.94,30.08,35.22,40.36];
p=polyfit(t,x,1);
val=polyval(p,50);
r=[p,val];

十一、ode23和ode45

  1. 使用myode时,注意输入要把x,y都带上,先是自变量,后是因变量
  2. y(1)可表示y,y(2)可表示一阶导数
  3. 输出 [t,y] 分别是对应时间点和对应时间的输出值

在这里插入图片描述

y1 = 0:0.01:1;%取值范围
y0 = 1;初值
[t,y] = ode23(@myode, y1, y0);function dydx=myode(x,y)dydx=(2/3)*x/y^2;
function r = myfunx1 = 0:0.1:100;%定义范围x0 = [2; 0]; %定义初始条件[t, x] = ode23(@myode, x1, x0);%这里的x是一个向量
endfunction dxdt = myode(t, x)dxdt = [20*(1-x(1)^2)*x(2) + 0.5*x(1); 0];%x(1)代表x,x(2)代表dx/dt
end

十二、fzero

x0 是初始猜测的根的值。

代码如下(示例):

x = fzero(fun, x0)
x = fzero(fun, x0, options)

十三、meshgrid和surf

  1. 注意z的输出
  2. [x,y]=meshgrid(x取值,y取值)
  3. surf绘制表面图

代码如下(示例):

function z=myfun24
[x,y]=meshgrid(-3:0.1:3,-3:0.1:3);
z=x.*y.*exp(sin(x.^2+y.^2));
surf(x,y,z);

十四、fsolve

  1. 记住两个公式:root2和fsolve的使用
  2. folsve(@root2,x0)
  3. f=root2(x) f(1)= f(2)=

代码如下(示例):

% 示例数据点
function r=myfun29
x0=[0,0];
r=fsolve(@root2,x0);function f=root2(x)
f(1)=x(1)^2-x(2)-1;
f(2)=(x(1)-2)^2+(x(2)-0.5)^2-1;

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/diannao/42187.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

14-39 剑和诗人13 - 顶级大模型测试分析和建议

​​​​​ 随着对高级语言功能的需求不断飙升&#xff0c;市场上涌现出大量语言模型&#xff0c;每种模型都拥有独特的优势和功能。然而&#xff0c;驾驭这个错综复杂的生态系统可能是一项艰巨的任务&#xff0c;开发人员和研究人员经常面临选择最适合其特定需求的模型的挑战。…

哈弗架构和冯诺伊曼架构

文章目录 1. 计算机体系结构 2. 哈弗架构&#xff08;Harvard Architecture&#xff09; 3. 改进的哈弗架构 4. 冯诺伊曼架构&#xff08;Von Neumann Architecture&#xff09; 5. 结构对比 1. 计算机体系结构 计算机体系结构是指计算机系统的组织和实现方式&#xff0c…

Python | Leetcode Python题解之第220题存在重复元素III

题目&#xff1a; 题解&#xff1a; class Solution(object):def containsNearbyAlmostDuplicate(self, nums, k, t):from sortedcontainers import SortedSetst SortedSet()left, right 0, 0res 0while right < len(nums):if right - left > k:st.remove(nums[left]…

Python基础问题汇总

为什么学习Python&#xff1f; 易学易用&#xff1a;Python语法简洁清晰&#xff0c;易于学习。广泛的应用领域&#xff1a;适用于Web开发、数据科学、人工智能、自动化脚本等多种场景。强大的库支持&#xff1a;拥有丰富的第三方库&#xff0c;如NumPy、Pandas、TensorFlow等…

Sass 语法

文章目录 编译变量 \$嵌套 {} > \~导入 import注释 // /*\* \**/混入 mixin/include继承 extend数据类型运算控制 if/for/each/while函数 function媒体查询 media根发出 at-root警告warn/错误error/调试debug 编译 编译命令 单文件转换命令 sass input.scss output.css单…

数学基础 -- 反函数

反函数技术文档 反函数的定义 反函数&#xff08;inverse function&#xff09;是指一种将函数的输出反过来作为输入&#xff0c;从而恢复原来输入的函数。具体来说&#xff0c;如果有一个函数 f f f&#xff0c;它把一个值 x x x 映射到一个值 y y y&#xff0c;即 f ( …

68.WEB渗透测试-信息收集- WAF、框架组件识别(8)

免责声明&#xff1a;内容仅供学习参考&#xff0c;请合法利用知识&#xff0c;禁止进行违法犯罪活动&#xff01; 内容参考于&#xff1a; 易锦网校会员专享课 上一个内容&#xff1a;67.WEB渗透测试-信息收集- WAF、框架组件识别&#xff08;7&#xff09; 右边这些是waf的…

Mean teacher are better role models-论文笔记

论文笔记 资料 1.代码地址 2.论文地址 https://arxiv.org/pdf/1703.01780 3.数据集地址 CIFAR-10 https://www.cs.utoronto.ca/~kriz/cifar.html 论文摘要的翻译 最近提出的Temporal Ensembling方法在几个半监督学习基准中取得了最先进的结果。它维护每个训练样本的标签…

PCIe驱动开发(1)— 开发环境搭建

PCIe驱动开发&#xff08;1&#xff09;— 开发环境搭建 一、前言 二、Ubuntu安装 参考: VMware下Ubuntu18.04虚拟机的安装 三、QEMU安装 下载网站&#xff1a; https://download.qemu.org 下载文件&#xff1a;qemu-4.1.0-rc5.tar.xz 使用如下命令解压&#xff1a; tar …

opencv 设置超时时间

经常爬视频数据&#xff0c;然后用opencv做成图片 因此设置超时时间很重要 cap.set(cv2.CAP_PROP_FPS, timeout_ms) for idx, row in data.iterrows(): if idx < 400: continue try: # 打开视频文件 timeout_ms 5000 cap cv2.VideoCapture(row[PLAY_URL]) cap.set(cv2.C…

Linux下使用libiw进行无线信号扫描的实例

打开电脑连接wifi是一件很平常的事情,但这些事情通常都是操作系统下的wifi管理程序替我们完成的,如何在程序中扫描wifi信号其实资料并不多,前面已经有两篇文章介绍了如何使用ioctl()扫描wifi信号,但其实在Linux下有一个简单的库对这些ioctl()的操作进行了封装,这个库就是l…

深入追踪:IPython 中 %tb 命令的异常追踪栈使用指南

深入追踪&#xff1a;IPython 中 %tb 命令的异常追踪栈使用指南 在 IPython 的强大功能中&#xff0c;%tb 命令是一个调试工具&#xff0c;用于在出现异常时查看详细的异常追踪栈信息。这对于开发者来说是一个不可或缺的功能&#xff0c;因为它提供了对错误发生上下文的深入了…

Unity 中,常用的 UnityEngine.Events 中的几个重要的事件处理函数

在 Unity 中&#xff0c;常用的 UnityEngine.Events 中的几个重要的事件处理函数包括&#xff1a; UnityEvent UnityEvent 是 Unity 提供的一种事件系统&#xff0c;可以用来实现脚本与场景中的对象之间的互动。它可以用来定义和响应事件&#xff0c;如按钮点击、物体碰撞等。示…

GPT-5或重塑我们的工作与生活

引言 在人工智能发展的浪潮中&#xff0c;每一次技术的革新都如同潮水般涌来&#xff0c;带来前所未有的机遇与挑战。当新一代大语言模型GPT-5即将登场的消息传来&#xff0c;我们不禁要问&#xff1a;它将如何重塑我们的工作和日常生活&#xff1f;又将开启哪些崭新的应用场景…

故障模式与影响分析(FMEA)的概念

故障模式与影响分析&#xff08;FMEA&#xff09;的概念 故障模式与影响分析&#xff08;Failure Mode and Effects Analysis&#xff0c;FMEA&#xff09;是一种系统性评估方法&#xff0c;用于识别产品设计或过程中可能发生的潜在故障模式&#xff0c;以及这些故障模式对系统…

制作爬取4399游戏名称软件

def 爬取4399(): #发送请求并且拿到源代码 import requests 链接https://www.4399.com/ #网站链接 请求头{User-Agent:Mozilla/5.0 (Windows NT 10.0; Win64; x64; rv:127.0) Gecko/20100101 Firefox/127.0} #构造请求头用于爬取网站源代码使用 网站源代码…

MySQL远程登录

root是超级管理员&#xff0c;默认情况下&#xff0c;root不能作为远程登录的用户名&#xff0c;远程登录前&#xff0c;需要将登录的数据库在本地登录&#xff0c;修改权限&#xff0c;输入&#xff1a; update user set host & where user root ; 回车键&#xff0c…

clickhouse高可用可拓展部署

clickhouse高可用&可拓展部署 1.部署架构 1.1高可用架构 1.2硬件资源 部署服务 节点名称 节点ip 核数 内存 磁盘 zookeeper zk-01 / 4c 8G 100G zk-02 / 4c 8G 100G zk-03 / 4c 8G 100G clikehouse ck-01 / 32c 128G 2T ck-02 / 32c 128G 2T ck-03 / 32c 128G 2T ck-04 /…

[Qt] 控件的QSizePolicy属性选项

在Qt中&#xff0c;QSizePolicy是一个非常重要的枚举类&#xff0c;它定义了控件&#xff08;widgets&#xff09;在布局管理&#xff08;layout management&#xff09;中的大小调整策略。这些策略决定了当控件的父布局或窗口大小发生变化时&#xff0c;控件应该如何调整自身的…

java wait, notify, notifyAll三个方法

wait(), notify(), 和 notifyAll() 是 Java 中用于线程间通信和同步的方法&#xff0c;它们都是 Object 类中的方法&#xff0c;而非 Thread 类的方法。这些方法通常与 synchronized 关键字一起使用&#xff0c;用于实现线程之间的协作和互斥访问共享资源。 关于生产者-消…