台湾国立大学郭彦甫Matlab教程笔记(16) 数值微分 numerical differentiation

本文主要是介绍台湾国立大学郭彦甫Matlab教程笔记(16) 数值微分 numerical differentiation,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

台湾国立大学郭彦甫Matlab教程笔记(16) 数值微分 numeric differentiation

复习:diff()函数用来计算vector前后 entry的差异

数值微分继续

various step size 不同Δx大小

不同Δx大小,看近似值和真正函数值之间的差异大不大。

看题目,用不同的h值计算函数f(x)=sin(x)的导数,画出图形

The derivatives(导数) of f(x)=sin(x) calculated using various h values

codes:

g=colormap(lines);hold on;
for i=1:4%为了给出不同的h值x=0:power(10,-i):pi;%power(10,-i)表示10的-i次方
y=sin(x);
m=diff(y)./diff(x);%算是计算导数
plot(x(1:end-1),m,'color',g(i,:));%x(1:end-1)因为diff差值,vector维度少一
end
hold off;set(gca,'xlim',[0,pi/2]);%设置x轴横坐标范围
set(gca,'ylim',[0,1.2]);
set(gca,'xtick',0:pi/4:pi/2);%设置x轴坐标间隔
xticklabels({'0','\pi/4','\pi/2'});%显示pi
set(gca,'fontsize',18);%设置字体
h=legend('h=0.1','h=0.01','h=0.001','h=0.0001');%设置图像名称
set(h,'fontname','Times New Roman');
box on;%右上边框显示

代码的解释
colormap回顾:colormap输入参数的用法
在这里插入图片描述在这里插入图片描述

plot(x(1:end-1),m,‘color’,g(i,:));的解释:
%x(1:end-1)因为diff差值,vector维度少一,所以需要end-1。这里color用的是g(i,:),由于g是colormap中的lines颜色图,颜色图是一个矩阵matrix,这里应该是二维的,g(i,:)取的是这个矩阵的第i行元素,我的理解是一个行向量,然后赋值给了color显示出来。

运行结果:

在这里插入图片描述
作业题:
给定一个函数,用不同的h值计算近似导数

在这里插入图片描述
我的代码:

g=colormap(lines);%颜色
hold on;
%先来确定x的范围:
for i=1:3x=0:power(10,-i):2*pi;%不同的h 步长y=exp(-x).*sin(x.^2/2);%函数表达式需要点乘m=diff(y)./diff(x);%差分,求导plot(x(1:end-1),m,'color',g(i,:));
end
hold off;%下面是一些设置
box on;
set(gca,'fontsize',18);
set(gca,'xlim',[0,2*pi]);%x轴范围
set(gca,'ylim',[-0.3,0.3]);%y轴范围
set(gca,'xtick',0:pi/2:2*pi);%刻度
xticklabels({'0','\pi/2','\pi','3\pi/2','2\pi'});%显示刻度
legend('h=0.1','h=0.01','h=0.001');%函数标头

有个问题:这个函数表达式是两个函数相乘的形式,中间是用*还是用.*点乘?
实验结果:
y=exp(-x).*sin(x^2/2);(错误)
在这里插入图片描述
计算这个表达式正确的形式

y=exp(-x).*sin(x.^2/2);

我的练习结果:

在这里插入图片描述
下一个题目是

二次和三次微分second and third derivatives

the second derivative f’’ and third derivative f’’’ can be obtained using similar approaches

举例子:
在这里插入图片描述

例程代码:

