零阶一类贝塞尔函数曲线及其包络线 matlab源程序
废话不多说,请看源程序:
clear ,clc;
format long
x=(0:0.01:100)';
y_0=besselj(0,x);
% y_1=besselj(1,x); %一阶,这里只画了0阶
% y_2=besselj(2,x); %二阶
plot(x,y_0);grid on;
axis([0,100,-1,1]);
title('0阶贝塞尔函数曲线图');
xlabel('Variable X');
ylabel('Variable Y');
%画包络线
hold on;
[up,down] = envelope(x,y_0,'spline');
plot(x, up, 'r');
plot(x, down, 'r');
其中的envelope是个单独的.m文件,是个函数,用来画包络线。说白了就是求拐点然后插值:
function [up,down] = envelope(x,y,interpMethod)
%ENVELOPE gets the data of upper and down envelope of the known input (x,y).
%
% Input parameters:
% x the abscissa of the given data
% y the ordinate of the given data
% interpMethod the interpolation method
%
% Output parameters:
% up the upper envelope, which has the same length as x.
% down the down envelope, which has the same length as x.
%
% See also DIFF INTERP1
% Designed by: Lei Wang, <WangLeiBox@hotmail.com>, 11-Mar-2003.
% Last Revision: 21-Mar-2003.
% Dept. Mechanical & Aerospace Engineering, NC State University.
% $Revision: 1.1 $ $Date: 3/21/2003 10:33 AM $
if length(x) ~= length(y)
error('Two input data should have the same length.');
end
if (nargin < 2)|(nargin > 3),
error('Please see help for INPUT DATA.');
elseif (nargin == 2)
interpMethod = 'linear';
end
% Find the extreme maxim values
% and the corresponding indexes
%----------------------------------------------------
extrMaxValue = y(find(diff(sign(diff(y)))==-2)+1);
extrMaxIndex = find(diff(sign(diff(y)))==-2)+1;
% Find the extreme minim values
% and the corresponding indexes
%----------------------------------------------------
extrMinValue = y(find(diff(sign(diff(y)))==+2)+1);
extrMinIndex = find(diff(sign(diff(y)))==+2)+1;
up = extrMaxValue;
up_x = x(extrMaxIndex);
down = extrMinValue;
down_x = x(extrMinIndex);
% Interpolation of the upper/down envelope data
%----------------------------------------------------
up = interp1(up_x,up,x,interpMethod);
down = interp1(down_x,down,x,interpMethod);
效果:
源程序下载:http://www.pudn.com/downloads495/sourcecode/windows/other/detail2062657.html
csdn上传资源标签里竟然没有Matlab,桑不起啊。。。。
零阶一类贝塞尔函数曲线及其包络线 matlab源程序相关推荐
- matlab 第一类修正贝塞尔函数,零阶贝塞尔函数 在MATLAB中怎样画出零阶修正贝塞尔函数...
第二类修正贝塞尔函数的零阶和一阶,分别怎样用mat回忆终究只是回忆,它只代表一段过去,一段历史,回忆再美也只是曾经,告别过去,期待未来. 书上说 (0和1都是下标) K0(z), the zeroth ...
- 自己实现MATLAB支持向量化输入输出的零阶修正贝塞尔函数
function [ y ] = Bessel0( x ) % 此函数用于计算0阶修正贝塞尔函数值 % 支持向量化输入输出 % 作用相当于MATLAB自带的库函数besseli(0,x) theta ...
- MatLab画贝塞尔函数曲线
原文地址:http://blog.sina.com.cn/s/blog_a229a65b01012jc5.html 第一类贝塞尔函数 在MatLab中用besselj(NU,Z)来表示: 用MatLa ...
- 0阶贝塞尔函数_从零开始的matlab学习笔记——(7)导函数
上节我们简单记录了一下怎么求极限,这节就来记录一下怎么求导函数. 显函数求导: 废话不多说,我们还是先看看代码 syms x f; f=sqrt(1+exp(x)); ans1=diff(f,x) p ...
- matlab 非线性仿真,科学网—含高阶色散和高阶非线性项的非线性耦合仿真Matlab源程序 - 王又法的博文...
如下Matlab 程序计算 论文Youfa Wang and Wenfeng Wang, Study of ultrafast pulse coupling dynamics considering ...
- matlab 0阶贝塞尔函数,零阶贝塞尔函数
第二类修正贝塞尔函数的零阶和一阶,分别怎样用mat书上说 (0和1都是下标) K0(z), the zeroth order modified Bessel funct首先双击matlab软件图标,打 ...
- matlab 0阶贝塞尔函数,零阶贝塞尔函数 分享助,第一类零阶贝塞尔
分享助,第一类零阶贝塞尔函数积分和sinc函数之间有也就是下面的这个等式怎么得出来的. 第二类修正贝塞尔函数的零阶和一阶,分别怎样用mat 书上说 (0和1都是下标) K0(z), the zerot ...
- C语言求解贝塞尔函数方程,任意阶复宗量贝塞尔函数的数值计算.pdf
任意阶复宗量贝塞尔函数的数值计算,零阶贝塞尔函数,一阶贝塞尔函数,0阶贝塞尔函数,第一类零阶贝塞尔函数,零阶修正贝塞尔函数,0阶修正贝塞尔函数,n阶贝塞尔函数近似,贝塞尔函数,matlab贝塞尔函数 ...
- python贝塞尔插值公式_贝塞尔函数、Jinc函数与图像插值算法
Jinc函数插值是一种传统的图像插值算法,可以用来进行图像缩放(即Jinc Resize).在传统的非机器学习图像超分辨率/放大算法中,Jinc Resize是效果较为出众的算法.虽然在传统神经网络和 ...
最新文章
- 【golang程序包推荐分享】go-ini、viper、godoc
- python环境配置是什么意思_python环境搭建
- 关于在linux下出现stdio.h文件不存在等gcc标准库不能找到的解决的方法
- camel_Apache Camel 2.14中的更多指标
- web基础html元素制作web
- 【数据竞赛】百赛百试,十试九灵的特征筛选策略-Pearson Correlation
- QGraphicsWidget收不到鼠标、键盘消息解决
- Sentinel(二十)之Envoy RLS Token Server
- 深度学习(05)--典型CNN结构(VGG13,16,19)
- 红包随机分配问题php,自己写仿微信红包的随机分配,为什么平均数差距这么大?...
- matlab深度学习基本操作,深度学习讲解,源代码分享,用一用改一改随便发篇EI
- 软件工程实践——软件评测
- CentOS | 资产管理软件GLPI+OCS安装与配置
- html5好看的颜色代码,css好看的颜色配色.html
- 怎么在我们的App中集成条码扫描功能?
- 【C语言开源库】C语言必备实用第三方库Melon(包括数据结构算法)
- spring-cloud-starter-bus-kafka利用kafka消息总线实现动态刷新配置
- Python 招聘信息爬取及可视化
- 基于Python+kociemba+opencv的3阶魔方自动还原
- iOS快捷指令:执行Python脚本(利用iSH Shell)
热门文章
- 如何写个优秀的Github项目Readme文档?经典模版拿去不谢~
- 借助github搭建自己的CDN服务
- 知识是不会让人自满的。 少年,只有名和利才会。
- 考试系统mysql数据库设计_在线考试系统数据库设计(表)
- 编程之类的文案_有什么有逼格的四字文案?
- java抽象类案例(有注释)
- 00后步入社会,有哪些适合你的创业项目?
- 慕课软件质量保证与测试(第三章.课后作业)
- 51单片机 AT89C52 的中断/定时器/计时器 笔记
- 黄河水稻山东智慧 国稻种芯·中国水稻节:济南泉城米袋子