lamda=410:15:560;
v=((lamda./1e9).^(-1)).*(2.99792458*1e8).*1e-14;
v=v';
U0=-[1.654 1.572 1.492 1.408 1.336 1.252 1.160 1.030 0.952 0.872 0.804]';
Ud=-[1.648 1.560 1.488 1.400 1.312 1.224 1.140 1.000 0.940 0.864 0.788]';
X=[ones(11,1) v];
[b1]=regress(U0,X);
[b2]=regress(Ud,X);
z1=b1(1)+b1(2)*v;
z2=b2(1)+b2(2)*v;
plot(v,U0,'ko',v,z1,'r');
grid on
xlabel("光的频率v(10^1"+" ^4"+"Hz)");
ylabel("电流为0对应的电压U0");
title("U0~v拟合直线");
txt1 = ['拟合系数d1为:' num2str(b1(1)) ];
txt2 = ['拟合系数k1为:' num2str(b1(2)) ];
text(5.1,-1.5,txt1);
text(5.1,-1.56,txt2);
plot(v,Ud,'ko',v,z2,'r');
grid on
xlabel("光的频率v(10^14Hz)");
ylabel("截止电压Ud");
title("U0~v拟合直线");
txt3 = ['拟合系数d2为:' num2str(b2(1)) ];
txt4 = ['拟合系数k2为:' num2str(b2(2)) ];
text(5.1,-1.5,txt3);
text(5.1,-1.56,txt4);
grid on
xlabel("光的频率v(10^1"+" ^4"+"Hz)");
ylabel("截止电压Ud");
title("Ud~v拟合直线");