x=-2:0.005:2;
y=x.^3;
m=diff(y)./diff(x);%一次微分
m2=diff(m)./diff(x(1:end-1));%二次微分plot(x,y,x(1:end-1),m,x(1:end-2),m2);%作图
xlabel('x','fontsize',18);
ylabel('y','fontsize',18);
legend('f(x)=x^3','f''(x)','f'''(x)',4);%函数标头 
set(gca,'fontsize',18);

【注意】每做一次微分,entry就少一个。

运行结果:
在这里插入图片描述

总结一下:
数值微分主要使用diff()函数。需要注意的是做一次微分,entry减少一,就是vector维度少一个。
一次微分:m=diff(y)./diff(x);%一次微分
二次微分:m2=diff(m)./diff(x(1:end-1));%二次微分

这篇关于台湾国立大学郭彦甫Matlab教程笔记(16) 数值微分 numerical differentiation的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/466071

相关文章

Spring Security 从入门到进阶系列教程

Spring Security 入门系列 《保护 Web 应用的安全》 《Spring-Security-入门(一):登录与退出》 《Spring-Security-入门(二):基于数据库验证》 《Spring-Security-入门(三):密码加密》 《Spring-Security-入门(四):自定义-Filter》 《Spring-Security-入门(五):在 Sprin

Makefile简明使用教程

文章目录 规则makefile文件的基本语法:加在命令前的特殊符号:.PHONY伪目标: Makefilev1 直观写法v2 加上中间过程v3 伪目标v4 变量 make 选项-f-n-C Make 是一种流行的构建工具,常用于将源代码转换成可执行文件或者其他形式的输出文件(如库文件、文档等)。Make 可以自动化地执行编译、链接等一系列操作。 规则 makefile文件

【学习笔记】 陈强-机器学习-Python-Ch15 人工神经网络(1)sklearn

系列文章目录 监督学习:参数方法 【学习笔记】 陈强-机器学习-Python-Ch4 线性回归 【学习笔记】 陈强-机器学习-Python-Ch5 逻辑回归 【课后题练习】 陈强-机器学习-Python-Ch5 逻辑回归(SAheart.csv) 【学习笔记】 陈强-机器学习-Python-Ch6 多项逻辑回归 【学习笔记 及 课后题练习】 陈强-机器学习-Python-Ch7 判别分析 【学

系统架构师考试学习笔记第三篇——架构设计高级知识(20)通信系统架构设计理论与实践

本章知识考点:         第20课时主要学习通信系统架构设计的理论和工作中的实践。根据新版考试大纲,本课时知识点会涉及案例分析题(25分),而在历年考试中,案例题对该部分内容的考查并不多,虽在综合知识选择题目中经常考查,但分值也不高。本课时内容侧重于对知识点的记忆和理解,按照以往的出题规律,通信系统架构设计基础知识点多来源于教材内的基础网络设备、网络架构和教材外最新时事热点技术。本课时知识

SWAP作物生长模型安装教程、数据制备、敏感性分析、气候变化影响、R模型敏感性分析与贝叶斯优化、Fortran源代码分析、气候数据降尺度与变化影响分析

查看原文>>>全流程SWAP农业模型数据制备、敏感性分析及气候变化影响实践技术应用 SWAP模型是由荷兰瓦赫宁根大学开发的先进农作物模型,它综合考虑了土壤-水分-大气以及植被间的相互作用;是一种描述作物生长过程的一种机理性作物生长模型。它不但运用Richard方程,使其能够精确的模拟土壤中水分的运动,而且耦合了WOFOST作物模型使作物的生长描述更为科学。 本文让更多的科研人员和农业工作者

论文阅读笔记: Segment Anything

文章目录 Segment Anything摘要引言任务模型数据引擎数据集负责任的人工智能 Segment Anything Model图像编码器提示编码器mask解码器解决歧义损失和训练 Segment Anything 论文地址: https://arxiv.org/abs/2304.02643 代码地址:https://github.com/facebookresear

沁恒CH32在MounRiver Studio上环境配置以及使用详细教程

目录 1.  RISC-V简介 2.  CPU架构现状 3.  MounRiver Studio软件下载 4.  MounRiver Studio软件安装 5.  MounRiver Studio软件介绍 6.  创建工程 7.  编译代码 1.  RISC-V简介         RISC就是精简指令集计算机(Reduced Instruction SetCom

数学建模笔记—— 非线性规划

数学建模笔记—— 非线性规划 非线性规划1. 模型原理1.1 非线性规划的标准型1.2 非线性规划求解的Matlab函数 2. 典型例题3. matlab代码求解3.1 例1 一个简单示例3.2 例2 选址问题1. 第一问 线性规划2. 第二问 非线性规划 非线性规划 非线性规划是一种求解目标函数或约束条件中有一个或几个非线性函数的最优化问题的方法。运筹学的一个重要分支。2

【C++学习笔记 20】C++中的智能指针

智能指针的功能 在上一篇笔记提到了在栈和堆上创建变量的区别,使用new关键字创建变量时,需要搭配delete关键字销毁变量。而智能指针的作用就是调用new分配内存时,不必自己去调用delete,甚至不用调用new。 智能指针实际上就是对原始指针的包装。 unique_ptr 最简单的智能指针,是一种作用域指针,意思是当指针超出该作用域时,会自动调用delete。它名为unique的原因是这个

查看提交历史 —— Git 学习笔记 11

查看提交历史 查看提交历史 不带任何选项的git log-p选项--stat 选项--pretty=oneline选项--pretty=format选项git log常用选项列表参考资料 在提交了若干更新,又或者克隆了某个项目之后,你也许想回顾下提交历史。 完成这个任务最简单而又有效的 工具是 git log 命令。 接下来的例子会用一个用于演示的 simplegit