【MATLAB实验】数学实验实验求根问题三种方法、积分画图及古典密码设计(二分法、牛顿法、不动点迭代法)
目录
特征值与特征向量
求最简型
编程题
求积分问题(quad、trapz、int)
求根问题(二分法、不动点迭代法、牛顿法)
1.二分法:
2.不动点迭代法
3.牛顿切线法
画图题
画螺旋曲线图
画椭圆抛物面图
画马鞍面图
画三角曲面图
密码问题(古典密码设计)
编辑
微分方程数值解问题
特征值与特征向量
求矩阵A=[5 0 0;0 3 4;0 4 3]的特征值和特征向量。
矩阵特征值的数学定义:设A是n阶方阵,如果存在常数λ和n维非零列向量x,使得等式Ax=λ x成立,则称λ为A的特征值,x是对应特征值λ的特征向量。
[X,D]=eig(A):求矩阵A的全部特征值,构成对角阵D,并产生矩阵X,X各列是相应的特征向量。
求最简型
求矩阵A=[1 9 2 8;1 9 3 1;1 9 3 2;1 9 3 7;1 9 3 8;1 9 4 5]的行最简型。
编程题
1.编函数文件 y=sushu(n),若n为素数返回y的值为1,否则,返回y的值为0。并在命令窗口,调用y=sushu(33)。
function y=sushu(n)y=1;for i=2:n-1if mod(n,i) == 0y=0;endend
end
2.编写函数文件“sushu_n.m”,使得在命令窗口调用“y=sushu n(n)”求小于n的所有素数,即函数返回值y是一个矩阵,它的没个元素都是小于n的素数,并在命令窗口调用y=sushu n(20)。
function y=sushu_n(n)y = [];for i = 2:n - 1flag = 0;for j = 2:i - 1if mod(i, j) == 0flag = 1;endendif flag == 0y = [y, i];end
end
3.求小于等于n的数的阶乘之和。并在窗口调用y=jiecheng(10)。(如y=jiecheng(4)=1!+2!+3!+4!)
function y=jiecheng2(n)
y=0;
for i=1:nx=1;for j=1:ix=x*j;endy=y+x;
end
求积分问题(quad、trapz、int)
1.用quad计算的值。
2.用trapz计算的值。
3. 分别用梯形法与抛物线法,计算,将积分区间[1,2]作120等分,并尝试直接使用函数trapz(),quad()进行计算求解,比较结果的差异。
%trapz梯形法
x=1:2/120:2;
y=1./x;
trapz(x,y)
ans=0.6858%quad辛普森法
quad('1./x',1,2)
ans =0.6931%符号求积分
int('1/x','x',1,2)
ans =log(2)
我们发现梯形法和辛普森法都能够计算出结果,结果相近似但与真实的值有差异,有一定差异,而符号法求出的积分为log(2)为真实的值。
4.试计算定积分。(注意:可以使用trapz()和quad或者附录程序求解吗?)
注意这里不能够使用trapz()函数和quad函数 。
求根问题(二分法、不动点迭代法、牛顿法)
1.二分法:
用二分法求在[2,3]内一个近似根,满足。
clear,clc
a=2;b=3;x=0.5*(a+b);
f=@(x)x^3-x-21;
while abs(f(x)) > 0.000001if f(a)*f(x)<0b=x;x=0.5*(a+b);elsea=x;x=0.5*(a+b);end
end
disp('the root is'),x=x
str1=sprintf('f(x)=%d',f(x));
disp(str1)
2.不动点迭代法
求在[0,2]上的根。(fun函数是将x放到一边,如本题)
clear,clc
f=@(x)x^5+9*x-10;
fun=@(x)(10-x^5)/9;
x=1.5;
while abs(f(x))>0.000001x=fun(x);
end
fprintf('x=%f,f(x)=%f\n',x,f(x))
3.牛顿切线法
用牛顿法求方程在[1,1.5]内的实根,取迭代的初始值。
x=1.5;k=0;
f=@(x)x^3+5*x^2-15;
df=@(x)3*x^2+10*x;
[k,x,f(x)]
while abs(f(x))>0.000001x=x-f(x)/df(x);k=k+1;[k,x,f(x)]
end
画图题
画螺旋曲线图
绘制简单的螺旋曲线。
画出空间曲线 。
画椭圆抛物面图
用函数surf画出表示的图形。
画马鞍面图
1.用mesh函数画出表示的图形。
2.用mesh函数画出表示的图形。
画三角曲面图
绘制 的带等高线的三维曲面图(surfc,x,y的取值范围为[0,2pi])。
密码问题(古典密码设计)
1.
2.
微分方程数值解问题
【MATLAB实验】数学实验实验求根问题三种方法、积分画图及古典密码设计(二分法、牛顿法、不动点迭代法)相关推荐
- C语言求最大公约数三种方法详解
C语言求最大公约数三种方法详解 题目要求 常用写法(穷举法) 辗转相减法 辗转相除法 main函数 整体代码 题目要求 运行最大公约数的常用算法,并进行程序的调式与测试. 常用写法(穷举法) 从两个数 ...
- 求最小公倍数的三种方法(C语言)
求最小公倍数的三种方法 1.常规暴力求解法 2.辗转相除法 3.迭乘法 //1.常规暴力求解法 #include <stdio.h> int main() {int a = 0;int b ...
- java中隐函数求导法则_隐函数求导的三种方法
this.p={ m:2, b:2, loftPermalink:'', id:'fks_085075084086088070081083074065081087082066093087080', b ...
- python求平方根的三种方法
python求平方根的三种方法 题干描述 题目解答 题干描述 没啥好说的qwq,求根号下x,并舍弃小数部分,只保留整数 题目解答 方法一:不多bb,直接0.5次方(这应该是最没有营养的解法,面试官估计 ...
- C语言求幂的三种方法
用三种方法求幂值 一. 暴力递归 直接对x乘y次 int result(int x,int y) {int num=1;for (int i=1; i<=y; i++) {num*=x;}ret ...
- c语言程序π,C语言求圆周率π(三种方法)
题目1) 利用公式①计求π的近似值,要求累加到最后一项小于10^(-6)为止. 题目2) 根据公式②,用前100项之积计算π的值. 题目1)提供了一种解法,题目2)提供了两种解法,请看解析. 题目1) ...
- c语言求圆周率 . 4,C语言求圆周率π(三种方法)(4页)-原创力文档
C语言求圆周率π(三种方法) 题目1) 利用公式①计求π的近似值,要求累加到最后一项小于10^(-6)为止.题目2) 根据公式②,用前100项之积计算π的值.题目1)提供了一种解法,题目2)提供了两种 ...
- C语言求最小公倍数的三种方法
第一种方法:累加法求 思路是,两个数字,要求他们的最小公倍数,那么这个最小公倍数,至少不要比这两个要求的数小.我们首先判断出两个数中较大的一个,然后判断这个数是否是要求的两个数的最小公倍数.如果不是, ...
- 求 LCA 的三种方法
(YYL: LCA 有三种求法, 你们都知道么?) (众神犇: 这哪里来的傻叉...) 1. 树上倍增 对于求 LCA, 最朴素的方法是"让两个点一起往上爬, 直到相遇", &qu ...
最新文章
- 关于学习Python的一点学习总结(31->继承及多态)
- Jmeter之创建Kafka生产者和消费者进行性能测试
- django(7)modelform操作及验证、ajax操作普通表单数据提交、文件上传、富文本框基本使用...
- 入门设计模式之建造者模式
- mysql global index_Oracle中addsplit partition对globallocal index的影响
- Linux epoll 笔记(高并发事件处理机制)
- 我最喜欢的模板jade(pug)学习和使用
- leetcode题库:3.无重复字符的最长子串
- java spin lock_JAVA 各种锁机制
- oracle数据库stuff的用法,sql STUFF用法
- spdep | 除了莫兰指数,还有哪些指数可以衡量空间自相关性?
- Elasticsearch--Docker安装kibana---全文检索引擎ElasticSearch工作笔记003
- agilent3070软件安装
- Hive里的分区、分桶、视图和索引再谈
- php木马在线查杀_网站木马监控查杀扫描 ASP木马 PHP木马上传拦截
- c语言编程被7整除的4位数,四位数3a8c能被7整除,且是55的倍数,求a..._一建考试_帮考网...
- css设定元素左对齐,左对齐标签 - 右对齐选择元素(CSS)
- 计算机常用软件英文读音,学习英文在线发音的软件有哪些?宝妈推荐的这些超级实用...
- Java List 过滤重复数据
- Zero-Shot Transfer Learning for Event Extraction