(转)matlab各类数学公式

  matlab矩阵应用
clear
%建立矩阵的两种方式
A1 = [1 2 3 4 5; 6 7 8 9 10];
A2 = [
        1 2 3 4 5
        6 7 8 9 10
      ];
%一种是换行用引号,一种是自然写法

clear
A = [5 4 3 2 1; 6 7 8 9 10;1 2 3 4 5;24 24 24 24 24;25 25 25 25 25];
B = [5 4 3 2 1; 6 7 8 9 10;1 2 3 4 5;24 24 24 24 24;25 25 25 25 25];
A_sqare = [1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25];
k = 5;

 

X1 = A’;   %A取转置
X2 = A + B; %求A和B矩阵的和
X3 = A – B; %求A和B矩阵的差
X4 = k * A; %数K乘以A矩阵
X5 = det(A_sqare); %求矩阵A的行列式(注A_sqare必须为方阵)
X6 = rank(A); %求矩阵A的秩
X7 = inv(A_sqare); %求矩阵A的逆(注A_sqare必须为方阵)
X8 = B / A; X8 = B * inv(A_sqare);  %A右除B = B右乘A的逆
X9 = B / A; X9 = inv(A_sqare) * B;  %A左除B = B左乘A的逆
A.*B    % .* 是A的每个元素和B相乘,非矩阵相乘,同理  ./ .^ 
A(2,:)   %  取A矩阵的第2行 然后 A(2,:) = [5 5 5 5 5]就对该行进行赋值
A(:,2)   %  取A矩阵的第2列 然后 A(:,2) = [5 5 5 5 5]就对该列进行赋值
A(2:1:4,3:1:4)  %取A矩阵中的一块,其语法为A(起始行:步长:终止行,起始列:步长:终止列)
zeros(5)   %生成n阶零矩阵
eye(5)     %生成n阶单位矩阵
eig(A)     %矩阵A的特征值
[X,D] = eig(A)   %矩阵A的  特征向量矩阵X  特征值组成的对角阵
A([1,2],:)   %1,2行互换
A(:,[2,3])   %1,2列互换
A(2,:) = 5 * A(2,:)  %第2行乘以5,列上的操作以此类推
K = [A B;B A]     %由几个小矩阵合成一个大矩阵
orth(A)      %非奇异矩阵正交化
a1 = A(2,:);
a2 = A(3,:);
a1*a2′       %两个向量内积
rref(A)     %A的极大无关向量组
                     matlab画图应用
%while循环
clear
sum = 0;k = 1;
while  k<101
  sum = sum+k;
  k = k + 1;  
end
sum

%for循环
clear
sum = 0;n = 1;
for n = 1:100
    sum = n + sum;
end
sum

%plot绘图
x = -10:0.1:10;
y = 3*x.^4+x.^2-1;
figure            %开启新绘图页面
plot(x,y)

%fplot绘图
clear
figure
fplot(@f1, [-10 10])     %调用f1.m里面的那个function

%ezplot绘图
clear
syms x
figure
y = 3*x^4 + x^2 – 1;
ezplot(y)

%螺旋线绘图  题目:x = cost , y = sint , z = t t属于[0,6pi]
clear
t = 0:0.1:6*pi;
x = cos(t);
y = sin(t);
z = t;
figure
plot3(x,y,z)

%空间曲线绘图  题目: z = sqrt(1 – x^2 – y^2) , (x – 1/2)^2 + y^2 = (1/2)^2
clear
t = 0:0.1:6;
x = 0.5*cos(t)+0.5;
y = 0.5*sin(t);
z = sqrt(1-x.^2-y.^2);
figure
plot3(x,y,z)

%二次曲面绘图  题目: x^2 + y^2 = z
clear
s = -10:1:10;
t = -10:1:10;
[x,y] = meshgrid(s,t);
z = x.^2 + y.^2;
figure
mesh(x,y,z);

