废话不多说,请看源程序:

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源程序相关推荐

  1. matlab 第一类修正贝塞尔函数,零阶贝塞尔函数 在MATLAB中怎样画出零阶修正贝塞尔函数...

    第二类修正贝塞尔函数的零阶和一阶,分别怎样用mat回忆终究只是回忆,它只代表一段过去,一段历史,回忆再美也只是曾经,告别过去,期待未来. 书上说 (0和1都是下标) K0(z), the zeroth ...

  2. 自己实现MATLAB支持向量化输入输出的零阶修正贝塞尔函数

    function [ y ] = Bessel0( x ) % 此函数用于计算0阶修正贝塞尔函数值 % 支持向量化输入输出 % 作用相当于MATLAB自带的库函数besseli(0,x) theta ...

  3. MatLab画贝塞尔函数曲线

    原文地址:http://blog.sina.com.cn/s/blog_a229a65b01012jc5.html 第一类贝塞尔函数 在MatLab中用besselj(NU,Z)来表示: 用MatLa ...

  4. 0阶贝塞尔函数_从零开始的matlab学习笔记——(7)导函数

    上节我们简单记录了一下怎么求极限,这节就来记录一下怎么求导函数. 显函数求导: 废话不多说,我们还是先看看代码 syms x f; f=sqrt(1+exp(x)); ans1=diff(f,x) p ...

  5. matlab 非线性仿真,科学网—含高阶色散和高阶非线性项的非线性耦合仿真Matlab源程序 - 王又法的博文...

    如下Matlab 程序计算 论文Youfa Wang and Wenfeng Wang,  Study of ultrafast pulse coupling dynamics considering ...

  6. matlab 0阶贝塞尔函数,零阶贝塞尔函数

    第二类修正贝塞尔函数的零阶和一阶,分别怎样用mat书上说 (0和1都是下标) K0(z), the zeroth order modified Bessel funct首先双击matlab软件图标,打 ...

  7. matlab 0阶贝塞尔函数,零阶贝塞尔函数 分享助,第一类零阶贝塞尔

    分享助,第一类零阶贝塞尔函数积分和sinc函数之间有也就是下面的这个等式怎么得出来的. 第二类修正贝塞尔函数的零阶和一阶,分别怎样用mat 书上说 (0和1都是下标) K0(z), the zerot ...

  8. C语言求解贝塞尔函数方程,任意阶复宗量贝塞尔函数的数值计算.pdf

    任意阶复宗量贝塞尔函数的数值计算,零阶贝塞尔函数,一阶贝塞尔函数,0阶贝塞尔函数,第一类零阶贝塞尔函数,零阶修正贝塞尔函数,0阶修正贝塞尔函数,n阶贝塞尔函数近似,贝塞尔函数,matlab贝塞尔函数 ...

  9. python贝塞尔插值公式_贝塞尔函数、Jinc函数与图像插值算法

    Jinc函数插值是一种传统的图像插值算法,可以用来进行图像缩放(即Jinc Resize).在传统的非机器学习图像超分辨率/放大算法中,Jinc Resize是效果较为出众的算法.虽然在传统神经网络和 ...

最新文章

  1. 【golang程序包推荐分享】go-ini、viper、godoc
  2. python环境配置是什么意思_python环境搭建
  3. 关于在linux下出现stdio.h文件不存在等gcc标准库不能找到的解决的方法
  4. camel_Apache Camel 2.14中的更多指标
  5. web基础html元素制作web
  6. 【数据竞赛】百赛百试,十试九灵的特征筛选策略-Pearson Correlation
  7. QGraphicsWidget收不到鼠标、键盘消息解决
  8. Sentinel(二十)之Envoy RLS Token Server
  9. 深度学习(05)--典型CNN结构(VGG13,16,19)
  10. 红包随机分配问题php,自己写仿微信红包的随机分配,为什么平均数差距这么大?...
  11. matlab深度学习基本操作,深度学习讲解,源代码分享,用一用改一改随便发篇EI
  12. 软件工程实践——软件评测
  13. CentOS | 资产管理软件GLPI+OCS安装与配置
  14. html5好看的颜色代码,css好看的颜色配色.html
  15. 怎么在我们的App中集成条码扫描功能?
  16. 【C语言开源库】C语言必备实用第三方库Melon(包括数据结构算法)
  17. spring-cloud-starter-bus-kafka利用kafka消息总线实现动态刷新配置
  18. Python 招聘信息爬取及可视化
  19. 基于Python+kociemba+opencv的3阶魔方自动还原
  20. iOS快捷指令:执行Python脚本(利用iSH Shell)

热门文章

  1. 如何写个优秀的Github项目Readme文档?经典模版拿去不谢~
  2. 借助github搭建自己的CDN服务
  3. 知识是不会让人自满的。 少年,只有名和利才会。
  4. 考试系统mysql数据库设计_在线考试系统数据库设计(表)
  5. 编程之类的文案_有什么有逼格的四字文案?
  6. java抽象类案例(有注释)
  7. 00后步入社会,有哪些适合你的创业项目?
  8. 慕课软件质量保证与测试(第三章.课后作业)
  9. 51单片机 AT89C52 的中断/定时器/计时器 笔记
  10. 黄河水稻山东智慧 国稻种芯·中国水稻节:济南泉城米袋子