求关于MATLAB复杂积分的曲线绘图.被积函数如图中所示,我想绘制一个用ezmesh函数的三维曲面图以及一个剖面的曲线图,以下是程序,但是运行不对,曲面图:c1 = 3.7148*10^(-16);c2 = 0.014388;dna = 0.24;Q1=@
来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/26 04:38:14
求关于MATLAB复杂积分的曲线绘图.被积函数如图中所示,我想绘制一个用ezmesh函数的三维曲面图以及一个剖面的曲线图,以下是程序,但是运行不对,曲面图:c1 = 3.7148*10^(-16);c2 = 0.014388;dna = 0.24;Q1=@
求关于MATLAB复杂积分的曲线绘图.
被积函数如图中所示,我想绘制一个用ezmesh函数的三维曲面图以及一个剖面的曲线图,以下是程序,但是运行不对,
曲面图:c1 = 3.7148*10^(-16);
c2 = 0.014388;
dna = 0.24;
Q1=@(tb)quad(@(x) tb*c1./x.^5./(exp(c2./(x.*tb))-1),8*10^(-6),12*10^(-6));
Q2=@(tb)quad(@(x) 4*c1./x.^5./(exp(c2./(x.*tb))-1)+2.8*10^(4)./(exp(1430tb)-1)-9.13*10^(4)./(exp(1800tb)-1),8*10^(-6),12*10^(-6));
K=@(x,tb) Q1(tb)/Q2(tb)/dna^2;
ezmesh(K,[8*10^(-6),12*10^(-6),200,500] )
zlabel('K')
title('K因子曲面图')
曲线图:c1 = 3.7148*10^(-16);
c2 = 0.014388;
dna = 0.24;
hold on
for tb = 200:50:500
F=@(x) tb*c1./x.^5./(exp(c2./(x.*tb))-1),8*10^(-6),12*10^(-6);
Q1=quad(F,8*10^(-6),12*10^(-6));
F2=@(x) 4*c1./x.^5./(exp(c2./(x.*tb))-1)+2.8*10^(4)./(exp(1430tb)-1)-9.13*10^(4)./(exp(1800tb)-1),8*10^(-6),12*10^(-6));
Q2 = quad(F2,8*10^(-6),12*10^(-6));
K=@(x) Q1(tb)/Q2(tb)/dna^2;
x =8*10^(-6):5*10^(-8):12*10^(-6);
plot(x,K(x),'.')
end
xlabel('x')
ylabel('K')
('K因子曲线图')
请大家指教,程序不知道怎么调试不成,
求关于MATLAB复杂积分的曲线绘图.被积函数如图中所示,我想绘制一个用ezmesh函数的三维曲面图以及一个剖面的曲线图,以下是程序,但是运行不对,曲面图:c1 = 3.7148*10^(-16);c2 = 0.014388;dna = 0.24;Q1=@
改了一下.你应再全面检查检查.
曲面图:
close all
clear,clc
c1 = 3.7148*10^(-16);
c2 = 0.014388;
dna = 0.24;
Q1 = @(tb) quad( @(x) tb*c1./x.^5./(exp(c2./(x.*tb))-1),8*10^(-6),12*10^(-6));
Q2 = @(tb) quad( @(x) 4*c1./x.^5./(exp(c2./(x.*tb))-1) + 2.8*10^4./(exp(1430.*tb)-1)-9.13*10^4./(exp(1800.*tb)-1),8*10^(-6),12*10^(-6));
M0 = @(x,tb) c1./x.^5./(exp(c2./(x.*tb))-1);
K = @(x,tb) Q1(tb)/Q2(tb)/dna^2./M0(x,tb);
ezmesh(K,[8*10^(-6), 12*10^(-6), 200, 500] )
zlabel('K')
title('K因子曲面图')
曲线图:
close all
clear,clc
c1 = 3.7148*10^(-16);
c2 = 0.014388;
dna = 0.24;
hold on
for tb = 200:50:500
F = @(x) tb*c1./x.^5./(exp(c2./(x*tb))-1);
Q1 = quad(F,8*10^(-6),12*10^(-6));
F2 = @(x) 4*c1./x.^5./(exp(c2./(x*tb))-1)+2.8*10^4./(exp(1430*tb)-1) - 9.13*10^4./(exp(1800*tb)-1);
Q2 = quad(F2,8*10^(-6),12*10^(-6));
M0 = @(x) c1./x.^5./(exp(c2./(x*tb))-1);
K = @(x) Q1/Q2/dna^2./M0(x);
x =8*10^(-6):5*10^(-8):12*10^(-6);
plot(x,K(x),'.')
end
xlabel('x')
ylabel('K')
title('K因子曲线图')