%旋转曲面绘图  题目: y = 1/x  围绕y轴旋转
clear
s = -10:0.1:10;
t = -10:0.1:10;
[x,y] = meshgrid(s,t);
r = 1./x;
[x,y,z] = cylinder(r);
figure
mesh(x,y,z)

%输入数
clear
K = input(‘请输入数’);    
                      matlab解方程应用
%常微分方程求解     详细参考PPT(7)
clear
y1 = dsolve(‘Dy=8-3*y’,'y(0)=2′)
y2 = dsolve(‘D2y=2*x*Dy/(1+x^2)’,'y(0)=1,Dy(0)=3′)

%微分方程组求解
clear
[X,Y] = dsolve(’2*Dx+4*x+Dy-y=exp(t),Dx+3*x+y=0′,’x(0)=1.5,y(0)=0′)

%线性方程组全部解
clear
format rat
A=[1 1 3 -1;0 1 -1 1;1 1 2 2;1 -1 1 -1];
B=[-2;1;4;0];
X = A/B

%方程求解
clear
X = solve(‘x-exp(-x)=0′,’x')

%区间里方程求解
%x=0:0.1:10;
%X = solve(’5*x^2*sin(x)-exp(-x)’,'x’)

%求微分方程的特解并且做出函数曲线
y0 = [1,0];
[t,x] = ode45(@vdp,[0,30],y0);   %从vdp.m这个文件里面读函数
y = x(:,1);
dy=x(:,2);
figure
plot(t,y,t,dy);

%解微分方程 
fun=inline(‘-2*y+2*x*x+2*x’);
[x,y]=ode23(fun,[0,0.5],1)
                     matlab微积分应用
                     %函数求导
clear
syms x y;     %将x y设为变量
f = cos(x)^3-cos(3*x);    %需要求导的式子
dy = diff(f,x);       %用diff 函数进行求导

%函数求极限
clear
syms x
f = x*log(1+x)/sin(x*x)
limit(f,’x',0,’left’)    %语法为  limit(求极限的式子,求极限的变量,需要逼近的数字,从左还是从右逼近)

%函数求积分
clear
syms x
f = sym(‘x*exp(x)/(1+x)^2′);
int(f)       %语法为  int(求极限的式子,下限,上限)  如果不加上限下限,就是函数式。
int(f,0,1)  
pretty(f)    %以自然函数形式呈现
%例-求z = x^2 + y^2 , z = 1 , z =2 围成的曲面
clear
syms x y z
z = x^2+y^2;
f = z;
I = int( int( f, y, sqrt(1-x^2), sqrt(2-x^2) ), x, 1, sqrt(2) )

%级数求和
clear
syms n
f = (n+1)/n*2^n;
j = symsum(f, n, 1, inf)    %级数求和,下限为1,上限为无穷大

%泰勒展开
clear
syms x
f = cos(x)
taylor(f, 10, x, pi/3)    %语法为taylor(待展开函数,取前几项, 变量名, 展开中心)

%求傅里叶系数
clear
syms x n
f = x^3+x^2;
n = 5;
a0 = int(f,x,-pi,pi)/pi
a1 = int(f*cos(1*x),x,-pi,pi)/pi
a2 = int(f*cos(2*x),x,-pi,pi)/pi
a3 = int(f*cos(3*x),x,-pi,pi)/pi
a4 = int(f*cos(4*x),x,-pi,pi)/pi
a5 = int(f*cos(5*x),x,-pi,pi)/pi
b1 = int(f*sin(1*x),x,-pi,pi)/pi
b2 = int(f*sin(2*x),x,-pi,pi)/pi
b3 = int(f*sin(3*x),x,-pi,pi)/pi
b4 = int(f*sin(4*x),x,-pi,pi)/pi
b5 = int(f*sin(5*x),x,-pi,pi)/pi

 

转载于:https://www.cnblogs.com/sallybin/p/3318585.html

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

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

相关文章

PS网页设计教程XVI——在PS中创建一个摩登实验室风格的网页设计

作为编码者&#xff0c;美工基础是偏弱的。我们可以参考一些成熟的网页PS教程&#xff0c;提高自身的设计能力。套用一句话&#xff0c;“熟读唐诗三百首&#xff0c;不会作诗也会吟”。 本系列的教程来源于网上的PS教程&#xff0c;都是国外的&#xff0c;全英文的。本人尝试翻…

好评率超高的9个公众号,值得收藏

不管世界变化多快&#xff0c;阅读仍是个不过时的动作。今天就为大家推荐几个优质公众号&#xff0c;从社会热点、读书成长、生活新知&#xff0c;你的碎片化时间&#xff0c;我们包了&#xff01;一起来看看呀~长按二维码&#xff0c;“识别图中二维码”即可订阅————电商头…

.NET 6新特性试用 | 文件范围的命名空间

前言在前面的《隐式using指令》文章中&#xff0c;我们介绍了global using&#xff0c;可以减少我们代码文件中的“噪音”&#xff0c;使代码更简洁。今天&#xff0c;我们介绍.NET 6中另一个简化代码的方式——文件范围的命名空间。Demo当我们向项目中添加一个类文件时&#x…

php的几种运行模式CLI、CGI、FastCGI、mod_php

1、CLI:就是命令行&#xff0c;例如可以在控制台或者是shell中键入命令: php -f index.php 然后获取输出 2、CGI:以下是不同的说法与理解 公共网关接口”(Common Gateway Interface)&#xff0c;HTTP服务器 与你的或其它机器上的程序 进行 “交谈”的一种工具 &#xff0c;其程…

Android之三大图片缓存原理、特性对比

一. 四大图片缓存基本信息 Universal ImageLoader 是很早开源的图片缓存&#xff0c;在早期被很多应用使用。Picasso 是 Square 开源的项目&#xff0c;且他的主导者是 JakeWharton&#xff0c;所以广为人知。Glide 是 Google 员工的开源项目&#xff0c;被一些 Google App 使用…

android宿舍管理系统源码,基于android操作系统的手机宿舍管理系统使用手册

《基于android操作系统的手机宿舍管理系统使用手册》由会员分享&#xff0c;可在线阅读&#xff0c;更多相关《基于android操作系统的手机宿舍管理系统使用手册(5页珍藏版)》请在人人文库网上搜索。1、基于an droid操作系统的手机宿舍管理系统使用手册 1 系统概述 本系统是在an…

产品需求文档 PRD

第一轮&#xff1a; 1&#xff0c;文档使用方&#xff1a;UI设计师 2、内容&#xff1a; 根据战略层定义出来产品功能范围&#xff0c; 说明此产品的目的&#xff0c;方便UI设计人员更好的理解产品 产品基本流程 详细的设计框架图&#xff0c;推荐用ax…

他读博期间连发3篇Science,28岁任武大教授后再发Nature!

全世界只有3.14 % 的人关注了爆炸吧知识来源&#xff1a;募格课堂整合自武汉大学、武汉晚报、半月谈等有这样一位学者&#xff0c;读博期间以一作身份发表3篇Science,28岁成为武汉大学化学与分子科学学院教授再发2篇Nature&#xff0c;35岁荣获杰青......今年9月&#xff0c;他…

如何限制并发的 异步IO 请求数量?

咨询区 Grief Coder&#xff1a;我的项目中有下面一段代码&#xff1a;// lets say there is a list of 1000 URLs string[] urls { "http://google.com", "http://yahoo.com", ... };// now lets send HTTP requests to each of these URLs in parallel …

Android之自定义 ActionBar 上的菜单(Menu)文字颜色

本文记录了修改 ActionBar 上 Menu 字体颜色的方法。第一种通过自定义主题实现&#xff0c;第二种通过复写onCreateOptionsMenu实现。 自定义Theme 修改MenuItem 文字颜色 自定义 ActionBar 上的菜单文字颜色需要通过 Activity 的 Theme 实现&#xff0c;而非 ActionBar 的 The…

删除远程桌面登陆痕迹

怎么删除本地远程桌面地址记录(连接地址输入框的下拉菜单)。 首先运行regedit&#xff0c;打开注册表&#xff0c;找到 HKEY_CURRENT_USER&#xff3c;Software&#xff3c;Microsoft&#xff3c;Terminal Server Client&#xff3c;Default 右边的那些键值就记录了IP或者域名&…

Java Web开发Tomcat中三种部署项目的方法

第一种方法&#xff1a;在tomcat中的conf目录中&#xff0c;在server.xml中的&#xff0c;<host/>节点中添加&#xff1a; <Context path"/hello" docBase"D:\eclipse3.2.2forwebtools\workspace\hello\WebRoot" debug"0" privileged&q…

Android 页面布局xd,Adobe XD强大的布局系列工具 助你事半功倍

数字设计有多种形式-从移动端、桌面网站到显示屏应用程序、智能助手、车载娱乐系统(例如CarPlay和Android Auto)等等。为当今的应用程序打造数字体验意味着需要针对各种屏幕尺寸、布局和交互类型进行设计。作为设计师&#xff0c;您的时间应该花在用户体验上&#xff0c;而不是…

JS中数组Array的用法{转载}

js数组元素的添加和删除一直比较迷惑&#xff0c;今天终于找到详细说明的资料了&#xff0c;先给个我测试的代码^-^var arr new Array();arr[0] "aaa";arr[1] "bbb";arr[2] "ccc";//alert(arr.length);//3arr.pop();//alert(arr.length);//2…

如果觉得午休时间太短怎么办?

1 假装自己戴了口罩&#xff08;素材来源网络&#xff0c;侵删&#xff09;▼2 正义不会迟到博主是男的&#xff08;素材来源网络&#xff0c;侵删&#xff09;▼3 阅读理解有多难&#xff1f;&#xff08;素材来源网络&#xff0c;侵删&#xff09;▼4 下雨天就不要爬树了&…

Android之ActionBar常用设计和使用总结

常用ActionBar总结 我们经常看见项目的顶部有一些文字显示在中间,或者左边、右边,而且还有一个返回的图标,很方便我们返回,接下来我总结下常用的ActionBar 第一步:写常见的ActionBar的布局文件 我这里写了3个actionbar文件 1、actionbar_return.xml …

WPF 实现温度计

WPF开发者QQ群&#xff1a; 340500857 | 微信群 -> 进入公众号主页 加入组织在WPF中没有现成的温度计控件&#xff0c;所以我们自己实现一个。微信群人数太多入群请添加小编微信号&#xff08;yanjinhuawechat&#xff09;或&#xff08;W_Feng_aiQ&#xff09;邀请入群&am…

100-6

2019独角兽企业重金招聘Python工程师标准>>> /**\第6题&#xff08;数组&#xff09;腾讯面试题&#xff1a; 给你10分钟时间&#xff0c;根据上排给出十个数&#xff0c;在其下排填出对应的十个数 要求下排每个数都是先前上排那十个数在下排出现的次数。 上排的十个…

linux系统644、755、777权限详解

在linux系统中644、755、777三种权限是非常重要的一些权限了&#xff0c;下面我来详细的介绍644、755、777三种权限的使用&#xff0c;希望对各位有帮助。常用的linux文件权限&#xff1a; 444 r--r--r-- 600 rw------- 644 rw-r--r-- 666 rw-rw-rw- 700 rwx------ 744 rwxr--r…

android获取自适应高度,Android中oncreate中获得控件高度或宽度的实现方法

Android中oncreate中获得控件高度或宽度的实现方法onCreate函数只是提供了数据初始化的机会&#xff0c;此时还没有正式绘制图形。在图形尚未渲染的情况下&#xff0c;getWidth()或getHeight()得到的值为0。如果需要自适应屏幕&#xff0c;需要一进入界面就得到父控件大小&…