-
大小: 31KB文件类型: .m金币: 1下载: 0 次发布日期: 2021-05-15
- 语言: Matlab
- 标签: 数值积分 matlab Gauss-laguer
资源简介
高斯-拉盖尔积分(Gauss-laguerre积分)matlab程序
可直接在函数名前设置输出为 [x,w]= 其中x为积分点,w为权重
代码片段和文件信息
function laguerre_rule ( order a b filename )
%*****************************************************************************80
% 可直接在函数名前设置输出为 [xw]= 其中x为积分点,w为权重,此处为求积分,因此w之和可能不为1。
%% LAGUERRE_RULE generates a Gauss-Laguerre rule.
%
% Discussion:
%
% This program computes a standard or exponentially weighted
% Gauss-Laguerre quadrature rule and writes it to a file.
%
% The user specifies:
% * the ORDER (number of points) in the rule;
% * A the left endpoint;
% * B the scale factor in the exponential;
% * FILENAME the root name of the output files.
%
% Licensing:
%
% This code is distributed under the GNU LGPL license.
%
% Modified:
%
% 10 February 2014
%
% Author:
%
% John Burkardt
%
timestamp ( );
fprintf ( 1 ‘\n‘ );
fprintf ( 1 ‘LAGUERRE_RULE\n‘ );
fprintf ( 1 ‘ MATLAB version\n‘ );
fprintf ( 1 ‘\n‘ );
fprintf ( 1 ‘ Compute a Gauss-Laguerre rule for approximating\n‘ );
fprintf ( 1 ‘ Integral ( a <= x < +oo ) exp(-b*(x-a)) f(x) dx\n‘ );
fprintf ( 1 ‘ of order ORDER.\n‘ );
fprintf ( 1 ‘\n‘ );
fprintf ( 1 ‘ The user specifies ORDER A B and FILENAME.\n‘ );
fprintf ( 1 ‘\n‘ );
fprintf ( 1 ‘ ORDER is the number of points.\n‘ );
fprintf ( 1 ‘ A is the left endpoint (typically 0).\n‘ );
fprintf ( 1 ‘ B is the exponential scale factor (typically 1).\n‘ );
fprintf ( 1 ‘ FILENAME is used to generate 3 files:\n‘ );
fprintf ( 1 ‘ * filename_w.txt - the weight file\n‘ );
fprintf ( 1 ‘ * filename_x.txt - the abscissa file.\n‘ );
fprintf ( 1 ‘ * filename_r.txt - the region file.\n‘ );
%
% Initialize the parameters.
%
alpha = 0.0;
beta = 0.0;
%
% Get ORDER.
%
if ( nargin < 1 )
order = input ( ‘ Enter the rule order ORDER: ‘ );
elseif ( ischar ( order ) )
order = str2num ( order );
end
%
% Get A.
%
if ( nargin < 2 )
fprintf ( 1 ‘\n‘ );
fprintf ( 1 ‘ A is the left endpoint typically 0.\n‘ );
fprintf ( 1 ‘\n‘ );
a = input ( ‘ Enter the value of A.‘ );
elseif ( ischar ( a ) )
a = str2num ( a );
end
%
% Get B.
%
if ( nargin < 3 )
fprintf ( 1 ‘\n‘ );
fprintf ( 1 ‘ B is the exponential scale factor typically 1.\n‘ );
fprintf ( 1 ‘\n‘ );
b = input ( ‘ Enter the value of B.‘ );
elseif ( ischar ( b ) )
b = str2num ( b );
end
%
% Get the quadrature file root name:
%
if ( nargin < 4 )
fprintf ( 1 ‘\n‘ );
fprintf ( 1 ‘ FILENAME specifies the ‘‘root name‘‘ of a set of quadrature files).\n‘ );
filename = input ( ‘ Enter the value of FILENAME as a quoted string: ‘ );
end
%
% Input summary.
%
fprintf ( 1 ‘\n‘ );
fprintf ( 1 ‘ ORDER = %d\n‘ order );
fprintf ( 1 ‘ A = %f\n‘ a);
fprintf ( 1 ‘ B = %f\n‘ b );
fprintf ( 1 ‘ FILENAME = “%s“.\n‘ filename );
%
% Construct the rule.
%
kind = 5;
[ x w ] = cgqf ( order kind alpha beta
相关资源
- Matlab路面裂缝识别69319
- 高灵敏度GPS接收机MATLAB仿真,附捕获
- 基于MATLAB的质点弹道计算与外弹道优
- 阵列天线的matlab仿真
- MATLAB 经典程序源代码大全
- MATLAB小波软阈值去噪代码33473
- 天线阵的波束形成在MATLAB仿真程序及
- 非线性SVM算法-matlab实现
- 《MATLAB 智能算法超级学习手册》-程序
- 组合导航matlab程序
- 读取txt文件内容matlab代码实现
- Matlab实现基于相关的模板匹配程序
- matlab优化工具箱讲解
- 基于MATLAB的快速傅里叶变换
- 光纤传输中的分布傅立叶算法matlab实
- 基于matlab的图像处理源程序
- matlab 椭圆拟合程序
- 算术编码解码matlab源代码
- optical_flow 光流法 matlab 实现程序
- 引导图像滤波器 Matlab实现
- 分形几何中一些经典图形的Matlab画法
- OFDM系统MATLAB仿真代码
- SVM工具箱(matlab中运行)
- 图像小波变换MatLab源代码
- LU分解的MATLAB实现
- 冈萨雷斯数字图像处理matlab版(第三
- 替代数据法的matlab程序
- 用matlab实现的多站定位系统性能仿真
- 通过不同方法进行粗糙集属性约简m
- k近邻算法matlab实现
评论
共有 条评论