K1_K4.m

clear all;
close all;

syms x K y dy;
y=3435/log((10000*x)/((5-x)*0.0992));
dy=diff(y);
K=symfun(abs((x*dy)/y),(x));
da=0.0000001;
a=da;
k=zeros(2000,2);

for i=1:length(k)
    
    k(i,1)=a;
    k(i,2)=K(a);
    a=a+da;
    
end

c=zeros(length(k),1);
c(:)=1;

scrsz=get(0,'Screensize');
figure('Position',[1 scrsz(4) scrsz(3) scrsz(4)])
subplot(2,1,1);
plot(k(:,1),k(:,2),k(:,1),c,'r');
xlim([0 0.00019]);
ylim([0 1.9]);
xlabel('Voltage / V');
ylabel('Condition');
title('Condition of the temerature algorithm');
legend('Condition','Boundary value','Location','Best');
grid on;

syms x K y dy;
y=3435/log(x);
dy=diff(y);
K=symfun(abs((x*dy)/y),(x));
dV=0.0000001;
V=dV;
k=zeros(2000,2);

for i=1:length(k)
    
    R=10000*V/((5-V)*0.0992);
    k(i,1)=V;
    k(i,2)=K(R);
    V=V+dV;
    
end

c=zeros(length(k),1);
c(:)=1;

subplot(2,1,2);
plot(k(:,1),k(:,2),k(:,1),c,'r');
xlim([0 0.00019]);
ylim([0 1.9]);
xlabel('Voltage / V');
ylabel('K4');
title('Condition of the 4th step of the temerature algorithm --- Range: 0.1 uV - 0.2 mV');
legend('Condition','Boundary value','Location','Northeast');
grid on;