Delta Sigma Toolbox: un=uvar(u,N); - File Exchange

Code covered by the BSD License

Highlights from
Delta Sigma Toolbox

dsdemo4fig() Window layout for the dsdemo4() function.
1; end end
DocumentNTF(arg1,osr,f0,quadr... axis_handle = DocumentNTF(ntf|ABCD|mod_struct,osr=64,f0=0,quadrature=0) Plot the NTF's poles and zeros as well as its frequency-response
ESLselect(v,sy,dw,df) sv = ESLselect(v,sy,dw,df) Select the elements of a multi-element
H=evalMixedTF(tf,f,df) H=evalMixedTF(tf,f)
H=evalTFP(Hs,Hz,f)
LCObj(x,param,dbg) LCObj Compute the objective function
LCObj1(x,param,max_radius,dbg) function [objective,constraints] = LCObj1(x,param,max_radius,dbg)
LCoptparam2tf(x,param) LCoptparam2tf(...) Convert optimization parameters to transfer functions.
LCplotTF(H,L0,param) LCplotTF(H,L0,param)
NTF=clans(order,OSR,Q,rmax,op...
NTF=clans5(order,OSR,Q,rmax,o... Version of clans for MATLAB 5 or lower
NTF=clans6(order,OSR,Q,rmax,o... Version of clans for MATLAB >=6
PlotExampleSpectrum(mod_struc... PlotExampleSpectrum(ntf|mod_struct,M=1,osr=64,f0=0,quadrature=0)
[ABCDs,umax,S]=scaleABCD(ABCD...
[H,L0,ABCD,k]=LCparam2tf(para...
[f,g]=dsclansObj(x,order,OSR,... Objective function for clans.m
[f1,f2,info]=designHBF(fp,del...
[f1_saved,f2_saved,info]=desi...
[f1_saved,f2_saved,info]=desi...
addPIS Add the PosInvSet subdirectory of the toolbox to the current path
axisLabels(range,incr) function s = axisLabels(range,incr)
bilogplot(V,f0,fbin,x,y,fmt) bilogplot(V,f0,fbin,x,y,fmt) Plot two side-by-side spectra
bplogsmooth(X,tbin,f0)
bquantize(x,nsd,abstol,reltol)
bunquantize(q) Calculate the value corresponding to a bidirectionally quantized quantity.
calculateQTF(ABCDr)
calculateSNR(hwfft,f,nsig) snr = calculateSNR(hwfft,f,nsig=1) Estimate the signal-to-noise ratio,
calculateTF(ABCD,k)
cancelPZ(zpk1, tol) zpk2 = cancelPZ(zpk1, tol=1e-6) Cancel zeros/poles in a zpk system
changeFig(fontsize,linewidth,... changeFig(fontsize,linewidth,markersize) Change the settings
circ_smooth(x,n)
delay(x,n)
designLCBP(n,OSR,opt,Hinf,f0,...
designLCBP6(n,OSR,opt,Hinf,f0... Modified designLCBP for use with latest optimization toolbox
dotplot(p,fmt) function dotplot(p,fmt)
ds_f1f2(OSR,f0,complex_flag);
ds_freq(osr,f0,quadrature) f = ds_freq(osr=64,f0=0,quadrature=0) Frequency vector suitable for plotting the frequency response of an NTF
ds_hann(n)
ds_optzeros( n, opt )
ds_quantize(y,n)
ds_synNTFobj1(x,p,osr,f0) y=ds_synNTFobj1(x,p,osr,f0) Objective function for synthesizeNTF()
dsclansNTF(x,order,rmax,Hz) Conversion of clans parameters into a NTF.
dscut(p1,y1,p2,y2)
dsdemo4(action) A GUI-based demonstration of delta-sigma with sound.
dsexample1 Design example for a lowpass modulator.
dsexample2 Design example for a bandpass modulator.
dsexample3 Design example for a continuous-time lowpass DS ADC
dsexample4 Example quadrature bandpass modulator
dsexpand(s,c,k,n,o)
dsisPlot(dbg,itn,order,x,s,e,... Show some pretty pictures
dsmap(u,ABCD,nlev,x,e,v)
dssplit2d(u,ABCD,p)
edgeplot(e,s,fmt) function edgeplot(e,s,fmt) Plot edges e with format fmt
evalF0(f1,z,phi) F0 = evalF0(f1,z,phi) Calculate the values of the F0 (prototype) filter
evalF1(f1,z,phi) F1 = evalF1(f1,z,phi) Calculate the values of the F1 filter
evalRPoly(roots,x,k)
evalTF(tf,z)
exampleHBF(n) Coefficients from toolbox v2.0 documentation
find2dPIS(u,ABCD,options) function s = find2dPIS(u,ABCD,options)
findPIS(u,ABCD,nlev,options)
findPattern(N,OSR,NTF,ftest,A... [data snr] = findPattern(N=1024,OSR=64,NTF,ftest,Atest,f0,nlev,quadrature,debug)
flattenStruct(s) flattenStruct(s) Puts each field of s into the current workspace.
frespF1(f1,f,phi) frespF1(f1,f) Plot/calculate the frequency response of the F1 filter
frespHBF(f,f1,f2,phi,fp,msg)
hull2d(p)
i=outsideConvex(x,n,o,tol) function i=outsideConvex(x,n,o,tol<0>). Determine which points are outside
impL1(arg1,n)
infnorm(H) [Hinf,fmax] = infnorm(H) Find the infinity norm of a z-domain transfer function.
l1norm(H) h1 = l1norm(H) Compute the l1-norm of a z-domain transfer function.
logsmooth(X,inBin,nbin,n)
lollipop(x,y,color,lw,ybot) lollipop(x,y,color,lw,ybot) Plot lollipops (o's and sticks)
mapABCD(ABCD,form)
mapCtoD(sys_c,t,f0)
mapQtoR(Z) A = mapQtoR(Z) Convert a quadrature matrix into its real equivalent.
mapR2Q( ABCDr )
mod1()
mod2()
outconvex2d(x,p) function out = outconvex2d(x,p)
padb(x, n, val)
padl(x, n, val)
padr(x, n, val)
padt(x, n, val)
partitionABCD(ABCD, m);
peakSNR(snr,amp) [peak_snr,peak_amp] = peakSNR(snr,amp) Find the snr peak by fitting
plotPZ(H,color,markersize,lis... function plotPZ(H,color='b',markersize=5,list=0)
plotSpectrum(X,fin,fmt) plotSpectrum(X,fin,fmt) Plot a smoothed spectrum
plotUsave(sv) Plot the elemet usage for a multi-elemet DAC
polyplot(p,fmt) function polyplot(p,fmt)
predictSNR(ntf,R,amp,f0)
printmif(file,size,font,fig) printmif(file,size,font,fig) Print graph to an Adobe Illustrator file
pulse(S,tp,dt,tfinal,nosum)
qhull(points,str)
realizeNTF(ntf,form,stf)
realizeNTF_ct( ntf, form, tda...
realizeQNTF(ntf,form,rot,bn)
rms(x,no_dc)
rmsGain(H,f1,f2,N)
selectQESL(v,sy)
simulateDSM(u,arg2,nlev,x0)
simulateESL(v,mtf,M,dw,sx0)
simulateHBF(x,f1,f2,mode) y = simulateHBF(x,f1,f2,mode=0) Simulate a Saramaki half-band filter.
simulateQDSM(u,arg2,nlev,x0)
simulateQESL(v,mtf,M,sx0)
simulateQSNR(ntf,R,amp,f0,nle...
simulateSNR(arg1,osr,amp,f0,n...
sinc_decimate(x,m,r)
stuffABCD(a,g,b,c,form)
synthesizeChebyshevNTF(order,... ntf = synthesizeChebyshevNTF(order=3,OSR=64,opt,H_inf=1.5,f0=0)
synthesizeNTF(order,osr,opt,H...
synthesizeNTF0(order,OSR,opt,... Determine the zeros.
synthesizeNTF1(order,osr,opt,... Determine the zeros.
synthesizeQNTF(order,OSR,f0,N...
thermometer(x,m)
un=uvar(u,N);
undbp(x)
undbv(x)
v=undbm(p,z)
y=dbm(v,R)
y=dbp(x)
y=dbv(x)
y=leftof(p,a,b) function y=leftof(p,a,b)
y=nabsH(w,H) nabsH computes the negative of the absolute value of H
y=sgn(x)
yRange(1):dy:yRange(2); s...
zinc(f,m,n) mag = zinc(f,m=64,n=1) Calculate the magnitude response
Contents.m
dsdemo1.m
dsdemo2.m
dsdemo3.m
dsdemo5.m
dsdemo6.m
dsdemo7.m
dsdemo8.m
dsistest.m
View all files

from Delta Sigma Toolbox by Richard Schreier
High-level design and simulation of delta-sigma modulators

un=uvar(u,N);

function un=uvar(u,N);
%un = uvar(u,N);	Compute a bounded sequence that has 90% of its 
%values at the extremes.

u1 = u(1); u2=u(2);
un = u1(ones(1,N));
r = rand(1,N);
ri2 = r > 0.55;
ri1 = r > 0.45 & ~ri2;
un(ri1) = u1 + (u2-u1)*10*(r(ri1)-.45);
un(ri2) = u2(ones(1,sum(ri2)));

Highlights from Delta Sigma Toolbox

Highlights from
Delta Sigma Toolbox