谁能提供复变量伽马函数的一些函数值?谁能利用计算软件提供复变量伽马函数Γ(x+yi)的一些函数值?范围是:0<y≤1≤x<2,(1)最低精度要求:x=1、1.1、1.2、1.3、1.4、1.5、1.6、1.7、1.8、1.9
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/02 00:30:55
![谁能提供复变量伽马函数的一些函数值?谁能利用计算软件提供复变量伽马函数Γ(x+yi)的一些函数值?范围是:0<y≤1≤x<2,(1)最低精度要求:x=1、1.1、1.2、1.3、1.4、1.5、1.6、1.7、1.8、1.9](/uploads/image/z/3788087-23-7.jpg?t=%E8%B0%81%E8%83%BD%E6%8F%90%E4%BE%9B%E5%A4%8D%E5%8F%98%E9%87%8F%E4%BC%BD%E9%A9%AC%E5%87%BD%E6%95%B0%E7%9A%84%E4%B8%80%E4%BA%9B%E5%87%BD%E6%95%B0%E5%80%BC%3F%E8%B0%81%E8%83%BD%E5%88%A9%E7%94%A8%E8%AE%A1%E7%AE%97%E8%BD%AF%E4%BB%B6%E6%8F%90%E4%BE%9B%E5%A4%8D%E5%8F%98%E9%87%8F%E4%BC%BD%E9%A9%AC%E5%87%BD%E6%95%B0%CE%93%EF%BC%88x%2Byi%EF%BC%89%E7%9A%84%E4%B8%80%E4%BA%9B%E5%87%BD%E6%95%B0%E5%80%BC%3F%E8%8C%83%E5%9B%B4%E6%98%AF%EF%BC%9A0%EF%BC%9Cy%E2%89%A41%E2%89%A4x%EF%BC%9C2%2C%EF%BC%881%EF%BC%89%E6%9C%80%E4%BD%8E%E7%B2%BE%E5%BA%A6%E8%A6%81%E6%B1%82%EF%BC%9Ax%3D1%E3%80%811.1%E3%80%811.2%E3%80%811.3%E3%80%811.4%E3%80%811.5%E3%80%811.6%E3%80%811.7%E3%80%811.8%E3%80%811.9)
谁能提供复变量伽马函数的一些函数值?谁能利用计算软件提供复变量伽马函数Γ(x+yi)的一些函数值?范围是:0<y≤1≤x<2,(1)最低精度要求:x=1、1.1、1.2、1.3、1.4、1.5、1.6、1.7、1.8、1.9
谁能提供复变量伽马函数的一些函数值?
谁能利用计算软件提供复变量伽马函数Γ(x+yi)的一些函数值?范围是:0<y≤1≤x<2,
(1)最低精度要求:x=1、1.1、1.2、1.3、1.4、1.5、1.6、1.7、1.8、1.9;y=0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8、0.9、1
(2)提供高精度函数值者,本人再加分.
谁能提供复变量伽马函数的一些函数值?谁能利用计算软件提供复变量伽马函数Γ(x+yi)的一些函数值?范围是:0<y≤1≤x<2,(1)最低精度要求:x=1、1.1、1.2、1.3、1.4、1.5、1.6、1.7、1.8、1.9
matlab中定义了gamma(x)函数,但x为实数;于是我用了下面的程序计算复数形式的gamma函数:
clear;clc
x=[1.1:0.1:2];y=0.1:0.1:1;%起点终点和步长可以自己改
[X,Y]=meshgrid(x,y);
syms t;
f=t.^(X+Y*i-1)*exp(-t) %定义积分函数
ga=int(f,t,0,inf) %积分
vpa(f,6) %查看结果
结果在x=1.1时会出现错误,当x=1.1时,积分函数为:
t^(i/10 + 1)/exp(t)
t^(i/5 + 1)/exp(t)
t^((3*i)/10 + 1)/exp(t)
t^((2*i)/5 + 1)/exp(t)
t^(i/2 + 1)/exp(t)
t^((3*i)/5 + 1)/exp(t)
t^((7*i)/10 + 1)/exp(t)
t^((4*i)/5 + 1)/exp(t)
t^((9*i)/10 + 1)/exp(t)
t^(i + 1)/exp(t)
即其中的实部分本应为0.1,这里却为1;这是产生错误的原因,从而用下面的试算x=1.1的情况:
b=0.1+Y(:,1)*i
b =
0.1000 + 0.1000i
0.1000 + 0.2000i
0.1000 + 0.3000i
0.1000 + 0.4000i
0.1000 + 0.5000i
0.1000 + 0.6000i
0.1000 + 0.7000i
0.1000 + 0.8000i
0.1000 + 0.9000i
0.1000 + 1.0000i
>> c=t.^b*exp(-t)
c =
t^(i/10 + 1/10)/exp(t) %这里的实部正确,为0.1=1/10;
t^(i/5 + 1/10)/exp(t)
t^((3*i)/10 + 1/10)/exp(t)
t^((2*i)/5 + 1/10)/exp(t)
t^(i/2 + 1/10)/exp(t)
t^((3*i)/5 + 1/10)/exp(t)
t^((7*i)/10 + 1/10)/exp(t)
t^((4*i)/5 + 1/10)/exp(t)
t^((9*i)/10 + 1/10)/exp(t)
t^(i + 1/10)/exp(t)
>> int(c,t,0,inf) %对x=1.1时积分
ans =
gamma(i/10 + 11/10)
gamma(i/5 + 11/10)
gamma((3*i)/10 + 11/10)
gamma((2*i)/5 + 11/10)
gamma(i/2 + 11/10)
gamma((3*i)/5 + 11/10)
gamma((7*i)/10 + 11/10)
gamma((4*i)/5 + 11/10)
gamma((9*i)/10 + 11/10)
gamma(i + 11/10)
>> vpa(ans,6) %查看积分结果
ans =
0.943739-0.0397233*i
0.921608-0.0760291*i
0.886905-0.106088*i
0.842395-0.128064*i
0.791143-0.141237*i
0.736052-0.145883*i
0.679561-0.143*i
0.623529-0.13399*i
0.569237-0.1204*i
0.517481-0.103723*i
这个结果是正确的了.这个结果与最上方的程序之间产生矛盾,为什么会这样我也没有搞懂.