[程序] RCf300_.m
%f=300,RC电路李萨如图形随电阻的变化
clear,f=300;c=0.22*10^(-6);%清除变量,频率,电容
Um=10;ph0=0;%电圧峰值,初相
w=2*pi*f;Xc=1/(w*c)%圆频率,容抗
R=3000;z=sqrt(Xc^2+R^2),ph=atan(Xc/R)%电阻,阻抗,相位差
phd=ph*180/pi %相位差化为度数
Im=Um/z%电流峰值
t=linspace(0,0.005);%时间向量
u=Um*sin(w*t+ph0);%总电压
uR=Im*R*sin(w*t+ph0+ph);%电阻电压
uC=Im*Xc*sin(w*t+ph0-pi/2+ph);%电容电压
A=Im*R*sin(ph);%截距u=0,uR的值
B=Im*R;%uR的最大值
phexp=asin(A/B) %实验测量值
figure,plot(t,u,t,uR,'--','LineWidth',2)%开创图形窗口,画电流电压曲线
grid on,legend('总电压u','电阻电压u_R')%加网格,图例
xlabel('t'),title('总电压和电阻电压随时间的变化规律')%横坐标,标题
text(0,-0.5,['R=',num2str(R)])%电阻的值
figure,plot(u,uR)%开创图形窗口,画总电压和电阻电压的关系曲线
xlabel('u'),ylabel('u_R')%横纵坐标
title('李萨如图形'),grid on%标题,加网格
text(0,-1,['\phi_t_h=arctan(X_C/R)=',num2str(ph),'=',num2str(phd),'\circ'])%显示理论计算相位差
text(0,A,['A=',num2str(A)])%截距u=0,uR的值
text(B,B+0.1,['B=',num2str(B)])%uR的最大值
text(-8,1,['\phi_e_x_p=arcsin(A/B)=',num2str(phexp)])%显示实验测量相位差
text(-8,-1,['R=',num2str(R)])%电阻的值
hold on,stem([0,Um*sin(pi/2-ph)],[A,B],'--')%保持图像,画杆图,第一个【】对应画杆线的x坐标,第二【】对应着杆线的长度
[程序] RCR300_.m
%R=300,RC电路李萨如图形随频率的变化
clear,f=3000;c=0.22*10^(-6);%清除变量,频率,电容
Um=10;ph0=0;%电圧峰值,初相
w=2*pi*f;Xc=1/(w*c)%圆频率,容抗
R=300;z=sqrt(Xc^2+R^2),ph=atan(Xc/R)%电阻,阻抗,相位差
phd=ph*180/pi %相位差化为度数
Im=Um/z%电流峰值
t=linspace(0,0.005,100000);%时间向量
u=Um*sin(w*t+ph0);%总电压
uR=Im*R*sin(w*t+ph0+ph);%电阻电压
uC=Im*Xc*sin(w*t+ph0-pi/2+ph);%电容电压
A=Im*R*sin(ph);%截距u=0,uR的值
B=Im*R;%uR的最大值
phexp=asin(A/B) %实验测量值
figure,plot(t,u,t,uR,'--','LineWidth',2)%开创图形窗口,画电流电压曲线
grid on,legend('总电压u','电阻电压u_R')%加网格,图例
xlabel('t'),title('总电压和电阻电压随时间的变化规律')%横坐标,标题
text(0,-0.5,['f=',num2str(f)])%电阻的值
figure,plot(u,uR)%开创图形窗口,画总电压和电阻电压的关系曲线
xlabel('u'),ylabel('u_R')%横纵坐标
title('李萨如图形'),grid on%标题,加网格
text(0,-1,['\phi_t_h=arctan(X_C/R)=',num2str(ph),'=',num2str(phd),'\circ'])%显示理论计算相位差
text(0,A,['A=',num2str(A)])%截距u=0,uR的值
text(B,B+0.1,['B=',num2str(B)])%uR的最大值
text(-8,1,['\phi_e_x_p=arcsin(A/B)=',num2str(phexp)])%显示实验测量相位差
text(-8,-1,['f=',num2str(f)])%电阻的值
hold on,stem([0,Um*sin(pi/2-ph)],[A,B],'--')%保持图像,画杆图,第一个【】对应画杆线的x坐标,第二【】对应着杆线的长度