i=0:0.1:45 ; i0=i*pi/180; r1=70*pi/180; m1=195; K=1500; L=3025; qwz=acot(cot(i0)-K/L); %期望值 P1=sin(r1+i0); Q1=(K/m1)*(2*cos(r1)-cos(r1+i0))-cos(2*r1); R1=sqrt((K/m1)^2+1-2*(K/m1)*cos(r1+i0)); yhq=real(r1-asin(P1./R1)-acos(Q1./R1)); %实际值 plot(i,qwz,'b-'); grid on; hold on; %画在同一张图上 plot(i,yhq,'r+'); hold on; xlabel('外轮转角'); ylabel('内轮转角'); legend('理论转角','实际转角'); title('理论转角VS实际转角');将其纵坐标单位弧度改为度

您可以在绘制图形前将纵坐标的单位从弧度转换为度，然后再使用`plot`函数绘制图形。以下是修改后的代码：

i=0:0.1:45;
i0=i*pi/180;
r1=70*pi/180;
m1=195;
K=1500;
L=3025;
qwz=acot(cot(i0)-K/L);    %期望值
P1=sin(r1+i0);
Q1=(K/m1)*(2*cos(r1)-cos(r1+i0))-cos(2*r1);
R1=sqrt((K/m1)^2+1-2*(K/m1)*cos(r1+i0));
yhq=real(r1-asin(P1./R1)-acos(Q1./R1));     %实际值

% 将纵坐标的单位从弧度转换为度
qwz = qwz * 180 / pi;
yhq = yhq * 180 / pi;

plot(i,qwz,'b-');
grid on; hold on;
plot(i,yhq,'r+'); hold on;
xlabel('外轮转角'); 
ylabel('内轮转角（度）'); 
legend('理论转角','实际转角'); 
title('理论转角VS实际转角');

在这个示例中，我们将纵坐标的单位从弧度转换为度，并将纵坐标的标签修改为“内轮转角（度）”。然后，我们使用`plot`函数绘制修改后的数据。