GNUPLOT - 在一个范围内绘制一个函数
GNUPLOT - plotting a function in a range
我对我的情节有疑问。我想在一个范围内绘制一些函数，但我不知道如何定义 x 轴，因为有不同的变量。
这是我的 gnuplot 绘图命令：
reset
set grid
#set autoscale xy
set yrange [0:150]
set xrange [0:100]
set ytics 10 nomirror
sig = abs(1 - (x1 * x2 / xh**2))
print sig
set dummy s
p = 2
fn = 50.0
f1 = 10.0
f2 = 25.0
f3 = 35.0
f5 = 60.0
f6 = 75.0
f7 = 100.0
un = 400.0
u1 = f1/fn * un
u2 = f2/fn * un
u3 = f3/fn * un
u5 = un
u6 = un
u7 = un
ns1 = f1 / p
ns2 = f2 / p
ns3 = f3 / p
ns4 = fn / p
ns5 = f5 / p
ns6 = f6 / p
ns7 = f7 / p
skn = r2 / (sig * x2)
sk1 = skn * fn / f1
sk2 = skn * fn / f2
sk3 = skn * fn / f3
sk5 = skn * fn / f5
sk6 = skn * fn / f6
sk7 = skn * fn / f7
s1(s) = s * sk1
s2(s) = s * sk2
s3(s) = s * sk3
s4(s) = s * skn
s5(s) = s * sk5
s6(s) = s * sk6
s7(s) = s * sk7
n1(s) = ns1 * (1-s1(s))
n2(s) = ns2 * (1-s2(s))
n3(s) = ns3 * (1-s3(s))
n4(s) = ns4 * (1-s4(s))
n5(s) = ns5 * (1-s5(s))
n6(s) = ns6 * (1-s6(s))
n7(s) = ns7 * (1-s7(s))
mkn = 3 * un**2 * (1 - sig) / (4 * pi * fn/p * sig * x1)
mk1 = mkn
mk2 = mkn
mk3 = mkn
mk5 = mkn * (u5 / un * fn / f5)**2
mk6 = mkn * (u6 / un * fn / f6)**2
mk7 = mkn * (u7 / un * fn / f7)**2
print("MkN")
print mkn
print mk1
print mk2
print mk3
print mk5
print mk6
print mk7
mi1(s) = mkn * 2 / (s1(s) / sk1 + sk1 / s1(s))
mi2(s) = mkn * 2 / (s2(s) / sk2 + sk2 / s2(s))
mi3(s) = mkn * 2 / (s3(s) / sk3 + sk3 / s3(s))
mi4(s) = mkn * 2 / (s4(s) / skn + skn / s4(s))
mi5(s) = mkn * 2 / (s5(s) / sk5 + sk5 / s5(s))
mi6(s) = mkn * 2 / (s6(s) / sk6 + sk6 / s6(s))
mi7(s) = mkn * 2 / (s7(s) / sk7 + sk7 / s7(s))
set samples 10E3
plot mi1(s) / (0<=s) / (s<=1.5) w l lw 1.5 lc "blue" title " M_{i10Hz}"
它看起来不像它应该的样子。
我们有一个适用于 mathlab 的 equal plot 命令。
这是：
clear all;
clc;
R2=1;
X1sigma=5;
X2sigma=3;
Xh=45;
X1=X1sigma+Xh;
X2=X2sigma+Xh;
sig=abs(1- (X1*X2)/(abs(Xh^2)));
%frequenzen
f1n=50;
f1 = 10;
f2 = 25;
f3 = 35;
f4 = 50;
f5 = 60;
f6 = 75;
f7 = 100;
%Synchrondrehzahlen
ns=25;
ns1 = 5;
ns2 = 12.5;
ns3 = 17.5;
ns4 = 25;
ns5 = 30;
ns6 = 37.5;
ns7 = 50;
%Kippschlupf
Skn=R2 / (sig * X2);
Sk1= Skn * f1n / f1;
Sk2= Skn * f1n / f2;
Sk3= Skn * f1n / f3;
Sk4= Skn * f1n / f4;
Sk5= Skn * f1n / f5;
Sk6= Skn * f1n / f6;
Sk7= Skn * f1n / f7;
%Schlupf
S1= 0:0.01:1.5 *Sk1;
S2= 0:0.01:1.5 *Sk2;
S3= 0:0.01:1.5 *Sk3;
S4= 0:0.01:1.5 *Sk4;
S5= 0:0.01:1.5 *Sk5;
S6= 0:0.01:1.5 *Sk6;
S7= 0:0.01:1.5 *Sk7;
%Spannungen
U1N=400;
U1=U1N/5;
U2=U1N/2;
U3=U1N/1.43;
U4=U1N;
U5=U1N;
U6=U1N;
U7=U1N;
%Drehzahlen
n1= ns1 .* ( 1 - S1 );
n2= ns2 .* ( 1 - S2 );
n3= ns3 .* ( 1 - S3 );
n4= ns4 .* ( 1 - S4 );
n5= ns5 .* ( 1 - S5 );
n6= ns6 .* ( 1 - S6 );
n7= ns7 .* ( 1 - S7 );
%Kippmomente
Mkn= (3*(U1N^2)*(1-sig)) / (4 * pi * ns * sig * X1 );
Mk1=Mkn*(((U1/U1N)*(f1n/f1))^2);
Mk2=Mkn*(((U2/U1N)*(f1n/f2))^2);
Mk3=Mkn*(((U3/U1N)*(f1n/f3))^2);
Mk4=Mkn*(((U4/U1N)*(f1n/f4))^2);
Mk5=Mkn*(((U5/U1N)*(f1n/f5))^2);
Mk6=Mkn*(((U6/U1N)*(f1n/f6))^2);
Mk7=Mkn*(((U7/U1N)*(f1n/f7))^2);
%Drehmomente
Mn=33.705;
Mi1=(Mk1 .* 2.0 ) ./ ( ( S1 ./ Sk1 ) + ( Sk1 ./ S1 ) );
Mi2=(Mk2 .* 2.0 ) ./ ( ( S2 ./ Sk2 ) + ( Sk2 ./ S2 ) );
Mi3=(Mk3 .* 2.0 ) ./ ( ( S3 ./ Sk3 ) + ( Sk3 ./ S3 ) );
Mi4=(Mk4 .* 2.0 ) ./ ( ( S4 ./ Sk4 ) + ( Sk4 ./ S4 ) );
Mi5=(Mk5 .* 2.0 ) ./ ( ( S5 ./ Sk5 ) + ( Sk5 ./ S5 ) );
Mi6=(Mk6 .* 2.0 ) ./ ( ( S6 ./ Sk6 ) + ( Sk6 ./ S6 ) );
Mi7=(Mk7 .* 2.0 ) ./ ( ( S7 ./ Sk7 ) + ( Sk7 ./ S7 ) );
%plots
plot(2*n1, Mi1,'LineWidth', 2);
hold on;
plot(2*n2, Mi2,'LineWidth', 2);
hold on;
plot(2*n3, Mi3,'LineWidth', 2);
hold on;
plot(2*n4, Mi4,'LineWidth', 2);
hold on;
plot(2*n5, Mi5,'LineWidth', 2);
hold on;
plot(2*n6, Mi6,'LineWidth', 2);
hold on;
plot(2*n7, Mi7,'LineWidth', 2);
hold on;
axis([0 100 0 150]);
xlabel("n/1/s");
ylabel("M/Nm");
title("")
%Legende
hl = legend("10Hz", "25Hz", "35Hz", "50Hz", "60Hz", "75HZ", "100Hz" );
set(hl,'Interpreter','latex','Location', 'Northeast','FontSize',12);
grid on;
grid minor;
它应该是这样的。
也许有人可以帮助我。谢谢
如评论中所述，您的 gnuplot 脚本在多个定义中与您的 Matlab 脚本不匹配。
特别是在你的情况下，我强烈建议定义数组和函数。这将避免复制粘贴错误，并使代码更短、更清晰。
我希望你使用的是 gnuplot >=5.2，因为下面的代码使用了数组。
我从你的 Matlab 脚本开始，稍微清理了一下。在 gnuplot 的一行末尾使用分号 ; 不是必需的，但它也没有坏处。 Matlab 中的功能 ^ 是 gnuplot 中的 **。顺便说一句，在您的 Matlab 代码中，我看不到变量 Mn=33.705 的使用位置。
此外，如评论中所述，在 gnuplot 中使用除法时要注意，因为如果没有一个数字是浮点数，gnuplot 将使用整数除法。这是 gnuplot 获得意想不到的结果的最卑鄙的陷阱之一。
对于颜色，您还可以使用 0xRRGGBB.
格式的 RGB hexadecimal-numbers
代码：
### plotting using arrays and functions
reset session
R2=1;
X1sigma=5.0;
X2sigma=3.0;
Xh=45.0;
X1=X1sigma+Xh;
X2=X2sigma+Xh;
sig=abs(1- (X1*X2)/(abs(Xh**2)));
#%frequenzen
f1n=50.0;
array fArr[7] = [10., 25., 35., 50., 60., 75., 100.]
#%Synchrondrehzahlen
ns=25;
array nsArr[7] = [5., 12.5, 17.5, 25., 30., 37.5, 50.]
#%Kippschlupf
Skn=R2 / (sig * X2);
Sk(i) = Skn * f1n / fArr[i]
#%Schlupf
S(s,i) = s * Sk(i)
#%Spannungen
U1N=400;
array UdivArr[7] = [5.0, 2.0, 1.43, 1.0, 1.0, 1.0, 1.0]
U(i) = U1N / UdivArr[i]
#%Drehzahlen
n(s,i) = nsArr[i] * (1-S(s,i))
#%Kippmomente
Mkn= (3*(U1N**2)*(1-sig)) / (4 * pi * ns * sig * X1 );
Mk(i) = Mkn*(((U(i)/U1N)*(f1n/fArr[i]))**2);
#%Drehmomente
Mn=33.705;
Mi(s,i) = (Mk(i)*2.0 ) / ( ( S(s,i)/ Sk(i) ) + ( Sk(i)/S(s,i) ) )
array myColor[7] = ['blue', 'red', 'orange', 'dark-violet', 'green', 'light-blue', 'brown' ]
set title "M = f(n)"
set key reverse
set xlabel "n / 1/s"
set ylabel "M / Nm"
set grid xtics, ytics
plot for [i=1:7] [s=0:1.5:0.01] '+' u (2*n(s,i)):(Mi(s,i)) w l lw 2 lc rgb myColor[i] title sprintf("%g Hz",fArr[i])
### end of code
补充：（plot命令的一些解释）
plot for [i=1:7] 将 7 个图循环成一个索引为 i[s=0:1.5:0.01] '+' 此“生成”s 的值作为参数，从 0 到 1.5，步长为 0.01u (2*n(s,i)) 使用 x-value from 函数 s 和 i 作为参数:(Mi(s,i)) y-value 来自以 s 和 i 作为参数的函数 w l lw 2 线宽 2lc rgb myColor[i] 从数组 myColortitle sprintf("%g Hz",fArr[i]) 标题使用数组 fArr
检查关键字的 gnuplot 帮助，例如help plot for、help sample、help sprintf 等
结果：
我对我的情节有疑问。我想在一个范围内绘制一些函数，但我不知道如何定义 x 轴，因为有不同的变量。 这是我的 gnuplot 绘图命令：
reset
set grid
#set autoscale xy
set yrange [0:150]
set xrange [0:100]
set ytics 10 nomirror
sig = abs(1 - (x1 * x2 / xh**2))
print sig
set dummy s
p = 2
fn = 50.0
f1 = 10.0
f2 = 25.0
f3 = 35.0
f5 = 60.0
f6 = 75.0
f7 = 100.0
un = 400.0
u1 = f1/fn * un
u2 = f2/fn * un
u3 = f3/fn * un
u5 = un
u6 = un
u7 = un
ns1 = f1 / p
ns2 = f2 / p
ns3 = f3 / p
ns4 = fn / p
ns5 = f5 / p
ns6 = f6 / p
ns7 = f7 / p
skn = r2 / (sig * x2)
sk1 = skn * fn / f1
sk2 = skn * fn / f2
sk3 = skn * fn / f3
sk5 = skn * fn / f5
sk6 = skn * fn / f6
sk7 = skn * fn / f7
s1(s) = s * sk1
s2(s) = s * sk2
s3(s) = s * sk3
s4(s) = s * skn
s5(s) = s * sk5
s6(s) = s * sk6
s7(s) = s * sk7
n1(s) = ns1 * (1-s1(s))
n2(s) = ns2 * (1-s2(s))
n3(s) = ns3 * (1-s3(s))
n4(s) = ns4 * (1-s4(s))
n5(s) = ns5 * (1-s5(s))
n6(s) = ns6 * (1-s6(s))
n7(s) = ns7 * (1-s7(s))
mkn = 3 * un**2 * (1 - sig) / (4 * pi * fn/p * sig * x1)
mk1 = mkn
mk2 = mkn
mk3 = mkn
mk5 = mkn * (u5 / un * fn / f5)**2
mk6 = mkn * (u6 / un * fn / f6)**2
mk7 = mkn * (u7 / un * fn / f7)**2
print("MkN")
print mkn
print mk1
print mk2
print mk3
print mk5
print mk6
print mk7
mi1(s) = mkn * 2 / (s1(s) / sk1 + sk1 / s1(s))
mi2(s) = mkn * 2 / (s2(s) / sk2 + sk2 / s2(s))
mi3(s) = mkn * 2 / (s3(s) / sk3 + sk3 / s3(s))
mi4(s) = mkn * 2 / (s4(s) / skn + skn / s4(s))
mi5(s) = mkn * 2 / (s5(s) / sk5 + sk5 / s5(s))
mi6(s) = mkn * 2 / (s6(s) / sk6 + sk6 / s6(s))
mi7(s) = mkn * 2 / (s7(s) / sk7 + sk7 / s7(s))
set samples 10E3
plot mi1(s) / (0<=s) / (s<=1.5) w l lw 1.5 lc "blue" title " M_{i10Hz}"
它看起来不像它应该的样子。 我们有一个适用于 mathlab 的 equal plot 命令。 这是：
clear all;
clc;
R2=1;
X1sigma=5;
X2sigma=3;
Xh=45;
X1=X1sigma+Xh;
X2=X2sigma+Xh;
sig=abs(1- (X1*X2)/(abs(Xh^2)));
%frequenzen
f1n=50;
f1 = 10;
f2 = 25;
f3 = 35;
f4 = 50;
f5 = 60;
f6 = 75;
f7 = 100;
%Synchrondrehzahlen
ns=25;
ns1 = 5;
ns2 = 12.5;
ns3 = 17.5;
ns4 = 25;
ns5 = 30;
ns6 = 37.5;
ns7 = 50;
%Kippschlupf
Skn=R2 / (sig * X2);
Sk1= Skn * f1n / f1;
Sk2= Skn * f1n / f2;
Sk3= Skn * f1n / f3;
Sk4= Skn * f1n / f4;
Sk5= Skn * f1n / f5;
Sk6= Skn * f1n / f6;
Sk7= Skn * f1n / f7;
%Schlupf
S1= 0:0.01:1.5 *Sk1;
S2= 0:0.01:1.5 *Sk2;
S3= 0:0.01:1.5 *Sk3;
S4= 0:0.01:1.5 *Sk4;
S5= 0:0.01:1.5 *Sk5;
S6= 0:0.01:1.5 *Sk6;
S7= 0:0.01:1.5 *Sk7;
%Spannungen
U1N=400;
U1=U1N/5;
U2=U1N/2;
U3=U1N/1.43;
U4=U1N;
U5=U1N;
U6=U1N;
U7=U1N;
%Drehzahlen
n1= ns1 .* ( 1 - S1 );
n2= ns2 .* ( 1 - S2 );
n3= ns3 .* ( 1 - S3 );
n4= ns4 .* ( 1 - S4 );
n5= ns5 .* ( 1 - S5 );
n6= ns6 .* ( 1 - S6 );
n7= ns7 .* ( 1 - S7 );
%Kippmomente
Mkn= (3*(U1N^2)*(1-sig)) / (4 * pi * ns * sig * X1 );
Mk1=Mkn*(((U1/U1N)*(f1n/f1))^2);
Mk2=Mkn*(((U2/U1N)*(f1n/f2))^2);
Mk3=Mkn*(((U3/U1N)*(f1n/f3))^2);
Mk4=Mkn*(((U4/U1N)*(f1n/f4))^2);
Mk5=Mkn*(((U5/U1N)*(f1n/f5))^2);
Mk6=Mkn*(((U6/U1N)*(f1n/f6))^2);
Mk7=Mkn*(((U7/U1N)*(f1n/f7))^2);
%Drehmomente
Mn=33.705;
Mi1=(Mk1 .* 2.0 ) ./ ( ( S1 ./ Sk1 ) + ( Sk1 ./ S1 ) );
Mi2=(Mk2 .* 2.0 ) ./ ( ( S2 ./ Sk2 ) + ( Sk2 ./ S2 ) );
Mi3=(Mk3 .* 2.0 ) ./ ( ( S3 ./ Sk3 ) + ( Sk3 ./ S3 ) );
Mi4=(Mk4 .* 2.0 ) ./ ( ( S4 ./ Sk4 ) + ( Sk4 ./ S4 ) );
Mi5=(Mk5 .* 2.0 ) ./ ( ( S5 ./ Sk5 ) + ( Sk5 ./ S5 ) );
Mi6=(Mk6 .* 2.0 ) ./ ( ( S6 ./ Sk6 ) + ( Sk6 ./ S6 ) );
Mi7=(Mk7 .* 2.0 ) ./ ( ( S7 ./ Sk7 ) + ( Sk7 ./ S7 ) );
%plots
plot(2*n1, Mi1,'LineWidth', 2);
hold on;
plot(2*n2, Mi2,'LineWidth', 2);
hold on;
plot(2*n3, Mi3,'LineWidth', 2);
hold on;
plot(2*n4, Mi4,'LineWidth', 2);
hold on;
plot(2*n5, Mi5,'LineWidth', 2);
hold on;
plot(2*n6, Mi6,'LineWidth', 2);
hold on;
plot(2*n7, Mi7,'LineWidth', 2);
hold on;
axis([0 100 0 150]);
xlabel("n/1/s");
ylabel("M/Nm");
title("")
%Legende
hl = legend("10Hz", "25Hz", "35Hz", "50Hz", "60Hz", "75HZ", "100Hz" );
set(hl,'Interpreter','latex','Location', 'Northeast','FontSize',12);
grid on;
grid minor;
它应该是这样的。
也许有人可以帮助我。谢谢
如评论中所述，您的 gnuplot 脚本在多个定义中与您的 Matlab 脚本不匹配。
特别是在你的情况下，我强烈建议定义数组和函数。这将避免复制粘贴错误，并使代码更短、更清晰。 我希望你使用的是 gnuplot >=5.2，因为下面的代码使用了数组。
我从你的 Matlab 脚本开始，稍微清理了一下。在 gnuplot 的一行末尾使用分号 ; 不是必需的，但它也没有坏处。 Matlab 中的功能 ^ 是 gnuplot 中的 **。顺便说一句，在您的 Matlab 代码中，我看不到变量 Mn=33.705 的使用位置。
此外，如评论中所述，在 gnuplot 中使用除法时要注意，因为如果没有一个数字是浮点数，gnuplot 将使用整数除法。这是 gnuplot 获得意想不到的结果的最卑鄙的陷阱之一。
对于颜色，您还可以使用 0xRRGGBB.
代码：
### plotting using arrays and functions
reset session
R2=1;
X1sigma=5.0;
X2sigma=3.0;
Xh=45.0;
X1=X1sigma+Xh;
X2=X2sigma+Xh;
sig=abs(1- (X1*X2)/(abs(Xh**2)));
#%frequenzen
f1n=50.0;
array fArr[7] = [10., 25., 35., 50., 60., 75., 100.]
#%Synchrondrehzahlen
ns=25;
array nsArr[7] = [5., 12.5, 17.5, 25., 30., 37.5, 50.]
#%Kippschlupf
Skn=R2 / (sig * X2);
Sk(i) = Skn * f1n / fArr[i]
#%Schlupf
S(s,i) = s * Sk(i)
#%Spannungen
U1N=400;
array UdivArr[7] = [5.0, 2.0, 1.43, 1.0, 1.0, 1.0, 1.0]
U(i) = U1N / UdivArr[i]
#%Drehzahlen
n(s,i) = nsArr[i] * (1-S(s,i))
#%Kippmomente
Mkn= (3*(U1N**2)*(1-sig)) / (4 * pi * ns * sig * X1 );
Mk(i) = Mkn*(((U(i)/U1N)*(f1n/fArr[i]))**2);
#%Drehmomente
Mn=33.705;
Mi(s,i) = (Mk(i)*2.0 ) / ( ( S(s,i)/ Sk(i) ) + ( Sk(i)/S(s,i) ) )
array myColor[7] = ['blue', 'red', 'orange', 'dark-violet', 'green', 'light-blue', 'brown' ]
set title "M = f(n)"
set key reverse
set xlabel "n / 1/s"
set ylabel "M / Nm"
set grid xtics, ytics
plot for [i=1:7] [s=0:1.5:0.01] '+' u (2*n(s,i)):(Mi(s,i)) w l lw 2 lc rgb myColor[i] title sprintf("%g Hz",fArr[i])
### end of code
补充：（plot命令的一些解释）
plot for [i=1:7]将 7 个图循环成一个索引为 i 的图
[s=0:1.5:0.01] '+'此“生成”s的值作为参数，从 0 到 1.5，步长为 0.01u (2*n(s,i))使用 x-value from 函数s和i作为参数:(Mi(s,i))y-value 来自以s和i作为参数的函数w l lw 2线宽 2lc rgb myColor[i]从数组 myColor 中取出线条颜色
title sprintf("%g Hz",fArr[i])标题使用数组 fArr
检查关键字的 gnuplot 帮助，例如help plot for、help sample、help sprintf 等
结果：
