noiseBudgets.m

fileName = 'noiseBudgets';
%%
checkCircuit(fileName);
htmlPage('Circuit data');
img2html([fileName, '.svg'], 760);
netlist2html(fileName);
defPar('C_A', CA);
Cf = PR.R7.value*CA/PR.R5.value/(PR.R6.value + PR.R7.value);
defPar('C_f', Cf);
defPar('R_o', PR.R6.value);
defPar('R_ell', PR.R7.value);
params2html();
%%
htmlPage('Noise budgets');
text2html(['Determination of the budgets for the equivalent-input ',...
    'noise sources of the controller.']);
simType('numeric');
gainType('vi');
dataType('noise');
source('V1');
detector('V_ell');
sourceNoiseSpectrum = getInoise(execute());
text2html('The source-referred noise spectrum equals:');
syms('S_v_tot', 'F', 'S_vmax', 'S_imax', 'S_v', 'S_i');
eqn2html(S_v_tot, sourceNoiseSpectrum);
text2html(['A show stopper value $S_{\\mathrm{vmax}}$ for $S_v$ is found ' ...
    'by setting $S_i=0$ and solve $S_v$ from the specification of the ' ...
    'floor noise at a high frequency (30MHz).']);
Svmax = vpa(solve(subs(sourceNoiseSpectrum - PR.R8.value, [S_i, F], [0, 30E6])), 3);
text2html('In this way we obtain in $\\mathrm{V/\\sqrt{Hz}}$:');
eqn2html(S_vmax, sqrt(Svmax));
text2html(['A show stopper value $S_{\\mathrm{imax}}$ for $S_i$ is found ' ...
    'by setting $S_v=0$ and solve $S_i$ from the specification of the ' ...
    '1/f noise at a low frequency (100kHz).']);
Simax = vpa(solve(subs(sourceNoiseSpectrum - PR.R9.value, [S_v, F], [0, 100E3])), 3);
text2html('In this way we obtain in $\\mathrm{A/\\sqrt{Hz}}$:');
eqn2html(S_imax, sqrt(Simax));
script2html('noiseBudgets');
stophtml();