Pastebin launched a little side project called, check it out ;-) Want more features on Pastebin? Sign Up, it's FREE!

Szymon Bęczkowski's PhD dissertation

By: owiecc on Jan 22nd, 2014  |  syntax: Latex  |  size: 248.41 KB  |  views: 143  |  expires: Never
download  |  raw  |  embed  |  report abuse  |  print
Text below is selected. Please press Ctrl+C to copy to your clipboard. (⌘+C on Mac)
  1. %
  2. %  Power converter for LED based intelligent light sources
  3. %
  4. %  Created by Szymon Bęczkowski on 2008-10-15.
  5. %
  7. % preamble (fold)
  9. %!TEX TS-program = xelatex
  10. %!TEX encoding = UTF-8 Unicode
  12. \documentclass[12pt,a4paper,titlepage,onecolumn,openright,twoside]{report}
  14. % inverted text color begin
  15. % \usepackage{empheq} // not necessary
  16. % \everymath{\color{white}}
  17. % \everydisplay{\color{white}}
  18. % inverted text color end
  20. \usepackage[colorlinks=true,citecolor=black,linkcolor=black,pdftitle={Power converter for LED based intelligent light sources}]{hyperref} % clickable references
  22. \setlength{\parskip}{0pt} % no spacing between paragraphs
  24. \usepackage{amsmath}
  25. \usepackage{mathspec}
  26. \usepackage{microtype} % requires microtype 2.5
  27. \LoadMicrotypeFile{pmn}
  28. \usepackage{xltxtra,eukdate}
  29. \defaultfontfeatures{Mapping=tex-text}
  30. \setmainfont[Numbers=OldStyle]{Minion Pro}
  31. \setallsansfonts[Numbers={OldStyle,Proportional},Scale=MatchLowercase]{Minion Pro}
  32. \setallmonofonts[Numbers=OldStyle,Scale=MatchLowercase]{Minion Pro}
  33. \setmathsfont(Digits,Latin)[Scale=MatchLowercase]{Minion Pro}
  34. \setmathsfont(Greek)[Scale=MatchLowercase]{Minion Pro}
  35. \setmathrm{Minion Pro}
  36. \exchangeforms{phi}
  37. \setminwhitespace[750]
  39. % control widows and orphans
  40. \raggedbottom
  41. \widowpenalty=500
  42. \clubpenalty=500
  44. % map ligatures to plaintext for searching and copy&paste
  46. % Smaller captions
  47. \usepackage[hang,small]{caption}
  48. %\DeclareCaptionLabelFormat{bf-parens}{(\textbf{#2})}
  49. %\captionsetup{labelformat=bf-parens,labelsep=quad}
  51. % Set text colour
  52. \usepackage{color} % use: \textcolor{BrickRed}{text}
  54. % Title logo
  55. \usepackage{titlepic}
  57. % This is now the recommended way for checking for PDFLaTeX:
  58. \usepackage{ifpdf}
  60. %\newif\ifpdf
  61. %\ifx\pdfoutput\undefined
  62. %\pdffalse % we are not running PDFLaTeX
  63. %\else
  64. %\pdfoutput=1 % we are running PDFLaTeX
  65. %\pdftrue
  66. %\fi
  68. \ifpdf
  69. \usepackage[pdftex]{graphicx}
  70. \else
  71. \usepackage{graphicx}
  72. \fi
  74. % Hyphenation
  75. \hyphenation{pro-pos-ing con-gress in-ter-na-tion-al}
  77. % Small caps
  78. \newcommand{\dc}{\textsc{dc}}
  79. \newcommand{\ac}{\textsc{ac}}
  80. \newcommand{\uv}{\textsc{uv}}
  81. \newcommand{\ir}{\textsc{ir}}
  82. \newcommand{\cct}{\textsc{cct}}
  83. \newcommand{\led}{\textsc{led}}
  84. \newcommand{\hid}{\textsc{hid}}
  85. \newcommand{\rgb}{\textsc{rgb}}
  86. \newcommand{\rgbw}{\textsc{rgbw}}
  87. \newcommand{\cie}{\textsc{cie}}
  88. \newcommand{\xyz}{\textsc{xyz}}
  89. \newcommand{\cri}{\textsc{cri}}
  90. \newcommand{\tff}{\textsc{tff}}
  91. \newcommand{\ffb}{\textsc{ffb}}
  92. \newcommand{\pcb}{\textsc{pcb}}
  93. \newcommand{\pwm}{\textsc{pwm}}
  94. \newcommand{\dsp}{\textsc{dsp}}
  95. \newcommand{\pcm}{\textsc{pcm}}
  96. \newcommand{\am}{\textsc{am}}
  97. \newcommand{\ffbtff}{\textsc{ffb&tff}}
  98. \newcommand{\ccfb}{\textsc{ccfb}}
  99. \newcommand{\fwhm}{\textsc{fwhm}}
  100. \newcommand{\mosfet}{\textsc{mosfet}}
  101. \newcommand{\srgb}{{\footnotesize{}s}\textsc{rgb}}
  102. \newcommand{\srh}{\textsc{srh}}
  103. \newcommand{\matlab}{\textsc{matlab}}
  104. \newcommand{\ivxyz}{\textsc{ivxyz}}
  105. \newcommand{\dcdc}{\textsc{dc-dc}}
  107. \newcommand{\abbr}[1]{{\fontspec[Numbers=OldStyle]{Minion Pro}\textsc{#1}\fontspec[Numbers=Lining]{Minion Pro}}}
  108. \def\mathbi#1{{\fontspec[Numbers=Lining]{Minion Pro}\textbf{\em #1}}}
  110. % Math
  111. \newcommand{\ud}{\mathrm{d}} % straight "d" for integrals
  112. \newcommand{\deriv}[2]{\frac{\ud{}#1}{\ud{}#2}}
  113. \newcommand{\pderiv}[2]{\frac{\partial{}#1}{\partial{}#2}}
  114. \newcommand{\avg}[1]{\big< #1 \big>} % average
  116. % Typographic
  117. \newcommand{\hsp}{\hspace{1pt}}% hair space
  118. \frenchspacing % single space after sentence
  119. \setlength{\parindent}{1em}
  121. % Pretty sidenotes \marginpar{Sidenote}
  122. \usepackage{mparhack} % Fixes marginpar appearing on the wrong side of the page
  123. \let\oldmarginpar\marginpar
  124. \renewcommand\marginpar[1]{\-\oldmarginpar[\raggedleft\footnotesize\textsc{#1}]{\raggedright\footnotesize\textsc{#1}}}
  126. % TOC styling
  127. \usepackage[titles]{tocloft} % chage style of TOC, [titles] to control the style of TOC heading
  129. % distance between title and page number
  130. \renewcommand{\cftchapleader}{\hspace{8pt}}
  131. \renewcommand{\cftsecleader}{\hspace{8pt}}
  132. \renewcommand{\cftsubsecleader}{\hspace{8pt}}
  133. % subsections italic
  134. \renewcommand{\cftsubsecfont}{\it}
  135. % Trick to flush page numbers left within their box
  136. \makeatletter
  137. \renewcommand{\@pnumwidth}{0.1em}
  138. \makeatother
  139. % page number colours
  140. \renewcommand{\cftchappagefont}{\color[rgb]{0.715,0.05,0.07}}
  141. \renewcommand{\cftsecpagefont}{\color[rgb]{0.715,0.05,0.07}}
  142. \renewcommand{\cftsubsecpagefont}{\color[rgb]{0.715,0.05,0.07}}
  143. % no dots fill
  144. \renewcommand{\cftchapafterpnum}{\cftparfillskip}
  145. \renewcommand{\cftsecafterpnum}{\cftparfillskip}
  146. \renewcommand{\cftsubsecafterpnum}{\cftparfillskip}
  147. % indents
  148. \cftsetindents{chapter}{0pt}{24pt}
  149. \cftsetindents{section}{0pt}{24pt}
  150. \cftsetindents{subsection}{-1000pt}{1024pt}
  151. % vertical space
  152. \setlength\cftbeforechapskip{14pt}
  154. % Chapter, section and subsection headings styling
  155. \makeatletter
  156. % Chapter 1, Chapter 2, ... heading styling
  157. \def\@makechapterhead#1{%
  158.   { \parindent \z@ \raggedright \normalfont
  159.    \ifnum \c@secnumdepth >\m@ne  
  160.        \huge\thechapter\hspace{8pt}\fontspec[Numbers=Lining]{Minion Pro Subhead}#1\par\nobreak\vspace{28pt}
  161.        \par\nobreak
  162.    \fi
  163.    \interlinepenalty\@M
  164.  }
  165. }
  166. % List of abbreviations, Bibliography, ... heading styling
  167. \def\@makeschapterhead#1{%
  168.   { \parindent \z@ \raggedright
  169.    \normalfont
  170.    \interlinepenalty\@M
  171.    \huge\fontspec[Numbers=Lining]{Minion Pro Subhead} #1\par\nobreak\vspace{28pt}
  172.  }
  173. }
  175. \renewcommand{\section}{\@startsection%
  176.         {section}%                      % name
  177.         {1}%                            % level
  178.         {0pt}%                          % indent
  179.         {-\baselineskip}%       % before skip
  180.         {\baselineskip}%        % after skip
  181.         {\normalfont\sc}}       % the style
  182. \renewcommand{\subsection}{\@startsection%
  183.         {subsection}%           % name
  184.         {2}%                            % level
  185.         {0pt}%                          % indent
  186.         {-\baselineskip}%       % before skip
  187.         {\baselineskip}%        % after skip
  188.         {\normalfont\it}}       % the style
  189. \makeatother
  191. % Titlepage info
  192. \title{\color[rgb]{0.196,0.196,0.204}Control and driving methods for \led{} based intelligent light~sources}
  193. \author{\color[rgb]{0.196,0.196,0.204}Szymon Bęczkowski}
  194. \date{\color[rgb]{0.196,0.196,0.204} }
  195. \titlepic{\vspace{8cm}\includegraphics{graphics/ETlogo.eps}}
  196. % preamble (end)
  198. \begin{document}
  200. \ifpdf
  201. \DeclareGraphicsExtensions{.pdf, .jpg, .tif}
  202. \else
  203. \DeclareGraphicsExtensions{.eps, .jpg}
  204. \fi
  206. %%%%%%%%%%%%%%%
  207. %%   Title   %%
  208. %%%%%%%%%%%%%%%
  210. \maketitle
  212. %%%%%%%%%%%%%%%
  213. %% Abstract  %%
  214. %%%%%%%%%%%%%%%
  216. \newpage
  217. \thispagestyle{empty}
  218. \mbox{}
  219. \begin{abstract}
  220. High power light-emitting diodes allow the creation of luminaires capable of generating saturated colour light at very high efficacies. Contrary to traditional light sources like incandescent and high-intensity discharge lamps, where colour is generated using filters, \led{}s use additive light mixing, where the intensity of each primary colour diode has to be adjusted to the needed intensity to generate specified colour.
  222. The function of \led{} driver is to supply the diode with power needed to achieve the desired intensity. Typically, the drivers operate as a current source and the intensity of the diode is controlled either by varying the magnitude of the current or by driving the \led{} with a pulsed current and regulate the width of the pulse. It has been shown previously, that these two methods yield different effects on diode's efficacy and colour point.
  224. A hybrid dimming strategy has been proposed where two variable quantities control the intensity of the diode. This increases the controllability of the diode giving new optimisation possibilities. It has been shown that it is possible to compensate for temperature drift of white diode's colour point using hybrid dimming strategy. Also, minimisation of peak wavelength shift was observed for InGaN diodes.
  226. Control of trichromatic luminaires, dimmed with either pulse-width modulation or amplitude modulation, cannot be optimised. Introduction of hybrid dimming mechanism creates three additional degrees of freedom therefore luminaire parameters such as luminous flux, efficacy and colour quality can be maximised. Simulations show that the gamut of the device can be increased, especially in the cyan colour range for \rgb{} luminaires.
  228. A current-voltage model of light-emitting diode is presented. It utilises the fact that instantaneous values of diode's current and voltage correspond uniquely to a set of diode's colorimetric properties, like tristimulus values. This model can be used for colorimetric feedback in colour control loop. The model was created in thermal steady-state conditions and its validity has been tested with a diode driven with a pulsed current. The model can also be used to create highly accurate luminaire model.
  230. Finally, a dual interleaved buck converter has been proposed for driving high power light-emitting diodes. Interleaving two converters lowers the output ripple current thus lowering the requirement on the output capacitor. It has been shown that at the expense of cost and increased complexity an efficient design can be created for supplying high current to \led{}s without the need for electrolytic capacitors.
  231. \end{abstract}
  232. \newpage
  233. \mbox{}
  234. \thispagestyle{empty}
  235. \newpage
  237. %%%%%%%%%%%%%%%
  238. %%    TOC    %%
  239. %%%%%%%%%%%%%%%
  241. \setcounter{tocdepth}{2}
  242. \setcounter{page}{1}
  243. \tableofcontents
  245. %%%%%%%%%%%%%%%
  246. %%   Abbr.   %%
  247. %%%%%%%%%%%%%%%
  248. % abbreviations (fold)
  249. \chapter*{List of abbreviations, symbols and physical constants}
  250. \begin{tabular}{p{1.5cm}l} % List of abbreviations
  251. \ac & alternating current \\
  252. \am & amplitude modulation (same as \textsc{ccr}) \\
  253. \dc & direct current \\
  254. \textsc{ic}  & integrated circuit \\
  255. \textsc{iv} & current-voltage \\
  256. \textsc{pi}  & proportional–integral (controller) \\
  257. \textsc{pc}  & phosphor-converted (diode) \\
  258. \textsc{adc} & analog-digital conversion \\
  259. \textsc{ccr} & continous current reduction (same as \textsc{am}) \\
  260. \cct & correlated colour temperature \\
  261. \textsc{cff} & critical flicker frequency \\
  262. \cie & Commission Internationale de l'Eclairage \\
  263. \textsc{cpu} & central processing unit \\
  264. \cri & colour rendering index \\
  265. \textsc{dac} & digital-analog conversion \\
  266. \textsc{dcr} & \dc{} resistance \\
  267. \dsp & digital signal processor \\
  268. \textsc{emi} & electro-magnetic interference \\
  269. \ffb & flux feedback \\
  270. \hid & high-intensity discharge \\
  271. \led & light-emitting diode \\
  272. \textsc{pcb} & printed circuit board \\
  273. \textsc{pfc} & power factor correction \\
  274. \pwm & pulse width modulation \\
  275. \textsc{rgb} & red, green and blue \\
  276. \srh & Shockley-Read-Hall recombination \\
  277. \tff & temperature feed forward \\
  278. \textsc{tim} & thermal interface material \\
  279. \abbr{tcs} & test colour sample \\
  280. \textsc{vrm} & voltage regulation module \\
  281. \ccfb & colour coordinates feedback \\
  282. \fwhm & full width at half maximum \\
  283. \mosfet & metal-oxide-semiconductor field-effect transistor
  284. \end{tabular} \\ \newline  \newline
  285. \begin{tabular}{p{1.5cm}lc} % List of symbols
  286. $λ$ & wavelength & m \\
  287. $ν$ & frequency of optical radiation & Hz \\
  288. $η$ & efficiency & p.u. \\
  289. $η_{lum}$ & luminous efficacy & lm/W \\
  290. $t$ & time & s \\
  291. $τ$ & time constant & s \\
  292. $F$ & luminous flux & lm \\
  293. $A$ & area & $\phantom{^2}\mathrm{m^2}$ \\
  294. $C$ & capacitance & F \\
  295. $L$ & inductance & H \\
  296. $R$ & resistance & Ω \\
  297. $P$ & power & W \\
  298. $P_{rad}$ & radiometric power & W \\
  299. $i$ & current & A \\
  300. $v$ & voltage & V \\
  301. $d$ & duty cycle & p.u. \\
  302. $f$ & frequency & Hz \\
  303. $T$ & temperature & °C, K \\
  304. $R_i$ & colour rendering index for \textit{i}-th sample & — \\
  305. $R_a$ & general colour rendering index & — \\
  306. $E_g$ & bandgap energy & eV \\
  307. $C_{th}$ & thermal capacitance & J/K \\
  308. $R_{th}$ & thermal resistance & K/W \\
  309. \textsc{ct/cct} & colour temperature/corelated colour temperature & K \\
  310. \textsc{fwhm} & full width at half maximum & eV, nm \\
  311. $x,y$ & colour coordinates in \abbr{cie 1931} colour space & — \\
  312. $u,v$ & colour coordinates in \abbr{cie 1960} colour space & — \\
  313. $u'\!,v'$ & colour coordinates in \abbr{cie 1976} ($L'\!,u'\!,v'$) colour space & — \\
  314. $a^{*},b^{*}$ & colour coordinates in \abbr{cie 1976} ($L^{*}\!,a^{*}\!,b^{*}$) colour space & — \\
  315. $\overline{x},\overline{y},\overline{z} $ & colour matching functions & p.u. \\
  316. $X,Y,Z$ & tristimulus values & W, lm \\
  317. \end{tabular} \\ \newline  \newline
  318. \begin{tabular}{l@{~=~}ll} % List of constants
  319. $e$ & $1.6022\cdot{}10^{-19}\thinspace{}\mathrm{C}$ & elementary charge \\
  320. $c$ & $2.9979\cdot{}10^{8}\thinspace{}\mathrm{m/s}$ & speed of light in vacuum \\
  321. $h$ & $6.6261\cdot{}10^{-34}\thinspace{}\mathrm{Js}$ & Planck constant \\
  322. $h$ & $4.1357\cdot{}10^{-15}\thinspace{}\mathrm{eVs}$ & Planck constant \\
  323. $k$ & $1.3807\cdot{}10^{-23}\thinspace{}\mathrm{J/K}$ & Boltzman constant \\
  324. $k$ & $8.6175\cdot{}10^{-5}\thinspace{}\mathrm{eV/K}$ & Boltzman constant \\
  325. \end{tabular}
  326. % abbreviations (end)
  328. %%%%%%%%%%%%%%%
  329. %% Chapter 1 %%
  330. %%%%%%%%%%%%%%%
  331. \chapter{Introduction} % (fold)
  332. \label{cha:introduction}
  333. %
  334. \section[Background]{background}
  335. Intelligent lighting such as projectors and moving heads are used, whenever there is a need for a light source that can be tuned to individual needs during lamp operation. Typically these lamps can have programmable light colour, intensity, focus and sometimes the ability to project special shapes such as patterns or logos. Intelligent lighting is commonly used in theatres, concerts halls and clubs.
  336. \begin{figure}[h]% fig:Projectors
  337.   \centering
  338.   \includegraphics{graphics/Introduction/Projectors.pdf}
  339.  \caption{Martin smart\textsc{mac} moving head (left) and Exterior~1200 image projector (right).}
  340.  \label{fig:Projectors}
  341. \end{figure}
  343. \noindent{}Two examples of such a lamps are shown in figure~\ref{fig:Projectors}. Moving head lamps are designed with the ability to pan and tilt the light beam. Projectors, on the other hand, are fixed at specific target during installation. Both lamp types can be designed to provide a narrow or diffused beam of light, depending on the need.
  345. Intensity of the light provided by a luminaire must be high, so typically high-intensity discharge lamps are used as a light source. A compact passive-cooled luminaire like the smart\textsc{mac} uses a 150\thinspace{}\textsc{w} \hid{} lamp while the Exterior~1200 uses a 1200\thinspace{}\textsc{w} metal halide lamp yielding ca.~10 and 92\thinspace{}klm, respectively. Because the output light beam shape is very different from the light source profile, the optics inside the lamp must collect and collimate. A portion of the light is lost in the process. The output light flux for aforementioned luminaires is ca.~2700\thinspace{}lm for the smart\textsc{mac} and 20–35\thinspace{}klm for the Exterior 1200 (depending on the configuration). As can be seen, the efficiency of the luminaire is around 30\thinspace{}\%, not including the power needed to drive the control electronics.
  347. \begin{figure}[t]% fig:ProjectorsMACIII
  348.   \centering
  349.   \includegraphics{graphics/Introduction/ProjectorsMACIII.pdf}
  350.  \caption{Martin \textsc{mac iii} moving head luminaire and its blow-out diagram showing the 1.5\thinspace{}kW discharge light source, dimming mechanism and the colour wheel.}
  351.  \label{fig:ProjectorsMACIII}
  352. \end{figure}
  354. Above flux values are valid for generation of white light. Colour light is created by inserting a colour filter in the light beam (fig.~\ref{fig:ProjectorsMACIII}). The transmissivity of the filters depends on the desired colour and its saturation. Saturated red and blue colour filters can have transmission as low as few percent. This lowers the overall efficacy of the luminaire. Moreover, it is hard to create filters with sharp characteristics to create saturated colours. Transition from one colour to another requires mechanical change of the filter in the path of the light. Depending on the filter arrangement, it may be necessary to transition through a different colours to obtain the desired colour.
  356. Dimming of high power \hid{} lamps is done by either a high frequency lamp ballast (limited dimming range) or inserting a specially shaped shutter in the light path (fig.~\ref{fig:ProjectorsMACIII}). When the light is to be turned off for a short period of time the full-blocking shutter is placed in front of the light source. This is done because metal halide lamps cannot ignite unless they are cold because the high pressure obtained during normal operation prevents the lamp to restart. This requires a few minute cool down after the lamp was turned off.
  358. Lifetime of \hid{} lamps is limited to a few thousands hours. At the end of their lifetime high-intensity discharge lamps exhibit phenomenon called cycling, when the lamps periodically turn on and off because of the ageing effects. Sometimes the stress caused by off-on cycles or the increased pressure of the gas inside the lamp can lead to an explosion of the lamp tube, containing small amount of mercury.
  360. Constantly increasing intensity of \led{} sources already allows their use in low power (watts–tens of watts) applications such as replacement for conventional halogen and incandescent light sources in general lighting. Slowly, they enter medium power (hundreds of watts) applications like moving heads. \led{}s do not contain any toxic chemicals and have much longer lifetime than \hid{} lamps. Long lifetime decreases the service costs of the luminaire. % high mounting, keeping track of lifetime, lamp may disable to prevent blow
  362. Compared to incandescent and discharge lamps, \led{}s emit a narrow spectrum light that produces saturated colours. This means the white light cannot be directly generated by a single \led{} die. White light can be generated using either two or more pure colour \led{}s or a phosphor coated blue \led{}. Similarly to obtain any other colour in the device gamut, a combination of primary colours must be used. Generated colour range is not limited to the number of colour filters used, but on the accuracy of a colour control scheme. High bandwidth of intensity control allows for fast transition between colours and for strobing without the need for mechanical shutter. Removing mechanical parts like colour wheel or shutters decreases luminaire's weight and size and, at the same time, increases the reliability of the device.
  364. Using light-emitting sources in luminaires comes at a price of increased complexity in colour control, thermal management, optical system and power electronics controlling \led{} light sources. Colour control system depends on the accuracy needed in the luminaire. It can be as simple as open loop control or an advanced feedback loop compensating for diode's changing colour point. \led{}'s properties depend on the junction temperature therefore the diodes should be kept as cool as possible. High power \led{}s can dissipate even 80\thinspace{}\% of their input power in a few square millimetre area chip. The thermal management system should remove the heat from the \led{} structure as efficiently as possible or the junction temperature will increase and the light intensity and lifetime of the \led{} will drop.
  366. \marginpar{inled project}This work is part of an \textsc{inled} project (Intelligent Light Emitting Diodes), a collaboration between Martin Professional, Aalborg University (Institute for Energy Technology and Institute for Nano-Physics) and the Danish National Advanced Technology Foundation (Højteknologifonden). The objective of the project is to increase the knowledge in \led{} based luminaire technology. The knowledge gained from the project should be useful for designing substitutes for existing  technologies. The project focuses on creating: optics with nano scale coatings that would maximise the light transmission from \led{}s to the output port, efficient power electronic drivers supplying power to the diodes and a heat management system for the luminaire.
  368. \marginpar{motivation}The aim of the work shown in this dissertation is to research the driving of high power light-emitting diodes and to create an intelligent \led{} driver that would fully utilise the benefits of solid-state lighting technology.  
  370. \section[Related work]{related work} % Related work
  371. A research of the state of the art in solid-state lighting has been conducted to show current solution of the existing problems and show potentially beneficial research areas. Three main topics are covered: properties of \led{} devices, control of luminaires and methods of driving light-emitting diodes.
  373. \subsection{\led{} control}
  374. With the introduction of high power blue GaN \led{}s by Nakamura of Nichia Corporation \led{}s finally gained the potential to enter general illumination market. Four basic diode colours: red, green, blue and white became available on the market as high power devices in the mid 2000s. At this point scientist started investigating their driving properties. At that time  phosphor-converted white \led{}s had high enough efficacy for practical lighting applications. Lighting Research Center at Rensselaer Polytechnic Institute pioneered in colorimetric research of high-power devices. \marginpar{dimming properties}Dyble~et~al.~\cite{Dyble2005} analysed the chromaticity shifts of \textsc{pc}~white \led{} under pulse width modulation and amplitude modulation dimming. A year later research on red, green, blue and white diodes was conducted by Gu~et~al.~\cite{Gu2006} on the effects of different dimming schemes on diodes' spectra, luminous flux and efficacy. In 2007, Manninen and Orreveteläinen~\cite{Manninen2007} from Helsinki University of Technology published their work on spectral and thermal behaviour of AlGaInP diodes under \pwm{} dimming. General consensus from these three research papers was that \pwm{} dimming method provides linear dimming in 0–100\thinspace{}\% range with low colour shifts. \am{} on the other hand yielded higher efficacies at lower current levels but at the price of much higher colour shifts and limited dimming range.
  376. Work on more advanced dimming mechanisms started with Ashdown~\cite{Ashdown2006} proposing the use of pulse-code modulation to remove the need of hardware \pwm{} generators to dim the diodes.
  378. \marginpar{hybrid dimming}In 2009, Tse~et~al.~\cite{Tse2009b} from Hong Kong Polytechnic University proposed a general driving scheme of using two current levels to control \led{}s. Using low current magnitude together with high current \pwm{} signal both high efficiency and limited colour shifts can be achieved. A~variation of their general driving technique, where one of the currents is set to zero, was investigated in this dissertation.
  380. \marginpar{driving current}In 2007, Schmid~et~al.~\cite{Schmid2007} measured the effect of the typical driving currents on diode performance. The accuracy of the results is limited due to the measurement hardware used but it shows the effect of ripple content in the diode current on the light output.
  382. \subsection{Thermal properties}
  383. Thermal properties of high-power \led{}s have been a topic for many research programmes since mid 2000s. Scientists quantified the effect of junction temperature on such parameters as: luminous flux, colour point, spectrum shape and lifetime of the device. Three major research topics can be identified on the thermal properties of \led{}s: methods for measuring the junction temperature, estimating the junction temperature using a thermal model of \led{} luminaires and the research on the effect of temperature on the device properties.
  385. \noindent{}\marginpar{junction measurements}The measurement of the junction temperature is done mostly through  indirect methods, where the temperature dependent parameter is measured and using model or previous calibration the temperature is calculated.
  387. In 2004 Xi et~al.~\cite{Xi2004} proposed a method based on forward voltage measurement. In order to calibrate the measurement, the diode is measured when driven with very small duty cycle, typically 0.1\thinspace{}\%, to reduce the self heating effect. This method was used throughout the writing of this dissertation. Another method provided by Xi et~al.~\cite{Xi2005} is based on the relationship between the peak of the spectrum and the junction temperature.
  389. A method based on the diode spectrum was presented by Chhajed et~al.~\cite{Chhajed2005a}. It is based on the high energy slope of the diode's spectrum. When plotted in semilogarithmic scale, the slope is proportional to the carrier temperature.
  391. \marginpar{thermal modelling}
  392. Instead of measuring the junction temperature indirectly, a thermal model of the diode can be used to estimate it. Farkas~et~al.~\cite{Farkas2003} pioneered the use of diode thermal model. They applied a method described by Székely in 1988~\cite{Szekely1988}, used previously for high power silicon chips to measure and model the thermal path of the dissipated heat. This thermal model is used further in the dissertation to simulate the performance of both \led{}s and whole luminaires.
  394. \marginpar{temperature effects}
  395. The effect of junction temperature on InGaN and AlGaInP has been documented by Chhajed et~al.~\cite{Chhajed2005b} in 2005. They measure the effect of varying the junction temperature on colour diodes and on trichromatic luminaires composed of these diodes. \marginpar{spectrum model}Moreover, they described the \led{} spectrum model based on a Gaussian function. They used the model to find a good combination of primary diode colours to create a luminaire with high luminous efficacy and colour rendering index.
  397. In 2006, Man and Ashdown~\cite{Man2006} extended the model presented by Chhajed into double Gaussian model composed of sum of two Gaussian curves, accurately modelling the \led{} spectrum in wide range of temperatures.
  399. \marginpar{flux estimation}The linearity of \pwm{} dimming method is limited to a constant heatsink temperature case but this temperature changes with average driving current, bending the current-luminous flux relationship. Garcia~et~al.~\cite{Garcia2008} created a flux estimator basing on the \led{} case and ambient temperatures to deal with transient and steady-state behaviour of the luminous flux. This estimator can accurately predict the instantaneous value of luminous flux, therefore it allows a linear dimming behaviour and precise intensity control, but lack of the colour shift information limits the practical use to control of a phosphor-converter white \led{}s or a single colour diode. It is not sufficient for precise colour control of polychromatic luminaires.
  401. Hui and Qui~\cite{Hui2009} analysed the static performance of the \led{}-heatsink system in 2009. They showed that above certain thermal resistance in the heat flow path limiting the driving current will increase the maximum intensity of the diode. This phenomenon is very important for driving AlGaInP diodes which have very high temperature dependency.
  403. \noindent{}\marginpar{multi-domain models}Huang and Tang~\cite{Huang2006} presented a complete thermal-electrical-luminous model of an \led{} luminaire in 2009. Derivation of system coefficients was done using power perturbation on individual diodes and system identification methods. Unfortunately, the model disregards the colour shifts and nonlinearities present in \led{} based luminaires.
  405. \subsection{Lifetime}
  406. \marginpar{failure mechanisms}In 2004, Narendran et.~al.~\cite{Narendran2004} analysed the failure mechanisms in white \led{}s. Junction temperature was identified as major contributor in light output degradation phenomenon. Since then, many long time tests have been performed in order to create a reliable model of light output degradation. Jacob~et~al.~\cite{Jacob2006} analysed the failure modes of high power \led{}s in 2006. They mention a few separate failure mechanisms: shunting of the structure through parasitic dendrites growing on the structure, leakage paths in the active area or bond wire dislocations.
  408. \marginpar{gallium nitride diodes}In 2008, Hu~et~al.~\cite{Hu2008} performed a degradation test on GaN diodes. They report a change in thermal structure of the diodes. The increase in thermal resistance of the package increases the junction temperature and speeds up degradation.
  410. \marginpar{gallium phosphide diodes}Mathematical model of light output degradation for AlGaInP diodes was derived by Grillot et~al.~\cite{Grillot2006} in 2006. They proved that useful lifetime of the device is a linear function of current density and logarithmic function of stress time.
  412. \marginpar{driving current impact}Meneghini~et~al.~\cite{Meneghini2006} and Buso~et~al.~\cite{Buso2008} compared the \led{} output degradation speed under \dc{} and pulsed driving current. Their reliability tests showed that the device lifetime is similar under both driving conditions. However, different families of devices showed different effect on both driving conditions, indicating that the internal structure of the device and package may influence the reliability behaviour under various driving currents.
  414. \subsection{Luminaire control}
  415. Many different control techniques have been developed over the years to control light-emitting diodes based luminaires. Basic colour control loops have been presented by Deurenberg et al.~\cite{Deurenberg2005}.
  417. \marginpar{open loop}The simplest open loop control is very easy to implement as it requires only a single point calibration but it does not compensate for self heating of the diode structure, change of ambient temperature or the reduced intensity with time.
  419. \marginpar{colour feedback}Most common feedback mechanism implemented to control \led{} luminaires are the colour sensors. Subramanian et~al.~\cite{Subramanian2002} describe a simple colour control loop that reduces the colour shift of the luminaire compared to an open loop solution. The dynamics of this system is described later by Subramanian and Gaines~\cite{Subramanian2003}. Over the years many researchers used various colour sensors to control luminaires: Ackermann et~al.~\cite{Ackermann2006}, Lim et~al.~\cite{Lim2006} and Chang et~al.~\cite{Chang2009}. The consensus is that colour sensors provide easy feedback method with relatively good accuracy but \led{} spectra shifts can create colour errors. Additionally, the colour sensor stability is a major concern. To increase the accuracy of the colour sensors Robinson and Ashdown~\cite{Robinson2006} proposed an advanced feedback solution capable of estimating the spectral shifts.
  421. \marginpar{flux feedback}Subramanian et~al.~\cite{Subramanian2002} proposed using a single photodetector instead of colour sensor and measure each of the diodes separately using a timing scheme.
  423. \marginpar{temperature feedback}As temperature is the main cause of the variation of diodes' parameters, some control schemes are based on the temperature dependence models. Moisio et~al.~\cite{Moisio2005} measured the heatsink temperature and approximated the intensity of the \led{}s installed in the luminaire.
  424. Forward voltage of the diode is proportional to the junction temperature at a single current level. Qu et~al.~\cite{Qu2007} measured the diode voltage and converted it directly to tristimulus values.
  426. \marginpar{model based control}Huang et~al.~\cite{Huang2009b} described the system model of a multi-chip polychromatic luminaire. They used power perturbation and system identification methods for each colour string of diodes.
  428. \subsection{Optimisations}
  429. \marginpar{primaries selection}Luminaires can be optimised from a few different angles. The number of primary colours and the exact device choice to be installed in the luminaire can be subjected to optimisation procedures. Žukauskas et~al.~\cite{Zukauskas2001} analysed the effect of number of primaries in the luminaire on the colour quality and luminous efficacy achievable by the luminaire. Chhajed et~al.~\cite{Chhajed2005b} used a Gaussian \led{} spectrum model to choose the best combination of three colours yielding high \cri{} value in the required colour range.
  431. \marginpar{driving}When the luminaire consists of more than three primary colours, the colour control system becomes underdetermined and diode control can be optimised. Ries et~al.~\cite{Ries2004} showed the possibility to optimise such parameters of luminaires as luminous flux, efficacy and colour rendering index. Ou-Yang and Huang~\cite{Ou-Yang2007} showed a method to estimate the gamut of multi-primary \led{} based luminaires.
  433. In 2010, Lin~\cite{Lin2010} presented a method to optimise the \cri{} of a luminaire using mathematical formulations of the problem and numerical method based on the complex method.
  435. \subsection{\led{} driver}
  436. The work of Schmid~et~al.~\cite{Schmid2007} on the impact of driving current shape and various research on different dimming techniques indicates that in order to achieve high efficacy the diode should be driven with a \dc{} current. Van der Broeck et~al.~\cite{Broeck2007} analysed different converter topologies for driving \led{}s from a \dc{} voltage source using both \dc{} and pulsating current. They present different types of isolated and non-isolated converters describing their output current shape. Their conclusion is that pulsed output current converters contain, in general, less components (like output capacitor) and therefore may be beneficial for the lifetime of the converter.
  438. Many different topologies were applied for driving high current \led{}s. Depending on the application, solutions include using a buck by Torres et~al.~\cite{Torres2007}, boost by Xiaoru and Xiaobo \cite{Xiaoru2008}, Ćuk by de~Britto et~al.~\cite{deBritto2008}, flyback by Pan et~al.~\cite{Pan2007} or \textsc{sepic} by Zhongming et~al.~\cite{Zhongming2008} etc. Most of the research on \led{} drivers focus on high-power \led{}s supplied with a current up to 700\thinspace{}mA. Although some of the solutions presented for these diodes can be easily scaled to a diode driven with 13.5\thinspace{}\textsc{a}, the high current converter design calls for more advanced solution.
  440. Early processors were supplied with a few volts and tens of amps, similarly to high power \led{}s. For this kind of load Xunwei~et~al.~\cite{Xunwei1999} proposed using a multiphase interleaved buck converters.
  442. \section[Scientific contributions]{scientific contributions}
  443. Following contributions of this work have been identified as extending the current state of the art.
  445. \marginpar{hybrid pwm/am dimming}Introduction of hybrid \pwm{}/\am{} dimming technique (chapter~\ref{ssec:hybrid_pwm_am}). This dimming technique is a combination of pulse width modulation with variable peak current amplitude. It can be easily adapted to current designs as integrated \led{} drivers typically have the possibility to use both \pwm{} and \am{} dimming.
  447. Typically, one variable quantity is used to control the intensity of the diode: either a duty cycle or \dc{} current. By controlling a single \led{} with both peak current and duty cycle an additional degree of freedom appears and \led{} parameters can be optimised.
  449. Hybrid dimming mechanism can be used e.g. to reduce the colour point shifts of a white, phosphor converted diode under varying heatsink temperature or to minimise spectra shifts in InGaN diodes.
  451. Application of hybrid dimming mechanism to polychromatic luminaires creates opportunities to optimise luminaires' control, increase their efficacy and device gamut (chapter~\ref{ssec:Optimal_control_hybrid_dimming}). Using this technique, optimisation possibilities are given even for di- and trichromatic luminaires.
  453. \marginpar{current voltage led model}\led{} empirical model, described in chapter~\ref{sec:iv_model}, based on instantaneous values of current and voltage of the diode can form a basis of a luminaire colour control scheme (chapter~\ref{ssec:iv2xyzColourControl}). The model is valid for any form of diode dimming and accurately predicts colour shifts created by the changes of ambient temperature and by the dimming technique. The model can be used as a feedback mechanism for precise colour control.
  455. \noindent{}\marginpar{application of interleaved converters}An interleaved buck topology was proposed for driving high current \led{} light sources. Because of relatively low dynamic resistance the working point of the driver can be set so that the converter always works close to 50\thinspace{}\% duty cycle where the ripple cancellation effect is the strongest. This removes the need for high value, electrolytic output capacitors.
  457. \marginpar{publications}
  458. Some of the work presented in this dissertation has also been published in other sources: \\
  460. \reversemarginpar
  461. \noindent{}\marginpar{1}S.\thinspace{}Bęczkowski. Advanced dimming strategy for solid state luminaires. In \textit{Proceedings of the 10th International Conference on Solid State Lighting}.  pp. 1–6. \textsc{spie} International Society for Optical Engineering. 2010.\\
  463. \noindent{}\marginpar{2}S.\thinspace{}Bęczkowski, S.\thinspace{}Munk-Nielsen. \textsc{led} spectral and power characteristics under hybrid \textsc{pwm/am} dimming strategy. In \textit{Proceedings of the \textsc{ieee} Energy Conversion Congress and Exposition, \textsc{ecce} 2010}, pp. 731–735. \textsc{ieee} press. 2010.\\
  465. \noindent{}\marginpar{3}S.\thinspace{}Bęczkowski, S.\thinspace{}Munk-Nielsen. Dual interleaved buck converter for driving high power \textsc{led}s. In \textit{Proceedings of the 2011—14th European Conference on Power Electronics and Applications, \textsc{epe} 2011}, pp.1–6. 2011. \\
  467. \noindent{}\marginpar{4}L.\thinspace{}Török, S.\thinspace{}Bęczkowski, S.\thinspace{}Munk-Nielsen, J.\thinspace{}Gadegaard, T.\thinspace{}Kari, K.\thinspace{}Pedersen. High output \textsc{led}-based profile lighting fixture. (accepted for publication) \textit{Proceedings of \textsc{iecon 2011}—37th Annual Conference of the \textsc{ieee} Industrial Electronics Society, \textsc{iecon 2011}}, 2011.\\
  469. \noindent{}\marginpar{5}\textsc{dk\thinspace{}pa\thinspace{}2011\thinspace{}70529} patent application. S.\thinspace{}Bęczkowski. Method of controlling illumination device based on current-voltage model.
  471. \normalmarginpar
  472. \section[Methodology]{methodology}
  473. The manufacturing of light-emitting diodes is a complex process performed in clean room environment. Diode structures are created using epitaxial growth and the diode colour is controlled by the GaN/InN ratio or by the AlInP composition~\cite{Krames2007}. On account of local variations, colour point of diodes in a batch can differ from sample to sample. The same applies to intensity and forward voltage magnitude. Manufacturers split the diode batches into smaller bins with similar diode's parameters like forward voltage, intensity and colour point.
  475. Data shown in the dissertation has been obtained using measurements on particular diode sample. No statistical analysis has been done, therefore the results apply to the measured diodes and do not apply to the diode batch in general. However, the principles demonstrated on particular devices are applicable to the whole batch but may yield different results. Different internal diode structures may show different behaviour. All tests were performed on high power, state of the art devices.
  477. Tests performed at fixed temperatures are performed using constant heatsink temperature. The use of constant junction temperature eliminates the impact of thermal structure of the diode on results, but its measurement a complex process, therefore it was not used in the data gathering. Furthermore, the experiments can be easily recreated using inexpensive hardware.
  479. \subsection{Measurement hardware}\label{measurementHardware}
  480. \marginpar{spectrometer}Optical measurements were performed using an Instrument Systems \textsc{cas}~140~\textsc{ct} array spectrometer. This spectrometer uses a diffraction grating technique for measurement of the visible part of the light spectrum. Light is directed on an optical component with a periodic structure which diffracts the light. Diffracted light is analysed by a \textsc{ccd} detector. The spectrometer can analyse light within 360–830\thinspace{}nm range with 2.2\thinspace{}nm resolution and 0.5\thinspace{}nm data point interval.
  482. \begin{figure}[h]% fig:integrating sphere
  483.         \centering
  484.         \includegraphics{graphics/Introduction/MeasurementHardware/IntegratingSphere.pdf}
  485.         \caption{Instrument Systems \textsc{isp}250 integrating sphere used in the test setup (left) and the principles of integrating sphere operation (right). Input beam is reflected diffusely inside the sphere many times before hitting the output port. A baffle  inside the sphere blocks the input light from hitting the output port directly.}
  486.         \label{fig:integrating_sphere}
  487. \end{figure}
  489. \noindent{}\marginpar{integrating sphere}A 25\thinspace{}cm diameter Instrument Systems \textsc{isp}250 integrating sphere shown in figure~\ref{fig:integrating_sphere} is used to gather the light from device under test. Inner surface of the sphere is coated with diffusely reflecting material, barium sulphate ($\mathrm{BaSO}_{4}$), of a well defined reflectance. Light is diffusely reflected from the walls many times before it hits the output port of the sphere, therefore the measurement is not dependent on the shape of the input light beam. An Instrument Systems \textsc{isp}500-220 adapter plate is used to mount a \led{} test fixture.
  491. \marginpar{test fixture}Instrument Systems \textsc{led}-850 test fixture is used for \led{} assembly. It is possible to test various \led{} types from different manufacturers with special adapters. These adapters provide electrical contacts for four wire sensing and a thermal connection between \led{} and thermally controlled plate. Arroyo Instrument 5310 \textsc{tecs}ource temperature controller was used to control the temperature of the plate.
  493. Measurements of power converter were done using following hardware: Fluke 8845~\textsc{a} multimeter, 100\thinspace{}MHz 4~channel Tektronix \textsc{tds}~3014\textsc{c} digital phosphor oscilloscope, 50\thinspace{}MHz Tektronix \textsc{tcp}202 current probe, Delta Elektronika \textsc{sm}~52-\textsc{ar}-60 power supply and Zentro-Elektrik \textsc{el}1000 electronic load.
  495. \textsc{ni}~6215 data acquisition module containing 16 analog inputs (16~bit, 250\thinspace{}kS/s) was also used for bigger tests that needed high number of input channels.
  497. \begin{table}[t]\footnotesize
  498.         \caption{National Instruments compact\textsc{daq} modules used in the experiments.}
  499.         \label{tab:NIcDAQmodules}
  500.         \centering
  501.         \begin{tabular}{cccll}
  502.                 \textsc{module} & \textsc{resolution} & \textsc{sampling rate} & \textsc{range} & \textsc{comments} \\
  503.                 \hline
  504.                 \abbr{ni 9201} &        12 bit  &       500\thinspace{}kS/s & ±10\thinspace{}V                                         & 8-channel  \\
  505.                 \abbr{ni 9229} &        24 bit  &       \phantom{0}50\thinspace{}kS/s  & ±60\thinspace{}V                                              & 4-channel  (isolated) \\
  506.                 \abbr{ni 9227}  &       24 bit  &       \phantom{0}50\thinspace{}kS/s  & \phantom{±6}5\thinspace{}A \textsc{rms} & 4-channel \\
  507.                 \hline
  508.         \end{tabular}
  509. \end{table}
  510. For electrical measurements a National Instruments compact\textsc{daq} system was used. The system consists of a \textsc{ni}~c\textsc{daq}-9172 chassis with interchangeable measurement modules (table~\ref{tab:NIcDAQmodules}). Voltage output modules were used to control \led{} drivers. Current and voltage inputs were used to acquire electrical quantities.
  512. \section[Outline]{outline}
  513. \marginpar{chapter 1}The first chapter gives background to the topic of this dissertation, followed by review of previous scientific work done by other parties to show the state of the art in light-emitting diodes control. Finally methodology used in the work is described. %in detail.
  515. \marginpar{chapter 2}Colour theory chapter gives an insight of how the light is perceived by the human eye and how it is measured. Various colour spaces, used throughout the dissertation, are defined. The definitions of light quality and colour distance are given in order to quantify the performance of the colour control.
  517. \marginpar{chapter 3}The principles of operation and properties of light-emitting diodes are described in the third chapter. Firstly the physical structures of the diodes are discussed.  Optical and electrical properties are analysed with focus on the parameter change with respect to junction temperature variations. The lifetime of the devices is analysed. The impact of dimming the diodes with different methods is presented. Hybrid dimming technique is introduced and its properties in increasing the controllability of \led{}s are analysed. In the end, an empirical model of the light-emitting diodes is presented and verified.
  519. \marginpar{chapter 4}Luminaire control is presented in the fourth chapter. Firstly the review of state of the art control methods is given. A colour control scheme based on the empirical model, described in chapter three, is presented. Finally optimisation routines for designing and controlling polychromatic luminaires are described.
  521. \noindent{}\marginpar{chapter 5}Fifth chapter reviews available power converter topologies with respect to driving light-emitting diodes. An interleaved buck topology is analysed and its benefits and drawbacks are discussed. Power stage and control are designed. Laboratory results are presented and discussed.
  523. \marginpar{chapter 6}Final sixth chapter summarises the dissertation by providing an overview of the results and gives recommendations on the future work.
  525. % chapter introduction (end)
  527. %%%%%%%%%%%%%%%
  528. %% Chapter 2 %%
  529. %%%%%%%%%%%%%%%
  530. \chapter{Colour theory}% (fold)
  531. \label{cha:colour_theory}
  532. Colour science has its origins in late 17th century when Sir Isaac Newton conducted his experiments on the nature of light. Using two prisms Newton showed that white light can be divided into separate colours and that pure colours cannot be divided more. In order to prove that the prism is not colouring the light, he reconstructed white light using pure colours.
  533. \begin{figure}[h]% fig:newton_prism_diagram
  534.         \centering
  535.         \includegraphics{graphics/ColourTheory/NewtonPrismDiagram.pdf}
  536.         \caption{The diagram from Sir Isaac Newton's experiment. The first prism (\textsc{a}) splits the white light into its constituent colours. A single colour part of the spectrum hits the second prism (\textsc{f}) proving that a pure colour cannot be divided any more.\label{fig:newton_prism_diagram}}
  537. \end{figure}
  539. \noindent{}Colour is sensation that is dependent on both light and its properties and observer. In order to quantify it physical properties of light and observers have to be known.
  541. \section[Human vision]{human vision}
  542. \marginpar{human eye}The eye is an organ used for vision. Frontal part on an eye consisting of iris, pupil and lens acts similarly to the camera lens. Iris limits the amount of light going inside the eye thus protecting it from overexposure and damage. Lens focuses the image on the back, inner surface of an eye—retina. Chemical reactions in rods and cones in the retina produce electrical stimuli for the brain.
  544. \begin{figure}[t]% fig:eye_cross_section
  545.         \centering
  546.         \includegraphics{graphics/ColourTheory/eye.pdf}
  547.         \caption{Eye cross-section (left) and scanning electron microscope image of human retina (right, source: Science Photo Library, image \textsc{f001/0041}).\label{fig:eye_cross_section}}
  548. \end{figure}
  550. Rods and cones are photoreceptor cells in the retina responsible for low and high intensity light vision, respectively. Human eye uses photopic vision when luminance is above approximately 3\thinspace{}cd/m² and scotopic vision with luminance below approx. 0.03\thinspace{}cd/m². In between, the eye uses mesopic vision, a combination of the above. \marginpar{rods}Rods are sensitive enough to detect a single photon, but they are sensitive only in 400–640\thinspace{}nm wavelength range. Their maximum sensitivity is 1700\thinspace{}lm/W at 507\thinspace{}nm. Because only one type of rod cells exist, under photopic vision regime, there is no colour sensation. \marginpar{cones}Cone cells, on the contrary, are much less sensitive to the light, but their three variants (short, medium and long wavelength) shown in fig.~\ref{fig:Cones_Sensitivity} allow us to perceive colour through the opponent process of colour vision.
  552. \begin{figure}[ht]% fig:eye_cross_section
  553.         \centering
  554.         \includegraphics{graphics/ColourTheory/ConeSensitivity.pdf}
  555.         \caption{Sensitivity of cones after Moses and Hart \cite{PhysiologyOfTheEye}. \label{fig:Cones_Sensitivity}}
  556. \end{figure}
  558. \noindent{}The transition between scotopic and photopic vision is a very slow process. The speed depends on the direction of the transition and initial adaptation level. Full adaptation from light to dark conditions may take up to 45~minutes. Reverse transition is faster, but still may take up to few minutes. All calculations in this dissertation assume fully adapted photopic vision regime.
  560. \marginpar{eye sensitivity}Conversion between radiometric and photometric units is done using the \emph{eye sensitivity function}~$V\textrm{(λ)}$, also called \emph{luminous efficiency function}, shown in fig.~\ref{eye_sensitivity}. The function represents relative eye sensitivity to different wavelengths under photopic vision regime. The maximum eye sensitivity is in green spectral range at 555\thinspace{}nm wavelength. This point corresponds to the sensitivity of 683.002\thinspace{}lm/\textsc{w}.
  561. \begin{figure}[t]% fig:eye_sensitivity
  562.         \centering
  563.         \includegraphics{graphics/ColourTheory/EyeSensitivityFunctions.pdf}
  564.         \caption[Eye sensitivity function]{Eye sensitivity function \textit{V(λ)}: \abbr{1931} standard, \abbr{1978} Vos modification and most recent shape proposed in \abbr{2005} by Sharpe, Stockman, Jagla and Jägle.\label{eye_sensitivity}}
  565. \end{figure}
  566. Eye sensitivity function was obtained experimentally by the \emph{minimum flicker method}. The stimulus in the experiment was a circular area alternately illuminated by two different colours with 15\thinspace{}Hz frequency. At this frequency the two hues fuse into one colour however the brightnesses does not. By adjusting one light source properties, human subject was to minimise the visible flicker.  
  568. There have been several attempts to improve the original \abbr{1931} standard function: Judd in \abbr{1951}, Vos in \abbr{1978}~\cite{Vos1978} and recently, in~\abbr{2005}, Sharpe, Stockman, Jagla and Jägle~\cite{Sharpe2005} proposed a modified \textit{V*\textrm{(λ)}} function. The blue region of eye sensitivity has been underestimated in the \abbr{cie 1931} standard, but this sensitivity function is still commonly used as a basic observer. All measurements and calculations presented in this dissertation will use \abbr{cie 1931} standard unless explicitly mentioned otherwise.
  570. \section[Photometric quantities]{photometric quantities}
  571. Light can be characterised in two different ways. Radiometric quantities describe physical properties of an electromagnetic wave. Photometric units, on the other hand, characterise light as perceived by human eye. Radiometric measurements are very straightforward, but not useful if one wants to define the sensation that is caused by the light. For example, infra-red and ultra-violet light can be defined in terms of radiometric quantities, but is invisible to the human eye.
  573. \marginpar{luminous flux} Luminous flux is a measure of a perceived power of the light. The unit of luminous flux is \emph{lumen}~(lm). A monochromatic light source emitting an radiometric power of 1/683~watt at 555\thinspace{}nm has a luminous flux of 1\thinspace{}lm.
  575. \marginpar{luminous intensity} Luminous intensity is a measure of a perceived intensity of light per unit solid angle. The unit of luminous intensity is \emph{candela}~(cd), which is a base \textsc{si}~unit. The current definition is: a monochromatic light source emitting an optical power of 1/683\thinspace{}\textsc{w} at 555\thinspace{}nm into the solid angle of 1~steradian~(sr). In other words, a light source with luminous flux of 1\thinspace{}lm emitting into 1\thinspace{}sr has a luminous intensity of 1\thinspace{}cd. An isotropically emitting light source (emitting into 4π~sr) with the same luminous flux would have a luminous intensity of 1/4π~cd. The luminous intensity parameter is therefore closely related to light source output geometry. Any lens or other optical aid can greatly influence this parameter.
  577. \marginpar{luminance} Luminance is a measure of the luminous intensity of light traveling in a given direction per unit area. It describes the amount of light that is emitted from a particular area within a given solid angle. Luminance is typically used to characterise light sources with big emitting surfaces like \textsc{lcd} screens or \led{}s with active area greater than $1 \mathrm{mm}^{2}$.
  579. \marginpar{illuminance} The measure of the total luminous flux incident on a given area is called illuminance. A luminous emittance, similarly, is the total luminous flux emitted from a given area. Both are measured in lux~(lx) or lumen per square meter~($\mathrm{lm/m^{2}}$).
  581. The constant 1/683 is a legacy constant, used to equalise the current definitions of light intensity with the old ones. The ``standard candle'' (pure spermaceti candle weighing one sixth of a pound and burning at a rate of 120 grains per hour), used as a reference source in England prior to \abbr{1948}, had roughly 1/683 watts of power.
  583. \section[Colour spaces]{colour spaces}
  584. A colour space is a method of quantifying the sensation of colour. Humans describe the colour sensation by means of brightness, colourfulness and hue. During printing, colour is made by matching cyan, magenta, yellow and black inks (\textsc{cmyk}) on paper. Computer screen or digital photo camera describes the colour by the amounts of red, green and blue. Note that the reds, greens and blues of these devices, are not generally the same. To math the colours of the camera on the screen, a consistent \srgb{} colour space is used in both devices. Many different colour spaces exist: some are intuitive to use by humans, some are used when working with particular hardware (device dependent). Other may serve particular need, e.g.~perceptually uniform colour spaces, but they all describe the colour sensation in a quantitative way.  
  586. \marginpar{tristimulus theory of colour}Light sources can have a very complex spectral distributions, still a single colour can be described using only three scalar stimuli parameters. This phenomenon was first postulated in 19th century by Thomas Young and then developed further by Herrman von~Helmoltz. The tristimulus theory was experimentally verified in \abbr{1930} by the colour matching experiments. The human subjects were to match the colour of a monochromatic light to the light created by mixing three monochromatic lights (primaries). The \abbr{1930} experiments used light sources at 700\thinspace{}nm~(\textsc{r}), 546.1\thinspace{}nm~(\textsc{g}) and 438.1\thinspace{}nm~(\textsc{b}). Unfortunately, the colour primaries selection caused some of the colour matching functions' values to be negative. As negative values of light have no physical meaning instead of subtracting one of the tristimulus colours from the matching light, the colour was added to the sample monochromatic light. %Only then both colours could be matched.
  588. \marginpar{cie~\abbr{1931} \xyz{} colour space}
  589. In \abbr{1931} the \cie{} defined new primaries to overcome the problem with awkward negative weights in colour matching function. New $\overline{x}$, $\overline{y}$ and $\overline{z}$ functions had only positive values and $\overline{y}$ has been chosen to match eye sensitivity function $V(λ)$.
  591. Due to the nature of spatial distribution of cones inside the eye, the perception of a colours depends on subject's field of view. To eliminate the influence of this parameter, \cie{} introduced a 2° standard colorimetric observer. A 10° standard observer has been introduced in \abbr{1964} as a result of the work of Stiles and~Burch~\cite{Stiles1959} recommended for field of view more than 2°.
  592. \begin{figure}[t]% CIE_x_y_z_
  593.                 \centering
  594. \includegraphics{graphics/ColourTheory/ColourMatchingFunctions.pdf}
  595. \caption{\cie{} \abbr{1931} 2° and \abbr{1964} 10° observer $\overline{x}$, $\overline{y}$ and $\overline{z}$ colour matching functions, normalised to the maximum of \abbr{1931}, 2° observer $\overline{x}$ colour matching function.\label{fig:CIE_x_y_z_}}
  596. \end{figure}
  598. \marginpar{tristimulus values}In order to calculate tristimulus values \textit{X, Y} and \textit{Z} the spectrum of light source needs to be multiplied by respective colour matching functions $\overline{x}$, $\overline{y}$ and $\overline{z}$ and integrated over wavelength.
  599. \begin{equation} \label{eq:tristimulus}
  600.         X=\!\int_{380}^{780}\!\!\!\!\!\!\!P(λ)\overline{x}(λ)\ud{}λ\qquad{}
  601.         Y=\!\int_{380}^{780}\!\!\!\!\!\!\!P(λ)\overline{y}(λ)\ud{}λ\qquad{}
  602.         Z=\!\int_{380}^{780}\!\!\!\!\!\!\!P(λ)\overline{z}(λ)\ud{}λ
  603. \end{equation}
  604. Because the $\overline{y}$ matches the eye sensitivity function the $Y$ can be used as a measure of luminance of a light. Multiplying all tristimulus values by 683\thinspace{}lm/\textsc{w} gives the values in lumens therefore \textit{Y} becomes equal to the luminous flux. It is common to describe the colour of the light using \abbr{cie} $xyY$ colour space. Chromaticity coordinates, that can be represented on 2~dimensional $xy$-plane (fig.~\ref{fig:CIE_xy}) are obtained by normalising the tristimulus components.
  605. \begin{equation} \label{eq:XYZ2xyz}
  606.         x=\frac{X}{X+Y+Z}\qquad{}
  607.         y=\frac{Y}{X+Y+Z}\qquad{}
  608.         z=\frac{Z}{X+Y+Z}
  609. \end{equation}
  610. Because of this normalisation, only two values are needed to represent any possible colour.
  611. \begin{equation} \label{eq:xy_only}
  612.         x+y+z=1\qquad{}\Longrightarrow{}\qquad{}z=1-x-y
  613. \end{equation}
  614. Intensity of the light is not shown on the $xy$-plane. Any two light sources with the same hues and saturations but different lightness will project on the same point of a diagram. Using $x$, $y$ and the stimulus $Y$ one can easily calculate back the values of $X$ and $Z$ stimuli.
  615. \begin{equation}  \label{eq:xyY2XYZ}
  616.  X = \frac{Y}{y}x \qquad{}
  617.  Y =Y \qquad{}
  618.  Z = \frac{Y}{y}(1-x-y)
  619. \end{equation}
  621. \noindent{}When two lights with different spectral power distributions have the same colour chromaticity (same $xy$ coordinates) they are called metamers and the property itself is called metamerism\marginpar{metamerism}. This allows, for example, to mix a white light using blue and yellow (as is in white phosphor converted \led{}s) or red, green and blue. This can also lead to big differences between colour rendering of two lights with the same colour chromaticity. For example, white, phosphor based \led{} does not reproduce well saturated red colours, while white, trichromatic \rgb{} \led{} can have an excellent red colour rendering.
  622. \begin{figure}[t]% CIE_xy
  623.                 \centering
  624. \includegraphics{graphics/ColourTheory/CIExy.pdf}
  625. \caption{\cie{} \abbr{1931} $xy$ chromaticity diagram.\label{fig:CIE_xy}}
  626. \end{figure}
  628. \marginpar{xy chromaticity diagram}The $xy$-plane is a projection of a linear colour space and therefore colours can be linearly mixed on it. Therefore, every polychromatic source will have a gamut (available colour range) that is a convex hull on all primary colour points inside the horseshoe diagram.
  630. \noindent{}Monochromatic light sources (pure colours) lie on the curved boundary of the horseshoe. Straight line between blues and reds is called \textit{purple line}, colours laying on this line cannot be made using a monochromatic light source, they are created by mixing saturated red and blue colours.
  632. \emph{Equal energy point} or \emph{white point} (x=1/3, y=1/3) corresponds to the white light where all tristimulus values are equal. Any two colours that can be connected by a straight line crossing the white point are called complementary.
  634. \marginpar{perceptual uniformity}This colour space is not perceptually uniform. This means that the difference between colours of two points on the plane does not correspond to the geometrical distance between them. In \abbr{1942}, MacAdam analysed colour differences of closely spaced points~\cite{MacAdam1942}. He noted, that depending on a position in the $xy$-plane, different geometrical distance between two colours yields a noticeable colour difference. Similar colours, that appear identical to human eye, can be grouped on a $xy$-diagram in ellipsis shaped areas (fig.~\ref{fig:CIE_xy_MacAdam}). Ellipses in the green regions are very big compared to the ellipses in the blue region. This non-uniformity initiated the search for uniform chromaticity diagram. In result, in \abbr{1960}, the \cie{} introduced $uv$, and in \abbr{1976}, $u'v'$ uniform chromaticity diagram~\cite{Wyszecki2000}.
  636. \begin{figure}[t]% CIE_uv
  637.         \centering
  638. \includegraphics{graphics/ColourTheory/CIEuv.pdf}
  639. \caption{\cie{} \abbr{1976} $u'v'$ chromaticity diagram.\label{fig:CIE_uv}}
  640. \end{figure}
  641. \marginpar{cie~\abbr{1960} $Luv$ and cie~\abbr{1976} $L^{*}u^{*}v^{*}$ colour spaces}Coordinates for $uv$ and $u'v'$ can be calculated from tristimulus values
  642. \begin{equation} \label{eq:uv}
  643.         u=\frac{4X}{X+15Y+3Z}\qquad{}
  644.         v=\frac{6Y}{X+15Y+3Z}
  645. \end{equation}
  646. and
  647. \begin{equation} \label{eq:u'v'}
  648.         u'=\frac{4X}{X+15Y+3Z}\qquad{}
  649.         v'=\frac{9Y}{X+15Y+3Z}.
  650. \end{equation}
  651. These coordinates can also be calculated directly from $xy$ coordinates using equations:
  652. \begin{equation} \label{eq:x2u}
  653.         u=u'=\frac{4x}{-2x+12y+3}
  654. \end{equation}
  655. and
  656. \begin{equation} \label{eq:y2v}
  657.         v=\frac{6y}{-2x+12y+3}\qquad{}
  658.         v'=\frac{9y}{-2x+12y+3}.
  659. \end{equation}
  660. Reverse transformations are also possible, using:
  661. \begin{equation} \label{eq:uv2xy}
  662.         x=\frac{3u}{2u-8v+4}\qquad{}
  663.         y=\frac{2v}{2u-8v+4}
  664. \end{equation}
  665. and
  666. \begin{equation} \label{eq:u'v'2xy}
  667.         x=\frac{9u'}{6u'-16v'+12}\qquad{}
  668.         y=\frac{2v'}{3u'-8v+6}.
  669. \end{equation}
  671. \noindent{}\marginpar{colour space conversions}A conversions from one colour space to another requires a choice of a different set of primary stimuli. The only limitation is that the vectors created by the new set of stimuli $R'G'B'$ is linearly independent. Therefore no combination of $r$, $g$ and $b$ should make
  672. \begin{equation}
  673.         rR' + gG' + bB' = 0,
  674. \end{equation}
  675. except for trivial $r=g=b=0$. Each primary stimulus of the new set can be matched using the mixture of the old primaries therefore we can write relation between $RGB$ and $R'G'B'$ colour space as:
  676. \begin{equation} \label{eq:rgb2r'g'b'}
  677.         \begin{array}{r@{\:}c@{\:}l}
  678.         R'&=&a_{11}R + a_{12}G + a_{13}B \\
  679.         G'&=&a_{21}R + a_{22}G + a_{23}B \\
  680.         B'&=&a_{31}R + a_{32}G + a_{33}B
  681.         \end{array}
  682. \end{equation}
  683. A transformation matrix $\mathbi{A}$ is formed by the transformation coefficients $a_{ij}$
  684. \begin{equation} \label{eq:transformationMatrix}
  685. \mathbi{A} =
  686. \left[ \begin{array}{ccc}
  687. a_{11} & a_{12} & a_{13} \\
  688. a_{21} & a_{22} & a_{23} \\
  689. a_{31} & a_{32} & a_{33}
  690. \end{array} \right].
  691. \end{equation}
  692. Transition from $R'G'B'$ to $RGB$ based colour space requires inverted matrix $\mathbi{A}^{-1}$.
  694. %% Black-body radiation
  695. \section[Black body radiation]{black body radiation}
  696. \noindent{}\marginpar{black body definition}A black body is an idealised physical object that absorbs all electromagnetic radiation that falls on it. The ideal object would appear perfectly black because all the visible light is absorbed. Any light emitted from the object would only be the function of object's temperature.
  698. \noindent{}\marginpar{planckian radiation spectrum}A planckian black body radiation spectrum is a very useful standard for defining white light. It is characterised using only single variable: \emph{colour temperature}. The spectrum is given by equation
  699. \begin{equation} \label{eq:black-body_spectrum}
  700.         I(λ)=\frac{2hc^{2}}{λ^5\Big[\mathrm{exp}\Big(\frac{hc}{λkT}\Big)-1\Big]}
  701. \end{equation}
  702. derived by Max~Planck in \abbr{1900}. The ideal black body object heated up to 6500\thinspace{}\textsc{k} would emit white light with colour temperature of 6500\thinspace{}\textsc{k}. An incandescent light bulb is an example of non ideal black body source. The temperature of the filament has to be kept below melting point and therefore the maximum colour temperature of this kind of light source is around 3500\thinspace{}\textsc{k}. Low luminous efficacy of the classical light bulb approximately 15\thinspace{}lm/\textsc{w} is the effect of the black-body spectrum having its maximum in the infrared part of electromagnetic radiation.
  704. \begin{figure}[t]% Blackbody Locus on CIExy and CIEu'v' diagrams
  705.                 \centering
  706. \includegraphics{graphics/ColourTheory/BlackbodyLocus.pdf}
  707. \caption{Black body colour locus on \cie{} \abbr{1931} $xy$ (left) and \cie{} \abbr{1976} $u'v'$ (right) diagram. Dotted lines show the position of correlated colour temperature points. Note that \cct{} lines on $u'v'$ diagram are perpendicular to black body locus. \label{fig:BlackbodyLocus}}
  708. \end{figure}
  709. \marginpar{black body locus}The colour point locus of black body radiator starts in the red corner of \cie{} \abbr{1931} diagram, then moves through the orange and yellow to end in the white region. This corresponds to the colours of a real heated objects. Colour temperatures between 2000\thinspace{}\textsc{k} and 4000\thinspace{}\textsc{k} are commonly referred to as warm white. Colour temperatures above 7000\thinspace{}\textsc{k} are called cool white.
  711. If the colour coordinates of a light source is not on the black-body locus but in its proximity then the \emph{correlated colour temperature} parameter is used. It is defined as temperature of a black body whose colour is the closest to the source colour. On a $u'v'$ diagram, the \cct{} can be determined as a colour temperature of a closest geometrically point on a planckian curve. Correlated colour temperature cannot be determined geometrically on a $xy$ plane, because the plane is non-uniform.
  713. To calculate chromaticity coordinates an cubic interpolation can be used~\cite{US7024034}. The $x$ coordinate is calculated for 1666–4000\thinspace{}\textsc{k} and 4000–25000\thinspace{}\textsc{k} temperature ranges respectively as:
  714. \begin{align}
  715. x &= -0.2661239\cdot{}10^{9}/T^{3} - 0.2343580\cdot{}10^{6}/T^{2} + 0.8776956\cdot{}10^{3}/T + 0.179910 \nonumber \\
  716. x &= -3.0258469\cdot{}10^{9}/T^{3} + 2.1070379\cdot{}10^{6}/T^{2} + 0.2226347\cdot{}10^{3}/T + 0.24039
  717. \end{align}
  718. The corresponding $y$ coordinate is calculated for 1666–2222\thinspace{}\textsc{k}, 2222–4000\thinspace{}\textsc{k} and 4000–25000\thinspace{}\textsc{k} temperature ranges respectively as:
  719. \begin{align}
  720. y &= -1.1063814 x^3 - 1.34811020 x^2 + 2.18555832 x - 0.20219683 \nonumber \\
  721. y &= -0.9549476 x^3 - 1.37418593 x^2 + 2.09137015 x - 0.16748867 \\
  722. y &= +3.0817580 x^3 - 5.87338670 x^2 + 3.75112997 x - 0.37001483 \nonumber
  723. \end{align}
  724. %1667\text{K} \leq T \leq 4000\text{K}
  725. %4000\text{K} \leq T \leq 25000\text{K}
  726. %1667\text{K} \leq T \leq 2222\text{K}
  727. %2222\text{K} \leq T \leq 4000\text{K}
  728. %4000\text{K} \leq T \leq 25000\text{K}
  731. %% Illuminants
  732. \section[Standard illuminants]{standard illuminants}
  733. \cie{} has standardized several light illuminants: \textsc{a}, \textsc{b}, \textsc{c}, \textsc{d} series, \textsc{e} and \textsc{f} series. Examples of their chromaticities and colour temperatures can be seen in table \ref{tab:standard_illuminants}. The purpose of standardizing the illuminants is to be able to compare the colour sensation of various objects under typical light sources. The sources include the sun spectra at different times of a day, the direct sunlight, the average incandescent bulb and fluorescent lamps of various composition.
  734. \input{tables/illuminants.tex}
  736. \noindent{}Illuminants from \textsc{d} series were introduced in 1967 to replace the \textsc{b} and \textsc{c} illuminants acting as daylight simulators~\cite{Judd1964}. Judd et~al. analyzed typical daylight spectra and tabulated the results with 10\thinspace{}nm increments. Obtaining higher resolution requires interpolating this data. The \textsc{d} series illuminants are used in calculation of the colour rendering index—the measure of the quality of light.
  738. %% Quality of light
  739. \section[Quality of light]{quality of light}
  740. Although many light sources with different spectra can have the same perceived colour, the appearance of various colour objects under those illuminants can vary significantly.
  742. \marginpar{colour rendering index}In order to easily compare the ability to reproduce colours a \cri{} (colour rendering index) has been developed~\cite{CIE13.3}. Colour rendering is defined as a \emph{measure of the degree of colour shift objects undergo when illuminated by the light source as compared with the colour those same objects when illuminated by a reference source of comparable colour temperature}. This means we compare a test light source to natural light or an ideal source e.g. black-body radiator with the same colour temperature as the test source.
  744. A high \cri{} can be achieved by a broadband spectrum emitters e.g. tungsten halogen lamps (\cri{}~≈~100) which has a spectrum comparable to daylight, but its efficacy is only around 25\thinspace{}lm/\textsc{w}. On contrary, the low pressure sodium lamp (\cri{}~≈~25) can reach up to 200\thinspace{}lm/\textsc{w}, but they emit only two spectral lines at 589.0\thinspace{}nm and 589.6\thinspace{}nm. Therefore every object illuminated with this lamp will be either yellow either grey or black.
  746. \marginpar{cri calculation}Calculation of \cri{} starts with defining \cct{} of a test source. For \cct{}~<~5000\thinspace{}\textsc{k} a black-body radiator with the same temperature is used. Illuminant \textsc{d} is used for sources with \cct{}~≥~5000\thinspace{}\textsc{k}. The colour differences  $\Delta{}E_{i}$ on \cie{} \abbr{1964} $u^{*}v^{*}w^{*}$ colour space of 14 selected Munsell samples (see table~\ref{tab:test_colour_samples}) is measured or calculated when illuminated with test source and reference source
  747. \begin{equation} \label{eq:colour_difference_Ei}
  748.         ΔE_{i} = \sqrt{(u^{*}_{ref}-u^{*}_{test})^{2}+(v^{*}_{ref}-v^{*}_{test})^{2}+(w^{*}_{ref}-w^{*}_{test})^{2}}.
  749. \end{equation}
  750. For each colour sample (tab.~\ref{tab:test_colour_samples}) a particular \cri{} ($R_{i}$) can be calculated using
  751. \begin{equation} \label{eq:cri}
  752.         R_i = 100 - 4.6ΔE_{i} \qquad{} (i=1, \ldots ,14).
  753. \end{equation}
  754. A general \cri{} $R_{a}$ is the arithmetic mean of the particular $R_{i}$.
  756. \marginpar{comments on cri}Although \cri{} is widely used for calculating colour quality of light sources, it is known to have some drawbacks~\cite{Ohno2004,Ohno2006}. For example only medium saturated colours are taken into account when calculating the colour rendering index, particularly $R_{9}$ (saturated red) can be very low while overall \cri{} can still be quite high. Sometimes a high \cri{} lamp e.g. 2000\thinspace{}\textsc{k} black-body ($R_{a}=100$), can have very poor colour rendering due to its very reddish light. This is also true for lamps with very high values of \cct{}. White light with slight colour tint to green or purple (laying above or below planckian locus, respectively) can have a very high \cri{} index but are unacceptable for general lighting purposes. \cri{} is now being depreciated in favour of measurements based on colour appearance models e.g. \abbr{ciecam02}~\cite{Schanda2005}. Commission Internationale de l´Eclairage recommended the development of new colour rendering index that would supplement the current \cri{}~\cite{CIE177}. Freyssinier and Rea proposed a two metric system that would combine the present \cri{} value with gamut area index forming a two-metric solution~\cite{Freyssinier2010}.
  758. \input{tables/tcs.tex}
  760. %% Colour distance
  761. \section[Colour distance]{colour distance}
  762. MacAdam's research on the ability of the human eye to notice the difference in chromaticities of two light sources started the search for uniform colour space where the perceived colour difference can be measured geometrically on the colour plane.
  764. \cie{} defines colour difference as an Euclidean distance between two colour points in a \cie{}~\abbr{1976} $L^{*}a^{*}b^{*}$ colour space
  765. \begin{equation} \label{eq:colour_difference_Eab}
  766.         ΔE_{ab}^{*} = \sqrt{(L^{*}_{ref}-L^{*}_{test})^{2}+(a^{*}_{ref}-a^{*}_{test})^{2}+(b^{*}_{ref}-b^{*}_{test})^{2}}.
  767. \end{equation}
  768. A $ΔE_{ab}^{*}$ value of 2.15 corresponds to \marginpar{just noticeable difference}just noticeable difference between two colour points \cite{Stokes1992}. Colour distance between the colour command and generated colour is a good measure of the colour control loop accuracy.
  770. Some previous work use \cie{}~1976 colour space to define the colour distance.
  771. \begin{equation} \label{eq:colour_difference_duv}
  772.         Δu'v' = \sqrt{(u'_{ref}-u'_{test})^{2}+(v'_{ref}-v'_{test})^{2}}.
  773. \end{equation}
  774. Equation does not take into account the difference in light intensity therefore it is only a measure of chromaticity difference. Colour distance equal to 0.0035 is a just noticeable difference between two colour points~\cite{Wyszecki2000}.
  776. \begin{figure}[t]% CIE_xy_MacAdam elipses
  777.                 \centering
  778. \includegraphics{graphics/ColourTheory/CIExyMacAdam.pdf}
  779. \caption[MacAdam ellipses]{MacAdam ellipses on \cie{} \abbr{1931} $xy$ chromaticity diagram. Ellipses magnified~10×. A light source with colour coordinates laying inside the ellipse is indistinguishable frmm a light source with the coordinates in the centre of the ellipse. As seen from the various sizes of the ellipses the \abbr{cie 1931} colour space is not perceptually uniform.\label{fig:CIE_xy_MacAdam}}
  780. \end{figure}
  782. % chapter colour_theory (end)
  784. %%%%%%%%%%%%%%%
  785. %% Chapter 3 %%
  786. %%%%%%%%%%%%%%%
  787. \chapter{Light-emitting diodes}% (fold)
  788. \label{cha:light_emitting_diodes}
  789. \begin{figure}[!ht]
  790.                 \centering
  791. \includegraphics[scale=1.3]{graphics/LED/HJRound.pdf}
  792. \label{fig:hjround}
  793. \end{figure}
  794. \noindent{}\led{} is a solid-state device that emits light by means of electroluminescence. First observations of light emission from SiC (carborundum) crystals were made in \abbr{1907} by Henry Joseph Round. The first observed \led{} was a Schottky diode. The light was produced on metal-semiconductor junction. Under normal, forward bias conditions, the current flowing through the Schottky diode consists only of majority carriers. The minority carriers can be injected under strong forward bias conditions (or through avalanche effect) thus making the light emission possible.
  796. %% Bandgap
  797. \section[Bandgap]{bandgap}
  798. In semiconductor materials, valence and conduction energy bands do not overlap. The energy difference between the highest point in the valence band and the lowest point in conduction band is called a bandgap. The energy levels in the bandgap are forbidden for the electrons. Transition between two bands is possible only when electron has enough surplus energy or looses enough energy—the minimal amount of energy for transition equals the bandgap energy. \marginpar{direct bandgap}In direct bandgap semiconductors, the electrons from the conduction band minima recombine with holes at the valence band maxima with preservation of momentum. Contrary, in \marginpar{indirect bandgap}indirect bandgap semiconductors, the respective extrema of conduction and valence bands do not have the same value of quantum-mechanical wave vector. Electrons have very little momentum (in the range of $2π/a$, where a is lattice constant) therefore every transition is almost vertical. Emission of a \emph{phonon}, a quantised mode of vibration of a crystal lattice with $2π/λ$ wave vector, is essential to preserve the momentum~\cite{SSPhysics}). This transition is less likely to happen than direct transition. Moreover, the temperature in the active area increases, therefore the indirect bandgap semiconductors, like silicon, are very inefficient at generating light.  
  799. \begin{figure}[h]% Direct_Indirect_Semiconductors
  800.         \centering
  801. \includegraphics{graphics/LED/DirectIndirectSemiconductors.pdf}
  802. \caption[]{Photon generation in direct (left) and indirect (right) bandgap semiconductors. In indirect bandgap semiconductors, emission of phonon must assist the radiative recombination in order to preserve the momentum \cite{PhotonicDevices}.\label{fig:Direct_Indirect_Semiconductors}}
  803. \end{figure}
  805. %% Recombination
  806. \section[Recombination]{recombination}
  807. A pair of excess carriers—free electron and hole—can recombine producing a photon. This process is called \emph{radiative recombination}. However, the recombination can also occur in a non-radiative way where no light is generated. The result of radiative recombination is the transition of an electron from higher energy state to lower energy state with emission of radiative energy. In a non-radiative recombination process, the excess energy from electron transition is converted into thermal energy in form of crystal lattice vibrations. Maximising the radiative and minimising the non-radiative processes are the goal for maximising the efficiency of  \led{}s.
  808. % Photonic devices book, Jia-Ming Liu p. 817
  810. % Bimolecular radiative recombination
  811. \noindent{}\marginpar{bimolecular radiative recombination}In light emitting diodes almost all of the radiative recombination rate is contributed by the bimolecular recombination rate~\cite{PhotonicDevices}. The net recombination rate is given by
  812. \begin{equation}
  813.  R_{net} = Bnp - Bn_{0}p_{0}
  814. \end{equation}
  815. where $B$ is a bimolecular recombination coefficient, $n$ and $p$ are electron and hole concentrations and $n_{0}$ and $p_{0}$ are electron and hole concentration under equilibrium conditions.
  817. Under high excess carrier density, when $N \gg n_{0}$, $p_{0}$  ($N = n-n_{0}$, $p-p_{0}$) the radiative recombination rate can be expressed as
  818. \begin{equation}
  819.  R = BN^2.
  820. \end{equation}
  822. % Non-radiative recombinations
  823. \begin{figure}[h]% RecombinationTypes
  824. \centering
  825. \includegraphics{graphics/LED/RecombinationTypes.pdf}
  826. \caption[Types of recombinations]{Band diagram showing possible recombinations mechanisms: deep level Shockley-Read-Hall (a), Auger~(b,c), Bimolecular~(d) band to band radiative \cite{PhotonicDevices}.\label{fig:RecombinationTypes}}
  827. \end{figure}
  829. % Deep level
  830. \noindent{}\marginpar{deep level recombination}With the presence of e.g. a defect of the crystal lattice a non-radiative recombination can occur (fig.~\ref{fig:RecombinationTypes}a). These defects act as a very efficient non-radiative recombination centres, especially if the energy level is located close to the middle of the bandgap~\cite{SemiconductorPhysics}. An electron or hole is trapped in a forbidden region of an energy band by a defect in crystal lattice. These defects can be either unintentionally introduced or engineered to produce a, so called, deep-level \led{}.
  832. The statistics of the process were described by Shockley and Read and, independently, by Hall. Net recombination rate is described by
  833. \begin{equation} \label{eq:SRHRecombinationRate}
  834.         R_{net} = \frac{np - n_{i}^{2}}{τ_{p}(n+n_{1}) + τ_{n}(p+p_{1})}.
  835. \end{equation}
  836. In doped semiconductors, equation~\ref{eq:SRHRecombinationRate} simplifies to
  837. \begin{eqnarray}
  838.         R_{net} \approx \frac{n-n_{0}}{τ_{n}}\:for\:p\gg{}n \\
  839.         R_{net} \approx \frac{p-p_{0}}{τ_{p}}\:for\:n\gg{}p
  840. \end{eqnarray}
  841. there $τ_{n}$ and $τ_{p}$ are recombination lifetimes of electrons and holes, respectively. It can be seen that recombination rate becomes proportional to the excess carrier density
  842. \begin{equation}\label{eq:RadiativeRecombinationRate}
  843.  R=AN.
  844. \end{equation}
  846. % Auger
  847. \marginpar{auger recombination}\noindent{}During Auger transition (fig.~\ref{fig:RecombinationTypes}b and c), the energy created in a recombination process between an electron and a hole is transferred to  third carrier. No light is generated and the surplus energy pushes either the second electron higher in conduction band or the second hole downward in valence band. The energy passed on to the third carrier is then dissipated by emission of phonons. The process requires three carriers to be in the same place during recombination therefore the recombination rate is proportional to the third power of carrier density
  848. \begin{equation} \label{eq:AugerRecombinationRate}
  849.         R_A = CN^3,
  850. \end{equation}
  851. where $C$ is the \emph{Auger coefficient}.
  853. Auger recombination mechanism has been proposed as one of the possible reasons for efficiency droop under high current injection in InGaN structures~\cite{Shen2007}. Although the efficiency droop is similar to Auger recombination, some argue~\cite{Hader2008,Yen2009} that the numerical values of Auger coefficient are not high enough for this process to dominate i the total power losses and other loss mechanism like electron current leakage and non uniform hole distribution may be responsible.  
  855. Taking into account equation \ref{eq:SRHRecombinationRate}, \ref{eq:RadiativeRecombinationRate} and \ref{eq:AugerRecombinationRate} we can write the stationary balance equation for \led{}s as
  856. \begin{equation} \label{eq:StationaryBalance}
  857.         j/e = AN + BN^2 + CN^3 + DN^m % A- non-rad, B - rad, C - Auger, D^m - other
  858. \end{equation}
  859. where $j$ is the injected current density, $e$ is the elementary charge, $N$ is \abbr{2d} excess carrier density in the active area and coefficients $A$, $B$, $C$ and $D$ are describing non-radiative Shockley-Read-Hall recombination, bimolecular radiative process, Auger recombination and other non-radiative processes respectively. It is beneficial to supply \led{}s with a current level, at which the bimolecular radiative process dominates.
  861. %% Materials
  862. \section[Materials]{materials}
  863. Nowadays, elemental semiconductors for solid-state lighting (\led{}s and laser diodes) are the \textsc{iii}-\textsc{v} compounds (Al, Ga and In can be found in group \textsc{iii} of periodic table, N and P in group \textsc{v}). In order to create a light generating structure a heteroepitaxy growth technology is used.
  865. \marginpar{AlGaInP}The compound crystal is grown on a substrate. The $(\mathrm{Al}_{x}\mathrm{Ga}_{1-x})_{1-y}\mathrm{In}_{y}\mathrm{P}$ is lattice matched to $\mathrm{GaAs}$ (lattice constant 5.65\thinspace{}Å) for $y=0.48$. The semiconductor has a direct bandgap for emission wavelength higher than 555~nm therefore it is used to create red, orange, amber and yellow \led{}s~\cite{AlGalnPBook}.
  867. \marginpar{InGaN}$\mathrm{In}_{x}\mathrm{Ga}_{1-x}\mathrm{N}$ nitrides lack a native bulk substrate. InGaN high power structures are typically grown on lattice mismatched sapphire (lattice constant 2.75\thinspace{}Å, 15\thinspace{}\% mismatch) or SiC substrates (lattice constant 3.08\thinspace{}Å, 3.5\thinspace{}\% mismatch)~\cite{Roussel2006}.
  869. Lattice mismatch between compound and grown crystal yields in very high dislocation densities which act as non-radiative, deep-level recombination centres. Theses centres contribute to  decreased optical output~\cite{ledSchubert}. The \textsc{iii}-\textsc{v} phosphide materials used for red, orange, amber and yellow diodes are strongly affected by the crystal defects. The \textsc{iii}-\textsc{v} nitride compounds are less affected by the crystal defects due to carrier localisation that prevents from non-radiative recombination. Therefore it is possible to create high power blue \led{} even with high number of dislocations present in a diode.
  871. %% Structures
  872. \section[Structures]{structures}
  873. Three levels of \led{} structure will be discussed. Starting from the active region that generates light in chapter~\ref{ssec:activeRegion}, through the \led{} chip that is responsible for supplying the active region with free carriers and extracting the generated light in chapter~\ref{ssec:ledChip} to the highest level structure of \led{} packaging in chapter~\ref{ssec:packaging}.
  874. \subsection{Active region}\label{ssec:activeRegion}
  875. % homojunction
  876. \marginpar{homojunction}Homojunction is a basic structure for all diodes. A connection between two differently doped semiconductors with the same bandgap forms a p-n junction. Under forward bias, free carriers from both semiconductors are forced into the regions of opposite conductivity type.
  878. \begin{figure}[h]% Homojunction
  879.         \centering
  880. \includegraphics{graphics/LED/Homojunction.pdf}
  881. \caption[Homojunction]{Band diagram showing p-n junction under no bias (left) and forward bias (right) conditions. Under forward bias, the recombination occurs in the area around the junction defined by diffusion lengths. \label{fig:Homojunction}}
  882. \end{figure}
  884. \noindent{}Diffusion lengths ($L_n$ and $L_p$) of minority carriers depends on carrier mobility.
  885. \begin{equation} \label{eq:DiffusionLength}
  886.         L_{n}=\sqrt{D_{n}τ_{n}},\qquad{} L_{p}=\sqrt{D_{p}τ_{p}},
  887. \end{equation}
  888. where $τ_{n}$ and $τ_{p}$ are electron and hole minority carrier lifetimes. $D_n$ and $D_p$ are the diffusion constants
  889. \begin{equation} \label{eq:EinsteinRelation}
  890.         D_{n}={\frac{kT}{e}}\mu{}_{n},\qquad{} D_{p}={\frac{kT}{e}}\mu{}_{p}
  891. \end{equation}
  892. where $\mu{}_{n}$ and $\mu{}_{n}$ are electron and hole mobility dependent on the semiconductor.
  894. Free carriers recombine radiatively around the junction area within the diffusion distance. As recombination rate is given by the bimolecular recombination equation $R=Bnp$, where $n$ and $p$ are the electron and hole carrier concentrations and $B$ is bimolecular recombination coefficient. The light generation is limited in this structure by the low excess carrier density. The width of the $p$ region is critical for obtaining higher efficiencies. When the region is too shallow, the electrons can escape and recombine non-radiatively through crystal defects. Too thick region will easily reabsorb emitted photons thus reducing the overall efficiency. A heterojunction can be used to overcome both this problems at the same time.
  896. % double heterojunction
  897. \marginpar{double heterostructure}
  898. A third semiconductor is added between two, differently doped, semiconductors to form a confinement area for free carriers. This third semiconductor has lower bandgap than surrounding materials. As a result, the excess carriers get trapped between two barriers. Width of this barrier can be controlled during manufacturing process. Barrier width is much lower than diffusion lengths in the particular semiconductor therefore free carriers are confined in much smaller area as compared to homojunction structure. Increased carrier densities yield in highly increased light generation.
  899. \begin{figure}[h]% Double Heterostructure
  900.         \centering
  901. \includegraphics{graphics/LED/DoubleHeterostructure.pdf}
  902. \caption[Double heterostructure]{Band diagram showing double heterostructure under forward bias conditions. Free carriers are confined between two barriers. \label{fig:DoubleHeterostructure}}
  903. \end{figure}
  905. \noindent{}Free carriers are distributed within active region with Fermi-Dirac distribution. Therefore carriers with energy levels higher than barrier height can escape from the confinement area. Concentration of these carriers can be calculated as
  906. \begin{equation} \label{eq:LeakageElectrons}
  907.         n_{B} = \int_{E_{B}}^{\infty{}}\!\!\!\!ρ_{dos}\,f_{FD}\ud{}E
  908. \end{equation}
  909. where $ρ_{dos}$ is the density of states, $f_{FD}$ is the Fermi-Dirac distribution and $E_{B}$ is the height of the heterojunction barrier. This leakage contributes to decreased efficiency, especially at high temperatures. An electron blocking layer can be introduced into the structure to help confine free carriers.
  910. \begin{figure}[t]% Electron Blocking Layer
  911.         \centering
  912.         \includegraphics{graphics/LED/ElectronBlockingLayer.pdf}
  913. \caption[Double heterostructure]{Band diagram showing double heterostructure without (left) and with (right) electron blocking layer. Additional layer prevents the electrons escaping from the active region, decreasing the leakage current thus increasing the diode efficiency \cite{ledSchubert}. \label{fig:ElectronBlockingLayer}}
  914. \end{figure}
  916. \noindent{}Bigger width of an active area width than in quantum well structures makes the carrier density lower and therefore decreases the impact of Auger recombination on diode efficiency at higher currents.  % droop
  918. % MQW structures
  919. \marginpar{quantum well structures}When additional, thin (in the range of tens of angstroms) layer of semiconductor with smaller bandgap than the cladding layers is added a quantum well is formed. A quantum well structure confines the excess carriers to a very narrow area in order to increase the carrier concentration and therefore increase the recombination rate. Narrow width of the active area also reduces the self absorption effects in some \led{} structures~\cite{Hassan2005}.
  921. \marginpar{allowed energy levels}When the thickness of the well is close to the de~Broglie wavelength, the quantum effects appear. Solving the Schrödinger equation for a finite potential well produces the allowed energy levels values within the quantum well. The energy states are finite and are dependent on the well width. Therefore by varying the width of the well one can change the colour of an \led{}.
  923. \marginpar{saturation effect}One of the implication of using quantum wells is that they have limited capacity and, at high current injection, a carrier overflow may occur. When carrier concentration in the active area reaches certain level the Fermi energy rises above the top of quantum well. The carrier overflow the structure and further increase of carrier density in the active area is impossible. The light output becomes saturated. In order \marginpar{mqw structures}to overcome this problem multiple quantum well structures are used~\cite{ledSchubert} as shown in figure~\ref{fig:ElectronBlockingLayer}.
  925. \begin{figure}[h]% MQW
  926.         \centering
  927. \includegraphics{graphics/LED/MQW.pdf}
  928. \caption{Band diagram showing multiple quantum wells structure under forward bias conditions (left) and  details of single quantum well (right). A thick semiconductor layer with energy gap $E_{g2}$ is surrounded by a semiconductor with a higher bandgap $E_{g2}$ forming a quantum well. Allowable carrier energy values, subbands, shown inside the potential well \cite{ledSchubert}. \label{fig:MQW}}
  929. \end{figure}
  931. Carrier density on each of available energy states depends on the injection current density. At low values mostly the lowest energy states are occupied by electrons in conduction band well and the highest states in valence band.
  933. \subsection{\led{} chip}\label{ssec:ledChip}
  934. Due to low efficiency of homojunction structure, high power \led{}s use either double heterostructure or \textsc{mqw} structures for active regions. These active structures are embedded in a \led{} chip, whose task is to supply active region with free carriers, extract the generated light as efficiently as possible.
  936. \begin{figure}[h]% LED chip structures
  937.         \centering
  938. \includegraphics{graphics/LED/DiodesStructures.pdf}
  939. \caption{Common structures of high power \led{}s: truncated-inverted pyramid (left) for phosphide diodes~\cite{Krames1999}, flip-chip (middle) gallium nitride structure~\cite{Wierer2001} and thin-film flip-chip (right) GaN structure~\cite{Shchekin2006}.\label{fig:LEDChipStructures}}
  940. \end{figure}
  942. \noindent{}Supplying the current to the active region is done utilising metallic contacts  and conductive layers of semiconductors acting as current spreading layers.  
  944. \marginpar{light extraction}Light extraction is a very difficult task as the \led{} chips are built using materials with high refractive index causing total internal reflection phenomenon trapping the light inside the chip. Generated photons can be reabsorbed by the active area or by metallic contacts therefore they should be extracted from the chip as fast as possible. Special geometric shapes (fig.~\ref{fig:LEDChipStructures}) of the \led{} chip as well as micro patterning of edges allow the light to escape the structure~\cite{Krames2007}. Bottom contacts are made reflective to maximise the light output. For example, the thin-film flip-chip structure is made by growing the GaN structure on a sapphire substrate, flipping it and chemically removing the transparent substrate for increased optical transmittance. The revealed n-type GaN lager is then roughened through photochemical etching to increase light extraction~\cite{Shchekin2006,ledSchubert}.
  946. \subsection{Packaging}\label{ssec:packaging}
  947. \led{} packages very often include a lens that shape the light beam. The lens and plastic body of the diode is not thermally conductive, therefore the packaging must provide good thermal contact between \led{} chip and mounting surface. Some \led{} chips must be electrically isolated from the heatsink. This provides a big challenge as good electrical isolators are typically also good thermal isolators.
  949. %% Electrical characteristics
  950. \section[Electrical characteristics]{electrical characteristics}
  951. Light in the \led{} structures is created by electroluminescence—this means that the energy created during electron-hole recombination is directly converted to photons. Therefore forward voltage of an \led{} must be at least equal to the bandgap energy in the active region divided by the elementary charge.
  952. \begin{equation} \label{eq:ForwardVoltage}
  953.         V_{f} \geq{} \frac{}{e} \approx{} \frac{E_{g}}{e}
  954. \end{equation}
  955. The value of forward voltage can be increased by few phenomena. The series resistance will increase the diode drive voltage linearly with the current. Total diode series resistance is a sum of contact resistance, resistance caused by abrupt band structure and bulk resistance of \led{} materials. Additionally, carriers injected into quantum well loose some energy because discrete energy levels inside the well will be always smaller than the cladding layer energy. Therefore difference between free carrier energy in confinement layer and inside a well has to be dissipated by a phonon emission. The forward voltage of an \led{} can be expressed as
  956. \begin{equation} \label{eq:DetailedForwardVoltage}
  957.         V_{f} = \frac{E_{g}}{e} + IR_{s}+ \frac{\Delta{}E_{C}-E_{0}}{e} + \frac{\Delta{}E_{V}-E_{0}}{e}
  958. \end{equation}
  959. where $\Delta{}E_{C}-E_{0}$ and $\Delta{}E_{V}-E_{0}$ represent differences in energy levels between quantum well and conduction and valence confinement layers, respectively~\cite{Hassan2005}. \marginpar{shockley equation}Electrical characteristics of a p-n junction are described by Shockley equation \cite{PhotonicDevices}
  960. \begin{equation}
  961.  I = I_S \big( \mathrm{exp}(V/nV_T) - 1\big)
  962. \end{equation}
  963. where $I$ and $V$ are diodes current and voltage respectively and $n$ is the ideality factor. $I_S$ and $V_T$ are the reverse bias saturation current and thermal voltage described by
  964. \begin{equation}
  965.  I_S = eA \Bigg( \sqrt{ \frac{D_p}{τ_p} } \cdot \frac{n_i^2}{N_D} + \sqrt{ \frac{D_n}{τ_n} } \cdot \frac{n_i^2}{N_A}\Bigg), \qquad V_T = \frac{kT}{e}.
  966. \end{equation}
  967. Including parasitic series resistance and assuming forward-bias conditions, when $V \gg kT/e$ the \textsc{iv} characteristics can be rewritten as
  968. \begin{equation}\label{eq:iv_diode_model}
  969.  I = I_S\Big( \mathrm{exp}\big(e(V-IR_s)/nkT\,\big)\Big).
  970. \end{equation}
  972. \begin{figure}[t]% IV characteristics
  973.         \centering
  974. \includegraphics{graphics/LED/IV/IdealityFactorSeriesResistance/iv_rgbw.pdf}
  975. \caption{Measured current-voltage characteristics of red, green, blue and white Luxeon \textsc{k2} diodes at 20°\textsc{c}.\label{fig:IV_n_Rs}}
  976. \end{figure}
  978. \noindent{}Schockley's equations does not include the presence of quantum well(s) in the diode structure and tunneling injection. Despite this, it is a common practice to fit measurement data to Shockley equation even though it does not reflect all physical phenomena in the diode~\cite{Lee2006b}.
  980. In order to compare the theoretical electrical model for a single p-n junction with structures used in high intensity light-emitting diodes available on the market the characteristics of basic red, green, blue and white colour diodes were measured using equipment described in chapter \ref{measurementHardware}. Diodes were placed on a thermally controlled cold plate and charged capacitor was connected to their terminals. The capacitor was discharged through the diode providing short, low energy pulse (thus minimising the heating of the junction).
  982. Figure \ref{fig:IV_n_Rs} shows the current-voltage characteristics of four high power Luxeon \textsc{k2} diodes. The data was fitted to theoretical model described by equation~\ref{eq:iv_diode_model}. Identified values of ideality factor and series resistance are gathered in table \ref{tab:n_Rs}. Expected values of ideality factor are in the 1–2 range. However, measured values of the ideality factor are much higher, in the range of 3–7. This is attributed to tunneling injection.
  984. \begin{table}[!ht]\footnotesize
  985. \caption{Measured values of ideality factor and series resistance of Luxeon \textsc{k2} diodes.\label{tab:n_Rs}}
  986. \centering
  987. \begin{tabular}{lcc}
  988. \textsc{diode} & $n$ & $R_s$ \\
  989. \hline
  990. \abbr{lxk2-pd12-r00} (red)   & 3.166 & 2.334\thinspace{}Ω \\
  991. \abbr{lxk2-pm14-u00} (green) & 7.349 & 0.898\thinspace{}Ω \\
  992. \abbr{lxk2-pb12-k00} (blue)  & 6.014 & 0.775\thinspace{}Ω \\
  993. \abbr{lxk2-pw12-u00} (white) & 4.563 & 0.718\thinspace{}Ω \\
  994. \hline
  995. \end{tabular}
  996. \end{table}
  998. %% Optical characteristics
  999. \section[Optical characteristics]{optical characteristics}
  1000. Typical emission spectra of AlGaInP and InGaN diodes are shown in figures~\ref{fig:MeasuredRGBSpectra_nm} and \ref{fig:MeasuredRGBSpectra_eV}. Green gallium nitride diodes have the widest emission spectra due to high fluctuations of chemical composition in indium rich semiconductor~\cite{HandbookOfOpticalMetrology}.
  1002. \begin{figure}[!h]% Measured emission spectra of R, G, B diodes
  1003.         \centering
  1004.                 \includegraphics{graphics/LED/Spectrum/RGBspectra.pdf}
  1005.         \caption{Measured emission spectra of red AlGaInP and green and blue InGaN diodes with normalised intensities. \fwhm{} is equal to 20, 33 and 21.3\thinspace{}nm for red green and blue diode respectively.}
  1006.         \label{fig:MeasuredRGBSpectra_nm}
  1007. \end{figure}
  1009. \begin{figure}[!h]% Measured emission spectra of R, G, B diodes
  1010.         \centering
  1011.         \includegraphics{graphics/LED/Spectrum/RGBspectra_eV.pdf}
  1012.         \caption{\led{} emission spectra (fig.~\ref{fig:MeasuredRGBSpectra_nm}) shown in energy scale. \fwhm{} is equal to 2.41, 6.00 and 4.93\thinspace{}\textit{kT} for red green and blue diode respectively. Theoretical \fwhm{} value for \led{}s is equal to 1.8\thinspace{}\textit{kT}.}
  1013.         \label{fig:MeasuredRGBSpectra_eV}
  1014. \end{figure}
  1015. \clearpage % *** manual new page 43/44
  1016. % emission spectrum
  1017. \noindent{}\marginpar{theoretical emission spectrum}Theoretical emission spectrum of a bulk semiconductor \led{} is a product of the density of states and the distribution of carriers in the allowed bands characterised by a Boltzmann distribution \cite{ledSchubert}.
  1018. \begin{equation} \label{eq:TheoreticalemissionSpectrum}
  1019.         I(E) = ρ(E) \cdot f_{B}(E) \propto \sqrt{E-E_{g}} \cdot e^{-E/kT}
  1020. \end{equation}
  1021. \noindent{}This equation is to be used with care for InGaN diodes, since Boltzmann distribution is not applicable to localised carriers.
  1023. \begin{figure}[!ht]% Theoretical emission spectrum
  1024.         \centering
  1025.                 \includegraphics{graphics/LED/Spectrum/TheoreticalSpectrum.pdf}
  1026.         \caption{Theoretical emission spectrum as a product of Boltzmann distribution and the density of states \cite{ledSchubert} (left) and measured spectrum of a green \led{} (right), both plotted in linear (top) and log-lin scales (bottom).}
  1027.         \label{fig:TheoreticalSpectrum}
  1028. \end{figure}
  1030. \noindent{}\marginpar{full width at half maximum}The theoretical maximum of emission occurs at $E=E_{g}+kT/2$ and the value of the full width at half maximum is 1.8\thinspace{}$kT$ (figure~\ref{fig:TheoreticalSpectrum}). Practically, the spectral width of transitions in \led{} is wider, in the range of 2–4\thinspace{}$kT$ due to alloy broadening~\cite{ledSchubert}. This is clearly visible when comparing theoretical spectrum with a measured one (fig.~\ref{fig:TheoreticalSpectrum}). Practical diodes exhibit photon emission with energy much lower than predicted by the theory. Particularly, in InGaN diodes, where the bandgap energy is dependent on the indium composition and any fluctuations in its doping can cause broadening of the spectrum.
  1032. Plotting the theoretical emission spectrum on a semi logarithmic plot shows that the high energy side of the spectrum is linear and dependent on the junction temperature.
  1034. The shape of emission spectra is sometimes approximated by a gaussian distribution for modelling purposes \cite{Chhajed2005b}
  1035. \begin{equation} \label{eq:GaussianDistribution}
  1036.         I(λ) \propto{} \frac{1}{σ\sqrt{}} \mathrm{exp}\bigg(\!\!-\frac{(λ-λ_{peak})^2}{2σ^2}\bigg)
  1037. \end{equation}
  1038. where standard deviation relates to \fwhm{} as~\cite[p.~335]{ledSchubert}
  1039. \begin{equation} \label{eq:GaussianDistributionSigma}
  1040.         σ = \frac{Δλ}{2\sqrt{2ln2}} \approx \frac{Δλ}{2.335} .
  1041. \end{equation}
  1042. Any mismatch between the model of the spectrum and the actual \led{} spectrum may introduce significant colour point difference. The actual spectrum of an \led{} is not symmetrical therefore gaussian model can introduce significant errors. To overcome this problem, Man and Ashdown~\cite{Man2006} proposed a double gaussian model: a sum of two gaussian curves that provide much better empirical model than a single curve.
  1044. %% Thermal characteristics
  1045. \section[Thermal characteristics]{Thermal characteristics}
  1046. \led{}s show a very high temperature dependency. Because the junction temperature can vary significantly during normal operation (thermal cycling during low frequency \pwm{}, dimming of the diode, heatsink and ambient temperature change) it is necessary to investigate the phenomena thoroughly. Not only the efficiency drops with increased temperature, but the colour point can move far outside the MacAdam ellipsis from the original position. Colour point shift occurs due to the change in spectrum shape. In order to quantify the spectral shift, four basic colour diodes were measured at different heatsink temperatures. The resulting spectra are gathered in figure~\ref{fig:Spectra2Temperature}. The results show AlGaInP red diode experiences the biggest spectrum change while the spectrum shape of InGaN diodes is not affected by the change of temperature.
  1048. \begin{figure}[!ht]% Spectra vs. temperature
  1049.         \centering
  1050.                 \includegraphics{graphics/LED/Thermal/SpectraChangeRGBW.pdf}
  1051.         \caption{Measured change in spectral shape for Luxeon Rebel red (left), green, blue and white (right) diodes for varying heatsink temperature.}
  1052.         \label{fig:Spectra2Temperature}
  1053. \end{figure}
  1055. \noindent{}The decrease in luminous output is connected to several mechanisms: increase of non-radiative recombinations (deep level and surface) and carrier losses due to leakage over heterostructure barriers. InGaN diodes, however, experience high efficiency droop at high current value. The physical origin of the droop is not yet known~\cite{Schubert2007}. Electron leakage, lack of hole injection, carrier delocalisation, Auger recombination, defects, and junction heating were suggested as explanations for this phenomenon~\cite{Yen2009,Schubert2009}.
  1057. \begin{figure}[p]% Flux, Vf, Prad, FWHM, efficiency, efficacy vs. temperature
  1058.         \centering
  1059.         \includegraphics{graphics/LED/Thermal/parameter_variation_with_temperature_K2.pdf}
  1060.         \caption{Measured Luxeon \textsc{k2} diodes' parameters with varying heatsink temperature.}
  1061.         \label{fig:ParametersVsTemperatureK2}
  1062. \end{figure}
  1064. \begin{figure}[p]% Flux, Vf, Prad, FWHM, efficiency, efficacy vs. temperature
  1065.         \centering
  1066.         \includegraphics{graphics/LED/Thermal/parameter_variation_with_temperature_rebel.pdf}
  1067.         \caption{Measured Luxeon Rebel diodes' parameters with varying heatsink temperature.}
  1068.         \label{fig:ParametersVsTemperatureRebel}
  1069. \end{figure}
  1071. The rise of heatsink temperature corresponds to the rise of junction temperature. The relationship between junction and heatsink temperatures is not linear because the losses in the diode depend on the junction temperature. Because of the complexity of junction temperature measurement, the properties of various \led{}s have been analysed with respect to the heatsink temperature. Diodes' properties are plotted in figures~\ref{fig:ParametersVsTemperatureK2} and~\ref{fig:ParametersVsTemperatureRebel} for Luxeon \textsc{k2} and Luxeon Rebel diodes, respectively. The rate of change of these parameters can be approximated with a linear temperature dependencies. Temperature coefficients are gathered in table~\ref{tab:diodeParametersVsTemp}. In order to collect the necessary data the diodes were placed on thermally controlled plate and mounted to the integrating sphere. The current was set to 700\thinspace{}mA using a \dc{} regulated current source. The heatsink temperature was regulated in 20–60°\textsc{c} range. After each change in temperature the system was allowed to reach thermal steady state after which diode's parameters were recorded.
  1073. The emissivity of the diode depends strongly on the junction temperature. Figures \ref{fig:ParametersVsTemperatureK2}\thinspace{}b and~\ref{fig:ParametersVsTemperatureRebel}\thinspace{}b show the drop of radiometric power with the increase of heatsink temperature.
  1075. The spectrum of all diode shifts towards longer wavelengths (fig.~\ref{fig:ParametersVsTemperatureK2}\thinspace{}e and~\ref{fig:ParametersVsTemperatureRebel}\thinspace{}e). Both these phenomenon will influence the perceived flux (fig.~\ref{fig:ParametersVsTemperatureK2}\thinspace{}a and~\ref{fig:ParametersVsTemperatureRebel}\thinspace{}a). Shift towards longer wavelengths will increase the luminous flux in blue and green diodes but will decrease it in case of red diodes. This is because the peak of eye sensitivity is at 555\thinspace{}nm (fig.~\ref{eye_sensitivity}) and moving the spectrum towards this wavelength will increase the flux and moving away from it will decrease the flux.
  1077. Luminous flux figures show that the red AlGaInP diodes are mostly affected by the temperature change because of the aforementioned phenomena. Ca. 40\thinspace{}\% of the flux is lost with the heatsink temperature increase by 40°\textsc{c} for both diode families. Nitride diodes show little temperature dependence: less than 8\thinspace{}\% of the flux is lost. Phosphor converted white diodes experience bigger loss of intensity than pure GaN diodes. With the 40°\textsc{c} temperature rise around 13\thinspace{}\% was lost for both diode families.
  1079. \begin{table}[!ht]\footnotesize
  1080. \caption{Temperature coefficients of forward voltage, luminous flux, radiometric power, peak wavelength and \fwhm{} for the measured Luxeon diodes.\label{tab:diodeParametersVsTemp}}
  1081. \centering
  1082. \begin{tabular}{lccccc}
  1083. & $\ud{}V_{f}/\ud{}T$           & $\ud{}F/\ud{}T$               & $\ud{}P_{rad}/\ud{}T$                 & $\ud{}peak/\ud{}T$    & $\ud{}\textsc{fwhm}/\ud{}T$           \\
  1084. \textsc{diode}  & [mV/K] & [lm/K]       & [mW/K]        & [nm/K] & [nm/K]       \\
  1085. \hline
  1086. Luxeon \textsc{k2} red          &       −2.29 & −0.665 & −2.57 & 0.154 & 0.050 \\
  1087. Luxeon \textsc{k2} green        &       −2.68 & −0.122 & −0.35 & 0.052 & 0.059 \\
  1088. Luxeon \textsc{k2} blue         &       −4.71 & −0.033 & −0.86 & 0.028 & 0.055 \\
  1089. Luxeon \textsc{k2} white        &       −3.00 &       −0.266 & −0.80 & 0.034 & 0.111 \\
  1090. \hline
  1091. Luxeon Rebel red                                        &       −1.49 & −0.516 & −2.76 & 0.166 & 0.054 \\
  1092. Luxeon Rebel green                              &       −2.90 & −0.113 & −0.32 & 0.050 & 0.055 \\
  1093. Luxeon Rebel blue                                       &       −3.44 & −0.024 & −0.79 & 0.032 & 0.050 \\
  1094. Luxeon Rebel blue                                       &       −4.12 & −0.243 & −0.68 & 0.047 & 0.018 \\
  1095. \hline
  1096. \end{tabular}
  1097. \end{table}
  1098. \clearpage % *** manual new page 48/49
  1099. % Junction temperature measurements
  1100. \noindent{}\marginpar{junction temperature measurements}Many techniques may be used to obtain junction temperature on \led{}. Most common methods utilise forward voltage $V_{f}$, high energy slope and peak wavelength temperature dependency. High energy slope method requires knowledge of emission spectra of the diode (figure~\ref{fig:TheoreticalSpectrum}).
  1102. % Forward voltage temperature dependency
  1103. \subsection{Forward voltage}
  1104. The basic equation describing the temperature dependence of the forward voltage~\cite{ledSchubert,Xi2004,Xi2005} consists of three summands.
  1105. \begin{equation} \label{eq:ForwardVoltageVsTemp}
  1106.         \frac{\ud{}V_{f}}{\ud{}T} = \frac{eV_{f}/E_{g}}{eT} + \frac{1}{e}\frac{\ud{}E_{g}}{\ud{}T} - \frac{3k}{e}
  1107. \end{equation}
  1108. First describes the influence of intrinsic carrier concentration. Second applies to the temperature dependency of the bandgap. Last summand a temperature dependency of effective densities of states $N_{C}$ and $N_{V}$ in conduction and valence bands, respectively and acceptor and donor concentrations, $N_{A}$ and $N_{D}$, respectively~\cite{MicroelectronicsWhitaker}
  1109. \begin{equation}
  1110.         eV_{f}-E_{g} \approx kT\ln{\frac{N_{D}N_{A}}{N_{C}N_{V}}}.
  1111. \end{equation}
  1113. \noindent{}The bandgap energy can be approximated in certain temperature region by the Varnshi parameters~\cite{Varshni1967,Vainshtein1999}. The linear-quadratic relation describing temperature dependence of the bandgap uses two empirical parameters \textit{α}, \textit{β} and the width of a bandgap at absolute zero temperature $E_{0}$ (table~\ref{tab:varnshi_parameters})
  1114. \begin{equation} \label{eq:BandgapVarshni}
  1115.         E_{g} = E_{0} - \frac{α^{2}}{β+T}.
  1116. \end{equation}
  1118. \input{tables/Varshni_parameters.tex}
  1120. \noindent{}To obtain the temperature dependence of the energy gap in ternary semiconductor diodes (eg.~$\mathrm{In}_{x}\mathrm{Ga}_{1-x}\mathrm{N}$) a quadratic fitting is used
  1121. \begin{equation} \label{eq:tetraryBandgapEnergy}
  1122.         E_{g}(\mathrm{A}_{1-x}\mathrm{B}_{x}) = (1-x)E_{g}(\mathrm{A}) + xE_{g}(\mathrm{B}) - x(1-x)C.
  1123. \end{equation}
  1124. Similarly, the bandgap energy of quaternary $(\mathrm{Al}_{x}\mathrm{Ga}_{1-x})_{0.52}\mathrm{In}_{0.48}\mathrm{P}$ barrier layers can be expressed as
  1125. \begin{equation} \label{eq:quaternaryBangapEnergy}
  1126.         E_{g}(x) = xE_{g}(\mathrm{Al}_{0.52}\mathrm{In}_{0.48}\mathrm{P}) + (1-x)E_{g}(\mathrm{Ga}_{0.52}\mathrm{In}_{0.48}\mathrm{P}) - x(1-x)C,
  1127. \end{equation}
  1128. where $C=0.18\:\mathrm{eV}$ is a bowing parameter.
  1130. \begin{figure}[t]% IV characteristics
  1131.         \centering
  1132. \includegraphics{graphics/LED/IV/TemperatureDependency/iv_rgbw.pdf}
  1133. \caption{Current-voltage characteristics of red, green, blue and white diodes under varying heatsink temperature from 20°\textsc{c} to 50°\textsc{c}, every ten degrees.\label{fig:IV_temperature_dependency}}
  1134. \end{figure}
  1136. \noindent{}Because of the forward voltage temperature dependency \led{}s cannot be controlled easily with a voltage source. Self heating effect causes the drop in forward voltage moving the operating point to a higher current. Increased current yields an increase of temperature. This positive feedback may lead to thermal runaway and the destruction of the diode.
  1138. % Thermal model
  1139. \subsection{Thermal model}
  1140. In order to estimate junction temperature of \led{} manufacturers provide a junction to case thermal resistance. Based on its value and power dissipated in the device as heat one may calculate temperature difference between junction and diode's case. Assuming the thermal resistance of case to heatsink interface to be negligible or known the junction temperature can be calculated with the reference to measured heatsink temperature. However, the calculated temperature will only be valid in thermal equilibrium.
  1142. \marginpar{thermal impedance}Poppe~et~al.~\cite{Poppe2009,Poppe2010a,Poppe2010b,Poppe2010c} analysed transient response of a diode during a cool down process. By analysing this response, the thermal impedance of the \led{} thermal structure can be revealed.
  1144. \begin{figure}[!ht]% Thermal impedance curves
  1145. \centering
  1146. \includegraphics{graphics/LED/Thermal/luminus_rebel_thermal_impedance.pdf}
  1147. \caption{Thermal impedance curves for Luminus Rebel InGaN diode, after Parry and Rose. \label{fig:ThermalImpedanceCurves}}
  1148. \end{figure}
  1150. \noindent{}Each element of \led{} structure represents certain thermal resistance and capacitance on the heat flow path and is connected to a time constant in thermal system response. Overall thermal system response is a sum of individual element responses~\cite{Szekely1988}.
  1151. \begin{figure}[!ht]% RC network thermal model
  1152. \centering
  1153. \includegraphics{graphics/LED/Thermal/RC_Network.pdf}
  1154. \caption{Cauer (top) and Foster (bottom) \abbr{rc} network models of heat flow path. Current sources represent heat generation. Voltage source represents reference temperature e.g. ambient or heatsink temperature. \label{fig:RC_network_thermal_model}}
  1155. \end{figure}
  1156. Both thermal models presented in figure~\ref{fig:RC_network_thermal_model} are equivalent but the Cauer structure represents physical heat flow. By reducing the number of elements in the \textsc{rc} ladder, the physical structure of the diode can be modelled using thermal capacitances of the \led{} die, die attach and heat slug and thermal resistances of the interfaces between them. High power \led{}s are mounted on heatsinks using a thermal interface material that, due to its thin layer, can be modelled using only thermal resistance~\cite{Hui2009}. The heatsink itself typically has big thermal capacitance and resistance compared to \led{} structure.
  1158. \begin{figure}[!ht]% Luminus Rebel thermal model
  1159. \centering
  1160. \includegraphics{graphics/LED/Thermal/luminus_rebel_thermal_model.pdf}
  1161. \caption{Simplified thermal model of Luminus Rebel diodes. Component values can be derived from thermal impedance curves (fig.~\ref{fig:ThermalImpedanceCurves}). \label{fig:LuminusRebelThermalModel}}
  1162. \end{figure}
  1164. \noindent{}The value of junction to case thermal resistance has been found to be nonlinearly dependent on the ambient temperature and power dissipated in the junction area~\cite{Jayasinghe2007}. The magnitude and the direction of thermal resistance change depends on the structure of diode~\cite{ShengLiang2008}. This change is attributed to changing properties of thermal interface material and varying thermal properties of semiconductors~\cite{Jayasinghe2007}.
  1166. \marginpar{effect of thermal resistance on luminous flux}The heat flow path is partially determined by the internal structure of the diode and partially by the accompanying heatsink. In order to test the influence of the heatsink properties on the static properties of the diode a simulation was set up.
  1167. \begin{figure}[!ht]%
  1168. \centering
  1169. \includegraphics{graphics/LED/Thermal/Rth_variation_model.pdf}
  1170. \caption{Simulation model for quantifying the effect of heatsink-to-ambient thermal resistance on maximum luminous flux in steady-state conditions. \label{fig:Rth(hs-a)ChangeModel}}
  1171. \end{figure}
  1172. The static thermal model in the simulation shown in figure~\ref{fig:Rth(hs-a)ChangeModel} contains only thermal resistances shown in figure~\ref{fig:LuminusRebelThermalModel}. The capacitances are omitted as they contribute only to dynamic properties of the thermal system. Junction temperature and driving current are used to determine the forward voltage of the diode. Next, the current-voltage information is converted into flux and dissipated power using \led{} model described in detail in chapter~\ref{sec:iv_model}. The flux is plotted with respect to the driving current for different values of heatsink-to-ambient thermal resistance.
  1173. \begin{figure}[!ht]%
  1174. \centering
  1175. \includegraphics{graphics/LED/Thermal/effect_of_varying_Rth.pdf}
  1176. \caption{Effect of changing heatsink-to-ambient thermal resistance on luminous flux in steady-state conditions for InGaN (left) and AlGaInP (right) diode. \label{fig:Rth(hs-a)ChangeEffect}}
  1177. \end{figure}
  1179. \noindent{}As seen in figure~\ref{fig:Rth(hs-a)ChangeEffect}, the InGaN diodes are not as sensitive to the thermal resistance as AlGaInP diodes. Therefore the heatsink thermal resistance and \textsc{tim} thermal resistance will limit the amount of luminous flux in steady-state conditions~\cite{Hui2009}. In case of AlGaInP diodes too high thermal resistance value causes the $\ud{}F/\ud{}I$ to become negative. This creates dangerous effect where increasing the driving current lowers the luminous flux. In a colour control system that regulates the diode's flux this phenomenon can lead to thermal runaway and destruction of the luminaire. Results of the simulation show the behaviour of changing the thermal resistance of the heat flow path. This resistance can change as an effect of curing or ageing of \textsc{tim} materials~\cite{Prasher2006}, airborne contaminants that cover the surface of the heatsink~\cite{Cirolia2001}, change of fan speed in forced cooling or change of efficiency of convection cooling (increases at higher heatsink temperatures). All aforementioned processes will influence the junction temperature of the diodes therefore a proper safety margin in designing the cooling system should be incorporated.
  1181. %% Efficiency
  1182. \section[Efficiency]{efficiency}
  1183. An ideal source of light would produce one photon for every electron injected. Because of non-radiative recombination phenomenon, not every electron-hole pair produce a quantum of light. Therefore we can define \emph{internal quantum efficiency} as a ratio of the number of photons emitted from the active region of a \led{} per quantum of time to the number of electrons injected into active region per quantum of time
  1184. \begin{equation} \label{eq:InternalQuantumEfficiency}
  1185.         η_{int} = \frac{n_{\textit{created photons}}}{n_{\textit{electrons}}} = \frac{P_{int}/hν}{I/e}
  1186. \end{equation}
  1187. where $P_{int}$ is the optical power emitted from the active region and $I$ is the injection current.
  1189. The photons created in the active region have to be extracted into free space. This can be very challenging as the photons are being created in \textsc{iii}-\textsc{v} crystals with high refractive indices (n~$\approx{}3.5$ for AlGaInP and n~$\approx{}2.4$ for InGaN materials). Some light extracting features like texturing of a crystal surface has to be incorporated into \led{} design. At the same time, the light radiation characteristic has to meet the demands of the market. Therefore development of a \led{} packages with high light extraction efficiency $C_{ext}$ is still an ongoing process. Currently (January 2009) $C_{ext}$ reaches up to 60\% in AlGaInP and 80\% in InGaN high power devices.
  1191. The \emph{external quantum efficiency} is designed as a ratio of the number of photons emitted into the free space per quantum of time to the number of electrons injected into active region per quantum of time
  1192. \begin{equation} \label{eq:ExternalQuantumEfficiency}
  1193.         η_{ext} = \frac{n_{\textit{extracted photons}}}{n_{\textit{electrons}}} = \frac{P/hν}{I/e} = C_{ext} \cdot{} η_{int}
  1194. \end{equation}
  1196. \noindent{}A \emph{power efficiency}, often refereed to as \emph{wall plug efficiency} is defined as
  1197. \begin{equation} \label{eq:PowerEfficiency}
  1198.         η_{power} = \frac{P_{optical}}{P_{electrical}}
  1199. \end{equation}
  1200. where $P_{electrical}$ is electric power supplied to the diode and $P_{optical}$ is the optical power of light extracted from the structure.
  1202. A \emph{luminous efficacy} is a measure of how well the electromagnetic radiation is converted into luminous flux.
  1203. \begin{equation} \label{eq:LuminousEfficacy}
  1204.         η_{lum} = \frac{F_{lum}}{P_{optical}} =  \frac{683\int{}\!V(λ)P(λ)\ud{}λ}{\int{}\!P(λ)\ud{}λ}
  1205. \end{equation}
  1206. The highest possible luminous efficacy is defined by the eye sensitivity curve [fig.~\ref{eye_sensitivity}] and is equal to 683\thinspace{}lm/\textsc{w} when the light is a monochromatic source emitting at 555\thinspace{}nm.
  1208. In lighting applications it is necessary to know how well the input electrical power is converted into visible light. This is called \emph{luminous efficacy of a source} or \emph{wall plug efficacy}.
  1209. \begin{equation} \label{eq:LuminousEfficacyWallPlug}
  1210.         η_{lum,wp} = \frac{F_{lum}}{P_{electrical}} = η_{lum} \cdot{} η_{power} =  \frac{683\int{}\!V(λ)P(λ)\ud{}λ}{IV}
  1211. \end{equation}
  1212. Highest practically achievable wall-plug efficacies, for various light sources, has been collected in table~\ref{tab:light_sources'_efficacies}.
  1214. \input{tables/efficacies.tex}
  1216. When carriers are trapped inside a quantum well a part of their energy has to be dissipated by phonon emission in order to match the carriers' energy to one of the available states inside a well.
  1217. \clearpage % *** manual new page 54/55
  1218. \noindent{}\marginpar{phosphor-converted led{\scriptsize{}s}}A white \led{} can be produced by using a phosphor down-conversion of a blue or \textsc{uv} light (pump) e.g. mixing blue InGaN light with yellow light from YAG:Ce phosphor. Phosphor has to convert blue photons ($\approx{}460
  1219. \:\textrm{nm}$) into yellow photons ($\approx{}570\:\textrm{nm}$) so their energy has to drop. It is impossible to do the down-conversion without losing a part of the energy~\cite{Chhajed2005b}. The Stokes energy loss is given by
  1220. \begin{equation} \label{eq:StokesShift}
  1221.         ΔE=hν_{1}-hν_{2}=\frac{hc}{λ_{1}} - \frac{hc}{λ_{2}}
  1222. \end{equation}
  1223. thus the efficiency of down-convertion is
  1224. \begin{equation} \label{eq:DownConvertionEfficiency}
  1225.         η=1-\frac{hν_{1}-hν_{2}}{hν_{1}} = \frac{λ_{1}}{λ_{2}}.
  1226. \end{equation}
  1228. \section[Lifetime]{lifetime}
  1229. In contrast to traditional, incandescent and fluorescent light sources, \led{}s do not tend to fail catastrophically. Instead, the output flux decreases with time and operating junction temperature. The degradation mechanisms were the subject for many lifetime studies. The cause and the magnitude of the lumen change differs depending on the structure and materials of the diode and its driving conditions: shape of the driving current (pulsed or direct) and the temperature. Solder join failures and change in thermal properties of the thermal interface materials used for mounting the diodes can also influence diodes' properties.
  1231. \marginpar{nitride diodes}For white \led{}s, a colour shift towards lower colour temperatures can be seen due to the yellowing of the die~\cite{Yang2010}. Studies show that mechanisms of \led{} structure and package deterioration like generation of threading dislocations and yellowing and cracking of encapsulating lenses are connected to the structure temperature~\cite{Narendran2004,Narendran2005,Yang2010}. The increase of the leakage current results in a decreasing photon generation and the deterioration of \led{} encapsulant lens lowers the structure extraction efficiency. In extreme cases, the defects may short the device rendering it unable to generate light~\cite{ZanoniPhDThesis}.
  1233. \marginpar{phosphide diodes}According to a 60\hsp000 hour study presented by Grillot et~al. the main cause of lumen degradation in AlGaInP diodes are the increase in defects concentration and in leakage current~\cite{Grillot2006}. The overall degradation, defined as the relative change of the luminous flux with respect to the flux at the beginning of the test, was found to behave according to equation
  1234. \begin{equation}
  1235.         D = D_{1}+D_{2}\:j+(D_{3}+D_{4}\:j)\textrm{ln}(t)
  1236. \end{equation}
  1237. where $D_{1-4}$ are diode dependent coefficients, $j$ is the current concentration and $t$ is the time.
  1238. \clearpage % *** manual new page 55/56
  1239. \noindent{}\marginpar{industry standards}As \led{} do not tend to fail catastrophically, instead their light output slowly decreases over time, the \led{} lifetime if defined as the time to some predefined intensity drop. Manufacturers use this measure together with driving current and junction temperature to quickly asses the lifetime of the devices they are using (fig.~\ref{fig:RebelReliabilityData}).
  1241. \marginpar{\textsc{b} and \textsc{l} lifetimes }
  1242. The \textsc{assist} alliance recommends using two metrics \textsc{b50} and \textsc{l70} to approximate the useful lifetime~\cite{ledLifeReport}. \textsc{b50} corresponds to the time where 50\thinspace{}\% of the device population fails. \textsc{l70} denotes the time in which \led{}s lose 30\thinspace{}\% of the initial luminous flux. Therefore combined \textsc{b50/l70} metric predicts when more than 50\thinspace{}\% of the diodes will drop below 70\thinspace{}\% of the initial flux~\cite{ReliabilityManual}. Depending on the application other sets can be used, e.g.~\textsc{b10/l70} or \textsc{b10/l70}.
  1244. \begin{figure}[!ht]% Luxeon Rebel B50/L70
  1245.         \centering
  1246.         \includegraphics{graphics/LED/Lifetime/B50_L70_luxeon_rebel.pdf}
  1247.         \caption{\textsc{b50/l70} lifetime predictions for InGaN (left) and AlGaInP (right) Luxeon Rebel diodes~\cite{RebelReliabilityData}.}
  1248.         \label{fig:RebelReliabilityData}
  1249. \end{figure}
  1251. %% LED dimming characteristics
  1252. \section[Dimming]{dimming}
  1253. Intensity of an \led{} can be controlled in many different ways. Because of the small dynamic resistance and the temperature dependency of forward voltage \led{}s are typically current controlled. All dimming techniques vary the mean current supplied to the diode changing the diode's intensity. Different current waveforms with the same average value can have different effect on \led{} output light. The peak wavelength can change creating a colour shift in the mixed light. Flux, and therefore efficiency, can change depending on a dimming method. Gu et~al.~\cite{Gu2006} report changes up to 100\% in efficiency depending on the dimming method, diode colour and the operating point.
  1255. Following chapters describe four different control schemes and analyse their effect on light-emitting diodes: two basic control schemes used commonly—pulse width modulation and amplitude modulation; a hybrid \am{}/\pwm{} dimming technique proposed by the author and a pulse code modulation technique.
  1257. %\marginpar{requirements}
  1258. %\marginpar{resolution}
  1260. % PWM
  1261. \subsection{Pulse width modulation}Because \led{}s are semiconductor-based devices they can be turned on and off very quickly. The minimum frequency at which pulsed light sources are perceived as constant is predicted by Ferry-Porter law
  1262. \begin{equation}
  1263.         \textit{critical flicker frequency} = k(\mathrm{log}\,L - \mathrm{log}\,L_{0})
  1264. \end{equation}
  1265. where $L$ is the luminous intensity, $L_0$ is the threshold intensity and $k$ is a constant having a typical value of 12\thinspace{}Hz/decade. The \textsc{cff} is furthermore dependent on the size of the stimuli and its position relative to the eye~\cite{Tyler1990}. Minimum modulation frequency value used in the industry is 60\thinspace{}Hz, but for high intensity, fast moving light sources this frequency should be increased to 300–1000Hz range to avoid flicker~\cite{Ashdown2006}.
  1267. Brightness control is achieved by modulating current pulse width \includegraphics{graphics/System/Dimming/PWM/sparklinePWM.pdf}. A duty cycle $d$ is a ratio of a on pulse time to the period of modulation frequency. Average luminous intensity is a product of the duty cycle and the constant luminous intensity during the on pulse.
  1268. \begin{equation}
  1269.         I_{f,avg}=d \cdot{} I_{f,max}
  1270. \end{equation}
  1271. Light output can be dimmed down to 0\% and linear control can be easily obtained. However, peak emission shifts and bandwidth narrowing occur during \pwm{} dimming and have to be taken into account when designing control strategy for light mixing~\cite{Dyble2005, Manninen2007}.
  1273. \begin{figure}[!ht]% PWM dimming linearity
  1274.         \centering
  1275.         \includegraphics{graphics/System/Dimming/PWM/linearityRGBW.pdf}
  1276.         \caption{Measured luminous flux for \pwm{} dimming of red, green, blue and white (left to right) diodes as a function of duty cycle. Data points correspond to the duty cycle equal to: 0, 0.2, 0.3, 0.6, 0.8 and 1.0. \pwm{} experiences almost ideal dimming linearity.}
  1277.         \label{fig:DimmingPWMLinearity}
  1278. \end{figure}
  1280. \noindent{}Under decreasing duty cycle the peak wavelengths shifts towards shorter wavelengths for AlGaInP~\cite{Gu2006, Manninen2007} and InGaN~\cite{Gu2006} diodes. The shift is explained by cooling of the active area in the chip with decreased duty cycle and, as a result, broadening of the bandgap energy~\cite{ledSchubert}.
  1282. \begin{figure}[!ht]% PWM dimming colour shifts
  1283.         \centering
  1284.         \includegraphics{graphics/System/Dimming/PWM/colourPointShiftRGBW.pdf}
  1285.         \caption{Measured colour shifts for the \pwm{} dimming of red, green, blue and white (left to right) diodes. Arrow shows the direction of decreasing flux. Data points correspond to the duty cycle equal to: 1.0 (marked with \cie{}~1976 colour coordinates), 0.8, 0.6, 0.4 and 0.2.}
  1286.         \label{fig:DimmingPWMColourShifts}
  1287. \end{figure}
  1289. \noindent{}Because of small colour shifts and very good linearity this dimming method is the most widely used in industrial applications. \led{} driver operates in a pulsed mode either with maximal forward current of the diode or no current. Therefore it can be easily optimised for a single operating point.
  1291. % AM
  1292. \subsection{Amplitude modulation}Amplitude modulation, sometimes referred to as continuous current reduction (\textsc{ccr}), uses variable \dc{} current to control the intensity of the diode. Because the efficiency of \led{}s increases with decreasing current concentration in the active area, the \am{} dimming scheme is inherently nonlinear. The increase of luminous efficiency over pulse width modulation can be as high as 100\% for low intensities for green InGaN diodes~\cite{Gu2006}.
  1294. \begin{figure}[!ht]% AM dimming linearity
  1295.         \centering
  1296.         \includegraphics{graphics/System/Dimming/AM/linearityRGBW.pdf}
  1297.         \caption{Measured luminous flux for \am{} dimming of red, green, blue and white (left to right) diodes as a function of forward current. Data points correspond to following currents: 0, 0.14, 0.28, 0.42, 0.56 and 0.7\thinspace{}\textsc{a}.}
  1298.         \label{fig:DimmingAMLinearity}
  1299. \end{figure}
  1301. \noindent{}The dimming range for amplitude modulation is typically limited to 10–100\% range. This is because at very low current levels the \textsc{snr} is low therefore the accuracy of the current control is limited.
  1303. Very big colour shifts occur when using \am{} in both pure colour diodes and phosphor converted white diodes~\cite{Dyble2005}. Typically, AlGaInP diodes experience peak wavelength shift towards shorter wavelengths and InGaN \led{}s towards longer wavelengths.
  1305. \begin{figure}[!ht]% AM dimming colour shifts
  1306.         \centering
  1307.         \includegraphics{graphics/System/Dimming/AM/colourPointShiftRGBW.pdf}
  1308.         \caption{Measured colour shifts for \am{} dimming of red, green, blue and white (left to right) diodes. Arrow shows the direction of decreasing flux. Data points correspond to following currents: 0.7 (marked with \cie{}~1976 colour coordinates), 0.56, 0.42, 0.28 and 0.14\thinspace{}\textsc{a}. Note the 10× increase of scale compared to fig.~\ref{fig:DimmingPWMColourShifts}.}
  1309.         \label{fig:DimmingAMColourShifts}
  1310. \end{figure}
  1312. \noindent{}Due to big spectral shifts and output flux to driving current nonlinearity, use of \am{} dimming technique in colour crucial applications is limited. Although if these problems would be solved, amplitude modulation would provide much more efficient dimming scheme than pulse width modulation.
  1314. \led{} driver has to provide power with varying voltage and current therefore the driver cannot be optimised for single operating point. On the other hand the luminous efficacy of the diode increases significantly with decreased current and the efficacy of the driver-\led{} system can be much higher than in \pwm{} dimmed luminaires.
  1316. % PWM/AM
  1317. \subsection{Hybrid modulation}\label{ssec:hybrid_pwm_am}Tse et~al.~\cite{Tse2009a,Tse2009b,Tse2009c,Tse2009d} proposed a general driving technique that is somewhere between the \pwm{} and \am{}. Their technique bases on the possibility to control the high and low peak current that is simultaneously dimmed with \pwm{}. Changing the peak currents would influence the efficiency of the \led{} and \pwm{} signal would control the brightness without big colour shifts. A variation of this technique was studied where low peak current is set to zero and the intensity of the diode is controlled by varying the peak current and duty cycle \includegraphics{graphics/System/Dimming/HybridPWMAM/sparklineHybridPWMAM.pdf}. The \led{} driver is off during $1-d$ part of the dimming period therefore losses in the converter are reduced. \led{} is controlled by varying both duty cycle and peak current.
  1318. \begin{equation}
  1319.         I_{f,avg} = d \cdot{} I_{f,peak}
  1320. \end{equation}
  1321. It is always possible to roll back to either \pwm{} or \am{} control by setting either current or duty cycle, respectively, to their maximum value.
  1323. \begin{figure}[!ht]% PWM/AM dimming XYZ
  1324.         \centering
  1325.         \includegraphics{graphics/System/Dimming/HybridPWMAM/XYZ.pdf}
  1326.         \caption{Measured, normalised tristimulus values for \pwm{}/\am{} hybrid dimming as a function of relative peak current and duty cycle. Data gathered at constant heatsink temperature.}
  1327.         \label{fig:DimmingPWMAM_XYZ}
  1328. \end{figure}
  1330. Hybrid \pwm{}/\am{} modulation creates a degree of freedom in diode control. This gives an opportunity to optimise colorimetric and radiometric properties of an \led{}.
  1332. \noindent{}\marginpar{stabilising white led colour point}Noticing the fact that for white phosphor converted diode \pwm{} and \am{} dimming methods create opposite colour shifts (figures \ref{fig:DimmingPWMColourShifts} and \ref{fig:DimmingAMColourShifts}) and the direction of this shift is almost parallel to the shift caused by the temperature change (fig ref needed) one can use hybrid \am{}/\pwm{} to stabilise both colour point of a diode and luminous flux over a range of temperatures.
  1334. \begin{figure}[!ht]% PWM/AM dimming stabilizing white diode colour point
  1335.         \centering
  1336.         \includegraphics{graphics/LED/Dimming/HybridPWMAM/HybridDimmingWhiteDiodeTemperatureShift.pdf}
  1337.         \caption{Colour shifts of white \textsc{pc} \led{} caused by \pwm{} and \am{} dimming mechanisms and the directions of shift caused by increasing the heatsink temperatures (left). Relative driving current and duty cycle for obtaining stable colour point.}
  1338.         \label{fig:DimmingPWMAMStabilizingWhiteColourPoint}
  1339. \end{figure}
  1341. \noindent{}Figure \ref{fig:DimmingPWMAMStabilizingWhiteColourPoint} shows the directions of colour shift caused by two dimming schemes and the shift caused by the changing heatsink temperature. Shift caused by the temperature can be subdivided into two vectors: parallel and perpendicular to the shifts caused by dimming. The perpendicular part cannot be controlled using hybrid dimming technique but is very small in magnitude.
  1343. An experiment was set up to prove the hypothesis. A white, phosphor converted diode was placed in the test adapter. The heatsink temperature was set to 20°\textsc{c} and the diode was energised with nominal current. Colour coordinates and luminous flux were recorded when the diode reached thermal steady state. After increasing the heatsink temperature to 30°\textsc{c} the duty cycle and driving current were manually adjusted to bring the colour point as close to the starting colour point without changing the flux value. It was only possible if the current was increased over the nominal current value. Although the average current value did not exceed the nominal driving current diode's lifetime may be affected, which may not be recommended in some applications. In order not to exceed the maximum forward current, the experiment was repeated with lower initial driving current equal to around 60\thinspace{}\% of nominal current (which corresponds to ~70\thinspace{}\% of the nominal flux value). The flux was kept constant in the heatsink temperature range of 20–48°\textsc{c} while maintaining very low colour shift without overdriving the diode.
  1345. Decreasing the initial current (at $d=1$) or allowing for overdriving increases the range at which the colour point is kept stable.
  1347. \begin{table}[!ht]\footnotesize
  1348.         \caption{Current and voltage values needed to obtain stable colour point with the change of the heatsink temperature. Colour coordinates and the distance from initial colour point given in the \cie{}~1976 colour space.}
  1349.         \label{tab:experimentDataStableWhiteColourPointUnderHybridDimming}
  1350.         \centering
  1351.         \begin{tabular}{ccccccc}
  1352.                 $T_{hs} [^{\circ}\textrm{C}]$ & \textit{d} [p.u.] & $I_{peak}$ [A] & \textit{F} [lm] & $u'$ & $v'$ & $Δu'v'$ \\
  1353.                 \hline
  1354.                 20 & 1.000 & 0.411 & 53.3 & 0.2085 & 0.4981 & 0.00000 \\       
  1355.                 30 & 0.905 & 0.485 & 53.3 & 0.2085 & 0.4981 & 0.00000 \\       
  1356.                 40 & 0.796 & 0.597 & 53.3 & 0.2084 & 0.4982 & 0.00011 \\       
  1357.                 50 & 0.716 & 0.712 & 53.3 & 0.2083 & 0.4982 & 0.00021 \\       
  1358.                 \hline
  1359.         \end{tabular}
  1360. \end{table}
  1362. \noindent{}Measured data points with the colour distance from the initial colour point are presented in table~\ref{tab:experimentDataStableWhiteColourPointUnderHybridDimming}. Plotting the required current and duty cycle values (fig.~\ref{fig:DimmingPWMAMStabilizingWhiteColourPoint}) reveals that the relationship between the controlled values and heatsink temperature is almost linear and can be expressed using following formulas:
  1363. \begin{align}
  1364.         I_{f,peak} &= 0.0001 \cdot{} T_{hs}^{2} + 0.0030 \cdot{} T_{hs} + 0.3088 \nonumber \\
  1365.         d &= 0.0004 \cdot{} T_{hs}^{2} - 0.0122 \cdot{} T_{hs} + 1.2309
  1366. \end{align}
  1368. \noindent{}\marginpar{stabilising peak wavelength position}Another application of hybrid dimming can be at stabilising peak wavelength position during dimming. Colour sensors used in luminaires as a colour feedback are very sensitive to the change of spectral shape~\cite{Ashdown2007}. This phenomenon is caused by the very narrow shape of \led{} spectra. The spectra can shift towards shorter and longer wavelengths during dimming and as a result of heatsink temperature change. If one of this mechanisms of colour shifts could be eliminated, the accuracy of the colour sensors could in theory be improved.
  1370. \begin{figure}[!ht]% PWM/AM dimming peak wavelength shift
  1371.         \centering
  1372.         \includegraphics{graphics/LED/Dimming/HybridPWMAM/HybridDimmingPeakWavelengthShift.pdf}
  1373.         \caption{Peak wavelength shift with respect to the nominal position at 20°\textsc{c} for green (left) and blue (right) diodes under different dimming techniques.}
  1374.         \label{fig:DimmingPWMAMPeakWavelengthShift}
  1375. \end{figure}
  1377. InGaN diodes exhibit opposite peak wavelength shifts under \pwm{} and \am{}. A~hybrid modulation could therefore fix the peak wavelength position (fig.~\ref{fig:DimmingPWMAMPeakWavelengthShift}).
  1379. The initial peak wavelength position was recorded at nominal driving conditions at 20°\textsc{c} heatsink temperature. The duty cycle was then lowered by 10\thinspace{}\% and, by adjusting the peak current, the peak wavelength was moved to initial value. The process was repeated for duty cycles down to 20\thinspace{}\%. Pairs of peak current and duty cycle that satisfied the constant peak wavelength conditions were recorded also for 40 and 60°\textsc{c} to examine if the relationship is temperature dependent.
  1381. \begin{figure}[!ht]% PWM/AM dimming Ipeak vs d
  1382.         \centering
  1383.         \includegraphics{graphics/LED/Dimming/HybridPWMAM/HybridDimmingPeakWavelengthShiftPeakCurrentVsDuty.pdf}
  1384.         \caption{Driving conditions for green (left) and blue (right) diodes necessary to obtain stable peak wavelength position as shown in figure~\ref{fig:DimmingPWMAMPeakWavelengthShift}.}
  1385.         \label{fig:DimmingPWMAMCPeakCurrentVsDuty}
  1386. \end{figure}
  1388. \noindent{}Results (fig.~\ref{fig:DimmingPWMAMCPeakCurrentVsDuty}) show that the relationship between peak current and duty cycle needed to obtain stable peak wavelength is not dependent on heatsink temperature and can be described be the following equations
  1389. \begin{align}
  1390.         I_{f,peak,green} &= 0.1739 \cdot{} d^{2} + 0.0079 \cdot{} d + 0.5109 \nonumber \\
  1391.         I_{f,peak,blue}  &= 0.0410 \cdot{} d^{2} + 0.0727 \cdot{} d + 0.5851
  1392. \end{align}
  1394. \begin{figure}[!ht]% PWM/AM dimming colour point shift
  1395.         \centering
  1396.         \includegraphics{graphics/LED/Dimming/HybridPWMAM/HybridDimmingColourPointShift.pdf}
  1397.         \caption{Colour point shift for green (left) and blue (right) diodes under different dimming techniques. Heatsink temperature fixed to 20°\textsc{c}.}
  1398.         \label{fig:DimmingPWMAMColourPointShift}
  1399. \end{figure}
  1401. \noindent{}For both green and blue InGaN diodes colour point resulting from hybrid dimming moves towards the border of \cie{}~1931 diagram thus the colour of the light becomes purer. As seen in figure~\ref{fig:DimmingPWMAMColourPointShift} the colour shift is confined to a one step MacAdam ellipse.
  1403. % PCM
  1404. \subsection{Pulse code modulation}A pulse code modulation uses a current pattern that resembles the binary representation of the desired dimming value, e.g.~0.33 can be written as 01010100$_b$ and the corresponding current pattern looks like this \includegraphics{graphics/System/Dimming/PCM/sparklinePCM.pdf}. The dimming value is therefore encoded in the light pattern and can be decoded by other devices~\cite{Howell2002}. The first pulse in the pattern has a duration of $T$ and corresponds to the most important bit of the binary dimming value. The next pulse is two times shorter and corresponds to the next bit. Last pulse has the width of $T/2^{n-1}$, where $n$ is the resolution of the dimming signal. The total pulse pattern period has to be lower than 300\thinspace{}Hz to avoid visual flicker and therefore the dimming resolution is limited by the timing properties of the microcontroller and the length of the shortest pulse. Ashdown~\cite{Ashdown2006} proposed the use of temporal dithering in order to increase the dimming resolution in \pcm{}.
  1406. \begin{figure}[!ht]% PCM spectrum
  1407.         \centering
  1408.         \includegraphics{graphics/System/Dimming/PCM/pcm_spectrum.pdf}
  1409.         \caption{Spectra of \pwm{} (left) and \pcm{} (right) current shape. Base frequency, marked with red bar, is 200\thinspace{}Hz in both cases. The \dc{} value for both signals is the same. The \pcm{} signal contains more high frequency components in the spectrum.}
  1410.         \label{fig:PCMSpectrum}
  1411. \end{figure}
  1413. \noindent{}The biggest advantage of \pcm{} over \pwm{} is that the driving pattern can be generated for many channels using only a single timer in a microcontroller thus there is no need for hardware \pwm{} generators. The current spectra contains more high frequency components than under a \pwm{} driving scheme (figure~\ref{fig:PCMSpectrum}) therefore this dimming method is mostly suitable for driving high number of low power diodes like in a matrix displays.
  1415. % current-voltage model
  1416. \section[Current-voltage model]{current-voltage model}\label{sec:iv_model}
  1417. Recalling the fact that forward voltage of a diode is proportional to both junction temperature and forward current and diode colorimetric properties are dependent on junction temperature and forward current one can construct a model that estimates colorimetric properties basing only on instantaneous values of diode's current and voltage.
  1419. Diodes' tristimulus properties and forward voltage have to be measured at various junction temperatures and, if \am{} or hybrid dimming is to be used, various current levels in order to create a model.
  1421. \marginpar{data acquisition}A green InGaN test diode was placed on thermally controlled heatsink. Temperature was set in increments of 10°\textsc{c} in the 5–55°\textsc{c} range. The current was controlled in the 10–100\% of nominal current range allowing the diode to reach thermal steady state after each change. At that point current, voltage and tristimulus values were measured. Four wire setup was used to gather electrical parameters in order to avoid voltage drop on terminals and cables.
  1423. \noindent{}\marginpar{surface fitting}Collected data was fitted to a polynomial. Three separate functions were created for each of the tristimulus values and one function was created for the radiometric power. The higher order of polynomial used to represent the data the more computational power is required to use the model. Also the model may be overfitted in which case the measurement error would be visible in the model. By a trial and error approach the lowest order of polynomial still representing measured data accurately was found to be in a form of:
  1424. \begin{equation}
  1425.         a_{0} + a_{1}I + a_{2}I\,^{2} + a_{3}I\,^{3} + a_{4}V + a_{5}V\,^{2}
  1426. \end{equation}
  1427. %Model may be stored as a lookup table instead of mathematical function.
  1429. \begin{table}[!ht]\footnotesize
  1430. \caption{Model coefficients $a_{n}$ for the measured diode model. \label{tab:IV2XYZ_model_coefficients}}
  1431. \centering
  1432. \begin{tabular}{lccccccc}
  1433. &  $a_{0}$  &   $a_{1}$ &   $a_{2}$ &  $a_{3}$ &  $a_{4}$ &  $a_{5}$ & $R^{2}$\\
  1434.         \hline
  1435. \textsc{\textit{x}} &           -21.5417 & \phantom{0}49.8094 & \phantom{0}-38.8920 &           16.1661 & \phantom{00}17.3677 &           -3.4374 & 0.99997 \\
  1436. \textsc{\textit{y}} & \phantom{0}30.7711 &           150.6126 &           -119.9930 &           45.0165 & \phantom{0}-38.8748 &           11.2224 & 0.99991 \\
  1437. \textsc{\textit{z}} & \phantom{0}33.1338 & \phantom{0}14.2156 & \phantom{00}-6.8995 & \phantom{0}2.0560 & \phantom{0}-28.7647 & \phantom{0}6.2476 & 0.99986 \\
  1438.         $P_{rad}$  &           120.6807 &           189.0981 &           -145.9986 &           52.2570 &           -118.7007 &           28.9492 & 0.99989 \\
  1439. \hline
  1440. \end{tabular}
  1441. \end{table}
  1443. \begin{figure}[!ht]% IV2XYZ model
  1444.         \centering
  1445.         \includegraphics{graphics/LED/IV2XYZmodel/IV2XYZmodel.pdf}
  1446.         \caption{Model surfaces of $X$, $Y$ and $Z$ tristimulus values and optical power of the measured \led{}. Model is only valid in the vicinity to the measured data points. }
  1447.         \label{fig:IV2XYZmodel}
  1448. \end{figure}
  1450. \begin{figure}[!ht]% model verification. DC current
  1451.         \centering
  1452.         \includegraphics{graphics/LED/IV2XYZmodel/DCaccuracy.pdf}
  1453.         \caption{Measured ($\circ$) and modelled ($\cdot$) colour points on \textsc{cie1976} chromaticity diagram for various combinations of driving current and heatsink temperature.}
  1454.         \label{fig:IV2XYZmodelVerificationDC}
  1455. \end{figure}
  1457. \marginpar{model verification at dc current}In order to test the model and the data fitting, six measurement points were taken at random heatsink temperatures and driving currents. Measured current, voltage and tristimulus values are compared to estimated tristimulus values in table~\ref{tab:IV2XYZ_model_verification_DC} and the corresponding colour points for the measured and modelled data are shown in figure~\ref{fig:IV2XYZmodelVerificationDC}. For all measured data points the colour distance between the measured and estimated colour points is much lower than the just noticeable colour distance. This implies that the model predicts the chromaticity and intensity of \led{}s with high enough accuracy.
  1459. \begin{table}[!ht]\footnotesize
  1460. \caption{Measured current, voltage and tristimulus values compared to corresponding modelled tristimuls values. Error in model is visible as colour distance between estimated and measured colour points.\label{tab:IV2XYZ_model_verification_DC}}
  1461. \centering
  1462. \begin{tabular}{rrrrrrrrrrrrr}
  1463. & & & \multicolumn{3}{c}{\textsc{measured}} & & \multicolumn{3}{c}{\textsc{modelled}} & & & \\
  1464. \multicolumn{1}{c}{\textit{v}} & \multicolumn{1}{c}{\textit{i}} & & \multicolumn{1}{c}{\textsc{\textit{x}}} & \multicolumn{1}{c}{\textsc{\textit{y}}} & \multicolumn{1}{c}{\textsc{\textit{z}}} & & \multicolumn{1}{c}{\textsc{\textit{x}}} & \multicolumn{1}{c}{\textsc{\textit{y}}} & \multicolumn{1}{c}{\textsc{\textit{z}}} & & \multicolumn{2}{c}{$ΔE_{ab}^{*}$} \\
  1465. \cline{1-2}\cline{4-6}\cline{8-10}\cline{12-13}
  1466. 3.1388 & 0.5253 & & 16.917 & 71.915 & 10.266 & & 16.883 & 71.847 & 10.261 & & \multicolumn{2}{c}{0.0049} \\
  1467. 2.9731 & 0.3523 & & 13.142 & 54.639 &  7.097 & & 13.138 & 54.528 &  7.080 & & \multicolumn{2}{c}{0.0661} \\
  1468. 3.0566 & 0.4648 & & 15.820 & 65.558 &  8.933 & & 15.801 & 65.394 &  8.904 & & \multicolumn{2}{c}{0.0738} \\
  1469. 3.1348 & 0.6151 & & 18.846 & 76.887 & 10.942 & & 18.809 & 76.911 & 10.970 & & \multicolumn{2}{c}{0.0596} \\
  1470. 3.1809 & 0.6949 & & 20.216 & 82.294 & 12.013 & & 20.180 & 82.488 & 12.087 & & \multicolumn{2}{c}{0.3098} \\
  1471. 2.8123 & 0.2355 & &  9.935 & 39.530 &  4.628 & &  9.899 & 39.604 &  4.643 & & \multicolumn{2}{c}{0.2676} \\
  1472. \hline
  1473. \end{tabular}
  1474. \end{table}
  1476. \begin{figure}[!ht]% model verification. IV PWM waveforms
  1477.         \centering
  1478.         \includegraphics{graphics/LED/IV2XYZmodel/PWMwaveforms.pdf}
  1479.         \caption{Measured current and voltage waveforms for the modelled diode at various duty cycles.}
  1480.         \label{fig:IV2XYZmodelVerificationPWMwaveforms}
  1481. \end{figure}
  1483. \begin{figure}[!ht]% model verification. IV PWM waveforms
  1484.         \centering
  1485.         \includegraphics{graphics/LED/IV2XYZmodel/PWMaccuracy.pdf}
  1486.         \caption{Difference between measured colour shift for \pwm{} dimming (dashed line) and colour shift obtained using current-voltage model (solid line). Arrow shows the direction of decreasing duty cycle. Colour point at nominal driving current ($d=1$) marked with $u'v'$ coordinates. }
  1487.         \label{fig:IV2XYZmodelVerificationPWM}
  1488. \end{figure}
  1490. \noindent{}\marginpar{model verification at pulsed current}To test the accuracy of the model the test diode was driven with 200\thinspace{}Hz pulsed current with various duty cycles (fig.~\ref{fig:IV2XYZmodelVerificationPWMwaveforms}) at constant heatsink temperature. Before each measurement the diode was allowed to reach thermal steady state. Optical parameters were integrated over multiple \pwm{} periods. Current and voltage waveforms were recorded with 250\thinspace{}kS/s speed with 12~bit \textsc{adc}. Few periods were extracted from the waveforms and instantaneous current and voltage values were converted to tristimulus values using the previously created model. Resulting data was integrated and divided by the measured time period to obtain average tristimulus values and the resulting colour point.
  1492. \noindent{}Results, summarised in table~\ref{tab:IV2XYZ_model_verification}, show that the model accurately predicts the colour shift of the diode. It is therefore possible to use this model, as a colour feedback regardless of the dimming scheme used in the luminaire. The accuracy of this method lower with lowering duty cycle. This may be attributed to fast thermal and electrical transients that tend to dominate at very low duty cycle values. In order to measure these transients precisely, much higher voltage and current sampling would be required. Also lower number of samples are used to estimate the colorimetric properties therefore noise in the measurement may influence the accuracy of the model.
  1494. \begin{table}[!ht]\footnotesize
  1495. \caption{Measured and modelled tristimulus values of the test diode and the colour distance between the corresponding colour points. \label{tab:IV2XYZ_model_verification}}
  1496. \centering
  1497. \begin{tabular}{llrrrrrrrrrrr}
  1498. & & & \multicolumn{3}{c}{\textsc{measured}} & & \multicolumn{3}{c}{\textsc{modelled}} & & \\
  1499. \multicolumn{2}{c}{\textit{d}} & & \multicolumn{1}{c}{\textsc{\textit{x}}} & \multicolumn{1}{c}{\textsc{\textit{y}}} & \multicolumn{1}{c}{\textsc{\textit{z}}} & & \multicolumn{1}{c}{\textsc{\textit{x}}} & \multicolumn{1}{c}{\textsc{\textit{y}}} & \multicolumn{1}{c}{\textsc{\textit{z}}} & & \multicolumn{2}{c}{$ΔE_{ab}^{*}$} \\
  1500. \cline{1-2}\cline{4-6}\cline{8-10}\cline{12-13}
  1501. \multicolumn{2}{l}{0.2} & &  4.0447 & 17.3423 &  2.6773 & &  3.8695 & 16.6655 &  2.5735 & & \multicolumn{2}{c}{2.4091} \\
  1502. \multicolumn{2}{l}{0.4} & &  8.1630 & 34.4009 &  5.2286 & &  8.1669 & 34.6032 &  5.2597 & & \multicolumn{2}{c}{0.4614} \\
  1503. \multicolumn{2}{l}{0.6} & & 12.5078 & 51.9626 &  7.8009 & & 12.5364 & 52.2906 &  7.8480 & & \multicolumn{2}{c}{0.4869} \\
  1504. \multicolumn{2}{l}{0.8} & & 16.8984 & 69.2601 & 10.2805 & & 16.9262 & 69.5062 & 10.3069 & & \multicolumn{2}{c}{0.1671} \\
  1505. \multicolumn{2}{l}{1.0} & & 20.7419 & 83.9519 & 12.3199 & & 20.7569 & 84.0777 & 12.3363 & & \multicolumn{2}{c}{0.0319} \\
  1506. \hline
  1507. \end{tabular}
  1508. \end{table}
  1510. \noindent{}\marginpar{comments}The number of measurement points needed to describe diode's behaviour depends on the complexity of the model and the desired operating point. If the diode is to be dimmed using \pwm{} scheme, it is sufficient to measure the parameters at single forward current with varying heatsink temperature.
  1512. If the radiometric power of the emitted light is also modelled, it may be used together with the thermal model of the system to completely describe luminaire's behaviour. Subtracting the optical power from input electrical power gives the power losses in the diode structure. This loss will create a temperature rise with respect to heatsink temperature according to the thermal description of the heat flow path \ref{fig:RC_network_thermal_model}. Resulting temperature increase in the junction can be translated into forward voltage change resulting a new operating point to be fed back to the current-voltage model.
  1514. \marginpar{modelling led strings}Very often luminaire consists of strings of series connected light-emitting diodes. It is possible to model each diode in the string measuring separately individual diode voltages and string current but this approach is not practical. A single model created for the whole string would therefore be beneficial. Some restrictions in thermal design exist when using the current-voltage model with multiple \led{}s. Model is created with only two external stimuli changing: current and heatsink temperature.
  1516. % chapter light_emitting_diodes (end)
  1518. %%%%%%%%%%%%%%%
  1519. %% Chapter 4 %%
  1520. %%%%%%%%%%%%%%%
  1521. \chapter{Luminaire control}% (fold)
  1522. \label{cha:luminaire_control}
  1523. Polychromatic solid-state luminaires need a control system in order to keep the desired luminance and colour level stable. Diodes' parameters will change during normal operation due to self heating, ambient temperature variations and during ageing and these changes should be compensated by a proper control mechanism.
  1525. Firstly, the open loop control is discussed. The methods of choosing an appropriate operating point, depending on the number of primaries used in the luminaire, are given. For luminaires containing more than three primaries optimisation techniques may be used to optimise some of the lamp parameters like luminous flux or efficacy.
  1527. A review of existing colour control methods is presented, where the main control mechanisms: temperature feed forward, flux feedback and colour coordinates feedback are discussed. A colour control loop basing on the current-voltage diode model, presented in chapter~\ref{sec:iv_model} is shown.
  1529. Various optimisation possibilities are discussed. From optimally choosing a working point of a polychromatic luminaire to optimisation benefits for trichromatic luminaire using the hybrid dimming mechanism, presented in chapter~\ref{ssec:hybrid_pwm_am}.
  1531. %% Luminaire control strategies
  1532. \section[Colour control strategies]{colour control strategies}
  1534. Control scheme will differ depending on the number of available colours in luminaire. For polychromatic luminaires (three or more basic colours) control of colour point and luminance is possible by adjusting intensities of individual diodes. In trichromatic luminaires there is only one solution for the colour equations that yields the desired colour point. If the luminaire consists of four or more diodes the number of possible solutions may be infinite therefore the lamp operation can be optimised by choosing the operating point.
  1536. \subsection{Open loop}
  1537. Open loop system uses only the calibration matrix~\ref{eq:calibrationMatrix3Diodes} to calculate required duty cycle. The accuracy of this control method is limited by the self-heating of the diodes and the changes of ambient temperature because the values in calibration matrix are measured at a single temperature value. However this type of control can still be used in luminaires where colour accuracy is not important. Because the system is modelled as a set of linear equations and linear algebra principles are used to find the operating point the \pwm{} dimming mechanism should be used in the luminaires controlled with this method.
  1539. \begin{figure}[!ht]% Open loop
  1540.         \centering
  1541.         \includegraphics{graphics/System/open_loop.pdf}
  1542.         \caption{Open loop luminaire control system.}
  1543.         \label{fig:SystemOL}
  1544. \end{figure}
  1546. \noindent{}The optical characteristics of a trichromatic luminaire are characterised by a calibration matrix $\mathbi{C}$ containing diodes' tristimulus values at nominal conditions.
  1547. \begin{equation} \label{eq:calibrationMatrix3Diodes}
  1548. \mathbi{C} =
  1549. \left[ \begin{array}{ccc}
  1550. X_{1} & X_{2} & X_{3} \\
  1551. Y_{1} & Y_{2} & Y_{3} \\
  1552. Z_{1} & Z_{2} & Z_{3}
  1553. \end{array} \right]
  1554. \end{equation}
  1555. The calibration matrix is strongly temperature dependent as the diode experience peak wavelength shift and flux change with the change of the temperature.
  1557. Duty cycle needed to obtain desired color point and luminance level can be calculated by using inverse of the calibration matrix.
  1558. \begin{equation} \label{eq:openLoopColor}
  1559. \mathbi{d} = \left[ \begin{array}{c}d_{1} \\ d_{2} \\ d_{3}\end{array} \right]   =   \mathbi{C}^{-1} \cdot{} \!\left[ \begin{array}{c}X \\ Y \\ Z \end{array} \right]
  1560. \end{equation}
  1561. If any of the calculated duty cycle values is negative then the desired colour point lays outside the gamut of the device. If, on the other hand, any of the values is above~1 then the desired luminance level cannot be reached without overdriving the diode.
  1563. \marginpar{underdetermined system}When using more than three different diode colours, the system becomes underdetermined as there are three equations for each of the tristimuli value and $n>3$ unknowns—the diodes' relative output. The calibration matrix expands to a $3 \times{} n$ matrix composed of the diodes' tristimulus values.
  1564. \begin{equation} \label{eq:calibrationMatrixNDiodes}
  1565. \mathbi{C} =
  1566. \left[ \begin{array}{cccc}
  1567. X_{1} & X_{2} & \cdots & X_{n} \\
  1568. Y_{1} & Y_{2} & \cdots & Y_{n}\\
  1569. Z_{1} & Z_{2} & \cdots & Z_{n}
  1570. \end{array} \right]
  1571. \end{equation}
  1572. For an underdetermined system, there are typically infinitely many solutions of the problem. Every solution to underdetermined will  have the form $\mathbi{d}+c \cdot \mathbi{s}$ where $\mathbi{x}$ is a particular solution of $\mathbi{Cx}=\mathbi{b}$ and $c \cdot \mathbi{s}$ is a linear combination of solution to $\mathbi{Cs}=\mathbi{0}$.
  1573. \marginpar{row reduction}As an example a system composed of four diodes: red, green, blue and white, defined by calibration matrix
  1574. \begin{equation} \label{eq:calibrationMatrix4DiodesExample}
  1575. \mathbi{C} =
  1576. \left[ \begin{array}{cccc}
  1577. 60      & 6             & 17    & 22    \\
  1578. 25      & 26    & 5             & 20    \\
  1579. 0               & 4             & 100   & 20
  1580. \end{array} \right]
  1581. \end{equation}
  1582. generating $X=20$, $Y=40$ and $Z=20$ will be investigated. In order to solve the system equations an augmented matrix $(\mathbi{C}|\mathbi{b})$ is composed.
  1583. \begin{equation} \label{eq:augmentedMatrixNDiodes}
  1584. \left(\mathbi{C}|\mathbi{b}\right) =
  1585. \left[ \begin{array}{cccc|c}
  1586. 60      & 6             & 17    & 22 & 20       \\
  1587. 25      & 26    & 5             & 20 & 40       \\
  1588. 0               & 4             & 100   & 20 & 20
  1589. \end{array} \right]
  1590. \end{equation}
  1591. Using Gauss-Jordan elimination matrix \ref{eq:augmentedMatrixNDiodes} is transformed into following form
  1592. \begin{equation} \label{eq:augmentedMatrixGaussJordanEliminationNDiodes}
  1593. \left(\mathbi{C}|\mathbi{b}\right) =
  1594. \left[ \begin{array}{cccc|c}
  1595. 1 & 0 & 0 & 0.2677 & 0.1560 \\
  1596. 0 & 1 & 0 & 0.4770 & 1.3604 \\
  1597. 0 & 0 & 1 & 0.1809 & 0.1456
  1598. \end{array} \right]
  1599. \end{equation}
  1600. The corresponding system equations are given by
  1601. \begin{equation} \label{eq:equationsGaussJordanEliminationNDiodes}
  1602. \begin{array}{rcl}
  1603. d_{1} + 0.2677 \cdot d_{4} & = & 0.1560 \\
  1604. d_{2} + 0.4770 \cdot d_{4} & = & 1.3604 \\
  1605. d_{3} + 0.1809 \cdot d_{4} & = & 0.1456
  1606. \end{array}
  1607. \end{equation}
  1608. Introducing $c = d_{4}$ and rewriting equations in matrix form yields
  1609. \begin{equation} \label{eq:matrixGaussJordanEliminationNDiodes}
  1610. \mathbi{d} = \mathbi{d}_{re} + c \cdot \mathbi{s}_{re} =
  1611. \left[ \begin{array}{c}
  1612. d_{1} \\ d_{2} \\ d_{3} \\ d_{4}
  1613. \end{array} \right] =
  1614. \left[ \begin{array}{c}
  1615. 0.1560 \\ 1.3604 \\  0.1456 \\ 0
  1616. \end{array} \right] + c
  1617. \left[ \begin{array}{c}
  1618. -0.2677 \\ -0.4770 \\ -0.1809 \\ 1
  1619. \end{array} \right]
  1620. \end{equation}
  1621. Green diode in the initial solution $\mathbi{d}_{re}$ has to be overdriven by 36\%, therefore by varying the value of $c$ another solution has to be found.
  1624. \marginpar{pseudoinversion}By analogy to equation~\ref{eq:openLoopColor} solution to the underdetermined system can also be calculated using an inverted matrix. As inversion of the non-square matrix is not possible therefore in order to find an algebraic solution to the problem a pseudoinverted matrix $\mathbi{C}^{+}$ can be used~\cite{Greville1959}.
  1625. \begin{equation} \label{eq:openLoopColorNDiodes}
  1626. \mathbi{d}_{pinv} =\left[ \begin{array}{c}d_{1} \\ d_{2} \\ \vdots \\ d_{n} \end{array} \right] = \mathbi{C}^{+} \cdot{} \! \left[ \begin{array}{ccc}X \\ Y \\ Z \end{array} \right]
  1627. \end{equation}
  1628. Pseudoinversion yields the solution that has the smallest Euclidean norm $||d||_{2}$. This solution may or may not lay within the feasible solution space, where all duty cycles are in range of $0 \leq d_{n} \leq 1$. Luminaire described by calibration matrix \ref{eq:calibrationMatrix4DiodesExample} and pseudoinverted matrix
  1629. \begin{equation} \label{eq:pseudoinvertedCalibrationMatrix4DiodesExample}
  1630. \mathbi{C}^{+} =
  1631. \left[ \begin{array}{rrr}
  1632.    0.0189 &  -0.0076 &  -0.0032 \\
  1633.   -0.0164 &   0.0356 &   0.0004 \\
  1634.    0.0012 &  -0.0043 &   0.0098 \\
  1635.   -0.0026 &   0.0142 &   0.0011
  1636. \end{array} \right]
  1637. \end{equation}
  1638. generating $X=20$, $Y=40$ and $Z=20$, according to eq.~\ref{eq:openLoopColorNDiodes} should be driven as follows
  1639. \begin{equation} \label{eq:openLoop4DiodesExample}
  1640. \mathbi{d}_{pinv} =
  1641. \mathbi{C}^{+} \cdot{}\!
  1642. \left[ \begin{array}{r}
  1643.    20 \\
  1644.    40 \\
  1645.    20
  1646. \end{array} \right] =
  1647. \left[ \begin{array}{r}
  1648.    0.0226 \\
  1649.    1.0350 \\
  1650.    0.0413 \\
  1651.    0.5866
  1652. \end{array} \right]
  1653. \end{equation}
  1654. This particular solution lays outside available duty cycle values. In order to find other possible solutions a linear combination of the basis vectors for the kernel (null space) for calibration matrix $\mathbi{C}$ can be added to solution vector $\mathbi{d}$. The null space for matrix $\mathbi{C}$ is defined as set of vectors $\mathbi{s}$ for which $\mathbi{C}\cdot{}\mathbi{s}=\mathbi{0}$
  1655. \begin{equation} \label{eq:nullSpace4DiodesExample}
  1656. \mathbi{s} =
  1657. \left[ \begin{array}{r}
  1658.   -0.1929 \\
  1659.   -0.4704 \\
  1660.   -0.1508 \\
  1661.    0.8478
  1662. \end{array} \right]
  1663. \end{equation}
  1664. A feasible solution can be calculated by adding a linear combination of the basis vector $\mathbi{s}$. In this case e.g. $c=1/10$ will yield a feasible solution:
  1665. \begin{equation} \label{eq:openLoop4DiodesExamplePlusVector}
  1666. \mathbi{d} = \mathbi{d}_{pinv}+c \cdot{}\mathbi{s} =
  1667. \left[ \begin{array}{r}
  1668.    0.0226 \\
  1669.    1.0350 \\
  1670.    0.0413 \\
  1671.    0.5866
  1672. \end{array} \right] + \frac{1}{10}
  1673. \left[ \begin{array}{r}
  1674.   -0.1929 \\
  1675.   -0.4704 \\
  1676.   -0.1508 \\
  1677.    0.8478
  1678. \end{array} \right] =
  1679. \left[ \begin{array}{r}
  1680.    0.0033 \\
  1681.    0.9880 \\
  1682.    0.0262 \\
  1683.    0.6714
  1684. \end{array} \right]
  1685. \end{equation}
  1687. \noindent{}In general, finding a feasible solution by means of algebraic calculation will require finding an initial point and moving to the feasible range.
  1689. \begin{figure}[ht]% Graphical solution for 4 diode underdetermined system
  1690.         \centering
  1691.         \includegraphics{graphics/System/graphical_solution.pdf}
  1692.         \caption{Graphical solution for four diode underdetermined system. Feasible solution range is a subset of all solutions where all relative intensities are in the range of [0,1]. Initial solutions $\mathbi{d}_{re}$  and $\mathbi{d}_{pinv}$ described by equations \ref{eq:matrixGaussJordanEliminationNDiodes} and \ref{eq:openLoop4DiodesExample}, respectively, shown in the solution space.}
  1693.         \label{fig:openLoop4DiodesExampleGraphicalSolution}
  1694. \end{figure}
  1695. \marginpar{conclusions} Both techniques described above, pseudoinversion and row reduction, find the solution in form of initial solution plus a way to move the solution within the solution space. Graphical representation (fig.~\ref{fig:openLoop4DiodesExampleGraphicalSolution}) of the solution space shows for system described by calibration matrix \ref{eq:calibrationMatrix4DiodesExample} shows both initial solutions. The advantage of using row elimination is that initial solution will always yield a vector with $3-n$ zeroes and searching through solution space is limited by searching within $c_{re} \in [0,1]$. The performance of the two algorithms was measured in \abbr{matlab} and the results are summed up in table~\ref{tab:system_solution_performance}.
  1697. \begin{table}[b]\footnotesize
  1698. \caption{Initial answer calculation time using pseudoinversion and row reduction. \label{tab:system_solution_performance}}
  1699. \centering
  1700. \begin{tabular}{ccc}
  1701. \textsc{number of diodes} & \textsc{pseudoinversion} & \textsc{row reduction} \\
  1702. \hline
  1703. 3 & 0.234\thinspace{}ms & 0.617\thinspace{}ms \\
  1704. 4 & 0.251\thinspace{}ms & 0.836\thinspace{}ms \\
  1705. 5 & 0.249\thinspace{}ms & 0.983\thinspace{}ms \\
  1706. 6 & 0.233\thinspace{}ms & 1.166\thinspace{}ms \\
  1707. \hline
  1708. \end{tabular}
  1709. \end{table}
  1711. \noindent{}\marginpar{control accuracy}To test the accuracy of the open loop control system, a luminaire consisting of three, red green and blue, diodes was set up to deliver 50 lumen of white light ($x=1/3$ and $y=1/3$). Required tristimulus values calculated by eq.~\ref{eq:xyY2XYZ} are: $X=50$, $Y=50$ and $Z=50$. The luminaire was placed in the input port of the integrating sphere. Diodes' properties were measured using small duty cycle (\textit{d}~=~0.1\%) pulsed current so that the junction temperatures and heatsink temperature are almost equal to ambient temperature. Required duty cycles were calculated using equation~\ref{eq:openLoopColor} and previously measured calibration data.
  1713. \begin{figure}[!ht]% Open loop tristimulus
  1714.   \centering
  1715.   \includegraphics{graphics/System/OpenLoop/open_loop_tristimulus_deltaE.pdf}
  1716.  \caption{Measured tristimulus values of the open loop luminaire control system (left) and corresponding colour difference from the initial colour point (right). Self-heating of the \led{}s decreases their luminous flux therefore their tristimulus values decrease. The $X$ value experiences the biggest drop because of the high red diode temperature dependency.}
  1717.  \label{fig:OpenLoopTristimulus}
  1718. \end{figure}
  1720. \noindent{}The self-heating of the diodes yields a drop of the tristimulus values (fig.~\ref{fig:OpenLoopTristimulus}). After ca. two minutes, the colour point moved past the just noticeable difference level. The colour shift is dependent on the junction temperature therefore it is affected by the thermal properties of the heat flow path and the variations of the ambient temperature. Transient response of the luminaire is connected with the thermal capacities of the system components. Fast transient (time constant less than one min.) in the beginning of the experiment is connected with the heating of the diode structure and the aluminum \pcb, the diodes are mounted on. Slow transient (time constant ca. ten minutes) is the result of the heatsink thermal mass.
  1722. \begin{figure}[!ht]%
  1723.         \centering
  1724.         \includegraphics{graphics/System/IV2XYZ/luminaireModelOL.pdf}
  1725.         \caption{Simulation model of \rgb{} luminaire using a current-voltage diode model.}
  1726.         \label{fig:OpenLoopIV2XYZModel}
  1727. \end{figure}
  1729. \begin{figure}[!ht]%
  1730.         \centering
  1731.         \includegraphics{graphics/System/IV2XYZ/IV2XYZ_luminare_XYZ.pdf}
  1732.         \caption{Tristimulus values of the simulated \rgb{} luminaire. Diodes driven with nominal 0.7\thinspace{}\textsc{a} current.}
  1733.         \label{fig:OpenLoopIV2XYZTristimulus}
  1734. \end{figure}
  1738. % Temperature feed forward
  1739. \subsection{Temperature feed forward}\led{} characteristics are mostly affected by the junction temperature. Therefore a control scheme basing on temperature measurements is most convenient. The \led{}s junction (active, light-emitting region) temperature cannot however be measured directly. The practical approach is to measure heatsink or ambient temperature and approximate the junction temperature by means of \led{} thermal model. Garcia et~al. \cite{Garcia2008,Garcia2009} designed a flux estimator to control the output flux during steady state and transients conditions.
  1741. \begin{figure}[!ht]% Temperature feed forward
  1742.         \centering
  1743.         \includegraphics{graphics/System/tff.pdf}
  1744.         \caption{Temperature feed forward luminaire control system.}
  1745.         \label{fig:SystemTFF}
  1746. \end{figure}
  1748. \noindent{}In the \tff{} control loop controller adjusts the calibration matrix based on the measured or estimated junction temperature. The drawback of this method is that the light output (colour, intensity) is not the directly controlled variable. Every variation in \led{} model will result in long-time error~\cite{Subramanian2002}.
  1750. % Flux feedback
  1751. \subsection{Flux feedback}Flux feedback control scheme can maintain constant flux output for each of the colours. This can be done by placing photodiodes near each \led{} colour or by using one photodiode and a time sequence measurement of mixed light. This can easily overcome flux variations with ageing but cannot compensate for peak wavelength shifts.
  1753. \begin{figure}[!ht]% Flux feedback
  1754.         \centering
  1755.         \includegraphics{graphics/System/ffb.pdf}
  1756.         \caption{Flux feedback luminaire control system.}
  1757.         \label{fig:SystemFFB}
  1758. \end{figure}
  1760. \noindent{}Special care must be taken with placement of the sensor. Ambient light can distort flux measurement therefore either the placement of the sensor should minimise its influence or the ambient light can be measured independently and the subtracted from the measured \led{} flux.
  1762. % Flux feedback with temperature feed forward
  1763. \subsection{Flux feedback with temperature feed forward}A combination of temperature feed forward and flux feedback methods can produce a very accurate control scheme because it can compensate for both wavelength shift and flux change due to ageing and temperature variations. However, the method relies on model data describing relation between wavelength shift and temperature ($\ud{}λ/\ud{}T$).
  1765. \begin{figure}[!ht]% Flux feedback with temperature feed forward
  1766.         \centering
  1767.         \includegraphics{graphics/System/tffffb.pdf}
  1768.         \caption{Flux feedback with temperature feed forward luminaire control system.}
  1769.         \label{fig:SystemFFBTFF}
  1770. \end{figure}
  1772. % Colour coordinates feedback
  1773. \subsection{Colour coordinates feedback}A very high accuracy of colour reproduction can be achieved using colour coordinates feedback because the mixed light is a directly controlled variable. No model parameter variation or temperature difference can influence steady state output.
  1775. \begin{figure}[!ht]% Colour coordinates feedback
  1776.         \centering
  1777.         \includegraphics{graphics/System/ccfb.pdf}
  1778.         \caption{Colour coordinates feedback luminaire control system.}
  1779.         \label{fig:SystemCCFB}
  1780. \end{figure}
  1782. \noindent{}The challenge in designing such converter is choosing an appropriate sensor. There are many colour sensors available on the market. General trichromatic sensors consist of three photodiodes with spectral responses in red, green and blue color spectrum. Sometimes colour sensors cover also \ir{} and \uv{} parts of the spectrum.
  1784. \begin{figure}[!ht]% Colour sensor spectral responses
  1785.         \centering
  1786.         \includegraphics{graphics/System/ColorSensors/sensitivity.pdf}
  1787.         \caption{Spectral responses of colour sensors: Taos~\abbr{tcs~230}, MaZeT~\abbr{mcs3as}, MaZeT~\abbr{mtcs}i\abbr{ct} and MaZeT~\abbr{mmcs6} (top to bottom). Data taken from individual sensor datasheets.}
  1788.         \label{fig:ColourSensorsSpectralResponses}
  1789. \end{figure}
  1791. %The sensor spectral response should match the \cie{}~\abbr{1931} 2° colour matching functions. Any difference between sensors spectral response and respective colour matching function would result in error in calculating colour position on $xy$-plane.
  1793. The temperature and ageing problems are not entirely solved by using a colour sensor, because this problems can affect the sensor itself. The expected sensor lifetime should match the expected lifetime of the whole lamp.
  1795. In general case of n-chromaticity filter, the response of the filter has to be transformed from its colour space into working colour space of the luminaire. At least three measurements are necessary in order to calibrate the sensor. Typically, the calibration procedure includes more than three points to increase the accuracy of the calibration. Calibration matrix for n-chromaticity luminaire is a 3×\textit{n} matrix
  1796. \begin{equation} \label{eq:ColourSensorCalibrationMatrix}
  1797. \mathbi{C} =
  1798. \left[ \begin{array}{ccc}
  1799.    c_{11} & \ldots{} & c_{n1} \\
  1800.    c_{12} & \ldots{} & c_{n2} \\
  1801.    c_{13} & \ldots{} & c_{n3}
  1802. \end{array} \right]
  1803. \end{equation}
  1804. Taking $m$ sensor measurements and comparing them to $X, Y$ and $Z$ measured values yields two matrices
  1805. \begin{equation} \label{eq:ColourSensorMeasurements}
  1806. \mathbi{S} =
  1807. \left[ \begin{array}{ccc}
  1808.    s_{11} & \ldots{} & s_{m1} \\
  1809.    \vdots & \ddots{} & \vdots \\
  1810.    s_{1n} & \ldots{} & s_{mn}
  1811. \end{array} \right] \qquad
  1812. \mathbi{XYZ} =
  1813. \left[ \begin{array}{ccc}
  1814.    X_{1} & \ldots{} & X_{m} \\
  1815.    Y_{1} & \ldots{} & Y_{m} \\
  1816.    Z_{1} & \ldots{} & Z_{m}
  1817. \end{array} \right]
  1818. \end{equation}
  1819. Because of the finite resolution and accuracy of the measurement, the system can be approximated by equation
  1820. \begin{equation} \label{eq:ColourSensorSystem}
  1821. \mathbi{XYZ} \approx{} \mathbi{C} \cdot \mathbi{S}
  1822. \end{equation}
  1823. Calibration matrix \mathbi{C} can be approximated using least square method. First eq.~\ref{eq:ColourSensorSystem} is transposed yielding
  1824. \begin{equation} \label{eq:ColourSensorSystemTransposed}
  1825.  \mathbi{S}^{T} \cdot \mathbi{C}^{T} \approx{} \mathbi{XYZ}^{T}
  1826. \end{equation}
  1827. and the least square solution is calculated
  1828. \begin{align} \label{eq:ColourSensorCalibrationMatrixLeastSquareSolution}
  1829.  \widehat{\!\mathbi{C}}^{T}\!\!\!\!\; &= \left(\mathbi{S} \cdot \mathbi{S}^T\right)^{-1} \mathbi{S} \cdot \mathbi{XYZ}^T \nonumber \\
  1830.  \widehat{\!\mathbi{C}} &= \left(\mathbi{S} \cdot \mathbi{XYZ}^T\right)^T \cdot \left(\left(\mathbi{S} \cdot \mathbi{S}^T\right)^{-1}\right)^T \nonumber \\
  1831.  \widehat{\!\mathbi{C}} &= \mathbi{XYZ} \cdot \mathbi{S}^T \cdot \left(\mathbi{S} \cdot \mathbi{S}^T\right)^{-1}
  1832. \end{align}
  1834. \noindent{}A generic trichromatic colour sensor was calibrated with \rgb{} and \rgbw{} luminaire using random operating points (varying $x$ and $y$ colour coordinates and intensity values) and equation~\ref{eq:ColourSensorCalibrationMatrixLeastSquareSolution}. Results (fig.~\ref{fig:ColourSensorCalibration}) show that the \rgb{} system can be accurately measured using trichromatic sensor for \rgbw{} luminaire measurement yields high colour variations.
  1835. \begin{figure}[!ht]% Colour sensor calibration
  1836.         \centering
  1837.         \includegraphics{graphics/System/CCFB/color_sensor_calibration.pdf}
  1838.         \caption{Colour sensor calibration points for \rgb{} luminaire (left) and \rgbw{} luminaire (right). Circles show colour points measured by spectrometer. Dots show colour points measured by the sensor transformed using least square solution matrix~  $\widehat{\!\mathbi{C}}$.}
  1839.         \label{fig:ColourSensorCalibration}
  1840. \end{figure}
  1842. \subsection{Current-voltage model based colour control}\label{ssec:iv2xyzColourControl}
  1844. \begin{figure}[!ht]% % IV to XYZ color control
  1845.         \centering
  1846.         \includegraphics{graphics/System/iv2xyz/color_control_loop_IV2XYZ_model.pdf}
  1847.         \caption{Colour control loop utilising \ivxyz{} models. Measured current and voltage of three diode strings are converted into instantaneous tristimulus values. These values are used as a feedback for three current controllers.}
  1848.         \label{fig:color_control_IV2XYZ}
  1849. \end{figure}
  1850. A current-voltage model described in chapter~\ref{sec:iv_model} can be used in a \rgb{} luminaire colour control loop instead of optical feedback~(fig.~\ref{fig:color_control_IV2XYZ}). Voltages and currents of all colour strings present in the luminaire need to be measured. Electrical parameters are converted into tristimulus values using a model created during calibration. Respective tristimulus values are summed and subtracted from the command values. Resulting errors are fed into three \textsc{pi} controllers which control red, green and blue diode strings. The bandwidth of these controllers can be limited to few hundred hertz as human eye cannot see light changes above circa 100\thinspace{}Hz. \led{} colorimetric properties are dependent on the forward current and the junction temperature. The active area and internal \led{} structure heats up in the millisecond time range and is too fast for the human eye to perceive. \led{} junction will be also affected by the heatsink temperature. The time constant of a heatsink can be in the range of few seconds to few minutes and colour control loop should compensate for this change.
  1852. \begin{figure}[!ht] % (I)V to XYZ color control
  1853.         \centering
  1854.         \includegraphics{graphics/System/IV2XYZ/color_control_loop_(I)V2XYZ_model.pdf}
  1855.         \caption{Colour control loop where the current feedback is taken from the current controller.}
  1856.         \label{fig:color_control_IsetV2XYZ}
  1857. \end{figure}
  1858. The current information is present in the colour control loop in
  1859. the form of current command. This information can be fed back to \led{} model to decrease the number of sensors needed, as shown in fig.~\ref{fig:color_control_IsetV2XYZ}. Current command differs from actual current only during transients therefore this control scheme is best used with \am{} dimmed luminaires.
  1861. \begin{figure}[!ht] % IV2XYZ luminare model
  1862.         \centering
  1863.         \includegraphics{graphics/System/IV2XYZ/luminaireModel_IV2XYZ_PI.pdf}
  1864.         \caption{Simulation model of a trichromatic \rgb{} luminaire controlled by three \textsc{pi} controllers and using current-voltage diode model.
  1865.         \label{fig:IV2XYZ_luminare_model}}
  1866. \end{figure}
  1868. \section[Optimisations]{optimisations}
  1869. There are two possibilities for optimising \led{} based luminaires' colorimetric properties. The number of different colours present in the luminaire will affect achievable gamut, luminous flux and the light quality of the luminaire. Žukauskas et~al. analysed the properties of polychromatic \led{} luminaires~\cite{Zukauskas2001,Zukauskas2004a}. By increasing the number of different colours present in the luminaire from three to four the maximum colour rendering index can be increased from approximately 90 to approximately 98. Also the choice of the peak wavelengths of the particular diode colour will influence both the \cri{} and efficacy of the luminaire~\cite{Zukauskas2004a,Zukauskas2004b}. Chapter~\ref{ssec:linear_programming} shows an example how an optimisation procedure can influence the luminaire design.
  1871. Another possibility for optimising colorimetric properties is by the diode control. Optimisation is generally possible when a system is underdetermined systems. In polychromatic luminaires consisting of at least four different colours, like \textsc{rgba} or \textsc{rgbw} luminaires, one colour point can be reached by many different combinations of the diode intensities (figure~\ref{fig:openLoop4DiodesExampleGraphicalSolution}). However, a \rgb{} luminaire has three control variables (intensity of each diode colour) and three colour equation therefore there is no possibilities to optimise the control as there is exactly one solution to the system. A hybrid diode control described in chapter~\ref{ssec:hybrid_pwm_am} can be used to increase the number of variables in the system. Each diode is controlled with both duty cycle and peak current therefore the number of variables present in the system doubles. Chapter \ref{ssec:Optimal_control_hybrid_dimming} shows the benefits of using hybrid dimming in a trichromatic \rgb{} luminaire.
  1873. It is worth noting that even though the optimisation procedures will find many operating points that yield the same tristimulus values, the spectrum at these operating points will vary. It means that the direct light seen by the observer will not change with changing operating point but the light reflected from various objects may change due to metamerism. This creates a practical limit of how fast should the operating point change so that the observer will not notice the change created by the optimisation procedure. Also multiple lamps, each running an optimisation routine, generating light wi the same tristimulus values may produce different spectra.
  1875. \subsection{Linear programming}\label{ssec:linear_programming}
  1876. Linear programming is a technique of finding the optimum (maximum or minimum) of a objective function given a set of constrains on the function variables. Linear programs are typically solved by simplex algorithms or interior point methods. Simplex algorithms search the corners of the feasible solution area (limited by the constrains) as the solution, if exists, lays on one or more corner points. Interior point methods are more suitable to big optimisation problems and therefore will not be discussed.  
  1878. \marginpar{example}Given a luminaire consisting of four: red, green, blue and white colour diodes described by calibration matrix
  1879. \begin{equation} \label{eq:LPExampleCalibrationMatrix}
  1880. \mathbi{C} =
  1881. \left[ \begin{array}{cccc}
  1882. 60      & 6             & 17    & 74    \\
  1883. 25      & 26    & 5             & 77    \\
  1884. 0               & 4             & 100   & 92
  1885. \end{array} \right]
  1886. \end{equation}
  1887. find the maximum luminous flux at $x=0.25,\,y=0.25$ colour point. \marginpar{objective function}This problem can be solved by linear programming by maximising objective function
  1888. \begin{equation} \label{eq:LPExampleObjectiveFunction}
  1889. max: F = 25d_{1} + 26d_{2} + 5d_{3} + 77d_{4}
  1890. \end{equation}
  1891. where $d_{1}…d_{4}$ are relative intensities of the diodes with respect to their nominal intensities. The coefficients for equation \ref{eq:LPExampleObjectiveFunction} are taken from the second row of calibration matrix containing the $Y$ tristimulus values in lumens, and therefore, the luminous flux of the diodes. \marginpar{constrains}The solution should be constrained by the feasible relative fluxes of the diodes $0<d_{i}<1$ and the desired output colour point.
  1892. \begin{equation}
  1893. x = \frac{\sum_{i} d_{i}X_{i}}{\sum_{i} d_{i}(X_{i}+Y_{i}+Z_{i})} \qquad{}
  1894. y = \frac{\sum_{i} d_{i}Y_{i}}{\sum_{i} d_{i}(X_{i}+Y_{i}+Z_{i})}
  1895. \end{equation}
  1896. Computing first equation yields
  1897. \begin{align}
  1898. 0.25 &= \frac{ 60d_{1}+6d_{2}+17d_{3}+77d_{4} }{ d_{1}(60+25) +d_{2}(6+26+4)+d_{3}(17+5+100)+d_{4}(74+77+92) } \nonumber \\
  1899. 0.25 &= \frac{ 60d_{1}+6d_{2}+17d_{3}+77d_{4} }{ 85d_{1} +36d_{2}+122d_{3}+243d_{4} } \nonumber
  1900. \end{align}
  1901. \begin{align}
  1902. 85d_{1} +36d_{2}+122d_{3}+243d_{4} &= 240d_{1} +24d_{2}+68d_{3}+296d_{4} \nonumber \\
  1903. -155d_{1} +12d_{2}+54d_{3}-53d_{4} &= 0  \label{eq:LPExampleEqualityConstrain1}
  1904. \end{align}
  1905. and similarly for the second equation
  1906. \begin{equation} \label{eq:LPExampleEqualityConstrain2}
  1907. \phantom{0+}-15d_{1} -68d_{2}+102d_{3}-65d_{4} = 0 \phantom{+ 20d_{1} +26d_{2}+6d_{3}+2d_{4} } %phantoms to keep aligned with the first constraint
  1908. \end{equation}
  1909. Objective function \ref{eq:LPExampleObjectiveFunction} together with the equality constrains eq.~\ref{eq:LPExampleEqualityConstrain1}, eq.~\ref{eq:LPExampleEqualityConstrain2} and duty cycle constrains $0<d_{i}<1$ form a linear problem. The solution for this problem is 99.1\thinspace{}lm at $d_{1} = 0.0477$, $d_{2} = 0.5336 $ and $d_{3} = d_{4} = 1$.
  1913. %\marginpar{maximal flux within luminaires gamut}
  1914. Using linear programming to find the maximum output flux of the luminaire at given colour coordinates can be applied to find the limits of the device inside its gamut. Figure~\ref{fig:LPMaxFluxGamut} compares the performance of two luminaires.
  1915. \begin{figure}[!ht]% LPMaxFluxGamut
  1916.         \centering
  1917.                 \includegraphics{graphics/System/Optimizations/LinearProgramming/MaxFluxGamut.pdf}
  1918.         \caption{Comparison of maximum flux achievable by luminaires. Contour plotted every 5\thinspace{}lm. Left lumanaire described by calibration matrix \ref{eq:LPExampleCalibrationMatrix}. Right luminaire has increased number of green diodes. } \label{fig:LPMaxFluxGamut}
  1919. \end{figure}
  1920. Considering first system, white diode's colour point lays on the Planckian locus and the point corresponding to red, green and blue diode driven at maximal current lays below black-body locus. Therefore when all diodes are driven with maximum current the colour point will lay between these two points, under the Planckian locus. When generation of white light needs to be optimised, the colour point of summary of red, green and blue diodes should be as close to the locus as possible. In order to do this, the number of green diodes has been increased by a factor of 2.5× to form a second system.
  1922. The cross section of maximum flux plots across black-body locus is shown in figure~\ref{fig:LPWhiteLightSolutions}.
  1923. \begin{figure}[!ht]% LPWhiteLightSolutions
  1924.         \centering
  1925.         \includegraphics{graphics/System/Optimizations/LinearProgramming/LPExampleBlackbodyMaxFlux.pdf}
  1926.         \caption{Maximum luminous flux for white light with various colour temperatures [K]. Left graphs show the system defined by calibration matrix \ref{eq:LPExampleCalibrationMatrix}. Maximum flux for 6500\thinspace{}\textsc{k} white light is 115.2\thinspace{}lm. By increasing the number of green diodes by the factor of 2.5× (right graphs) the flux at the same colour point increases to 171.6\thinspace{}lm. Zero luminous flux indicates that the desired colour point lays outside gamut of the luminaire.} \label{fig:LPWhiteLightSolutions}
  1927. \end{figure}
  1928.         In the case of first system, the green diodes are clearly the limiting factor, as they are driven at their maximum luminance across all colour temperatures of the output light. Increased number of green diodes yields increased maximal flux in colour temperature range from around 2000\thinspace{}\textsc{k} to 10000\thinspace{}\textsc{k} while increasing the number of blue diodes would not change anything.
  1930. \subsection{Optimal control using hybrid dimming}\label{ssec:Optimal_control_hybrid_dimming}
  1931. Increasing the number of control variable yields a possibility of increased control of a trichromatic \led{} luminaire. An experiment was conducted in order to test the gains from hybrid dimming. A red, green and blue diode were chosen as this combination is most commonly present in trichromatic luminaires. Diodes' tristimulus values and electrical quantities were measured and fitted to a polynomial (figure~\ref{fig:DimmingPWMAM_XYZ}). This data represents every achievable colour point for each of the diodes. For simplicity constant heatsink temperature was used during data acquisition.
  1933. \marginpar{maximising luminous efficacy}Maximising luminous efficacy (lm/\textsc{w}) means minimising the input power for given flux level.
  1934. \begin{equation}
  1935.         \textrm{maximise: } f\!\;(d_{n},I_{peak,n}) = \sum_{n} P_{n}, \qquad{} n = \{\textrm{r, g, b}\}
  1936. \end{equation}
  1937. subject to
  1938. \begin{align}
  1939.         &X = \sum_{n}X_{n}, \qquad{}
  1940.          Y = \sum_{n}Y_{n}, \qquad{}
  1941.          Z = \sum_{n}Z_{n}  \nonumber \\
  1942.         &0 < I_{peak,min} < I_{peak,n} < I_{f,n,max} \\
  1943.         &0 < d_{min} < d_{n} < 1 \nonumber
  1944. \end{align}
  1946. \noindent{}Simulations (figure~\ref{fig:EfficacyIncreaseHybridDimming}) show an increase in luminous efficacy in hybrid case compared to both \pwm{} and \am{}. The increase depends on the value of the flux. The lower the flux value the more gain from the hybrid dimming approach.
  1947. \begin{figure}[!ht]% Efficacy increase at hybrid dimming
  1948.         \centering
  1949.         \includegraphics{graphics/System/Optimizations/HybridDimming/EfficacyIncrease.pdf}
  1950.         \caption{Luminous efficacy of hybrid, \pwm{} and \am{} dimmed luminaires. Efficacies calculated for colour point laying on a blackbody locus for 25 and 50 lumens. Zero efficacy means that this the colour point of this colour temperature cannot be generated by this luminaire.}\label{fig:EfficacyIncreaseHybridDimming}
  1951. \end{figure}
  1952. \clearpage % *** manual new page 84/85
  1953. \noindent{}\marginpar{maximising flux}Another possibility of optimisation is the maximisation of the luminous flux. An interior-point search algorithm was used to find a maximum flux value for a given $x$,~$y$ pairs. Using the same model three dimming approaches were compared. By setting the bounds on duty cycle to be equal to 100\thinspace{}\% an \am{} dimmed luminaire was simulated. Similarly when the peak current was forced to nominal current value for these diodes a \pwm{} dimmed luminaire was simulated. The optimisation procedure was run for each dimming technique. The optimisation problem was formulated as a maximisation of luminous flux
  1954. \begin{equation}
  1955.         \textrm{maximise: } f\!\;(d_{n},I_{peak,n}) = \sum_{n} Y_{n}, \qquad{} n = \{\textrm{r, g, b}\}
  1956. \end{equation}
  1957. subject to
  1958. \begin{align}
  1959.         &x = \frac{\sum_{n}X_{n}}{\sum_{n}(X_{n}+Y_{n}+Z_{n})}, \qquad{}
  1960.         y = \frac{\sum_{n}Y_{n}}{\sum_{n}(X_{n}+Y_{n}+Z_{n})} \nonumber \\
  1961.         &0 < I_{peak,min} < I_{peak,n} < I_{f,n,max} \\
  1962.         &0 < d_{min} < d_{n} < 1 \nonumber
  1963. \end{align}
  1964. where the $I_{f,n,max}$ is the maximum forward current of the $n$-th diode. The value of $I_{peak,min}$ is set to 0 or $I_{f,n,max}$ for \am{} and \pwm{} dimming, respectively. Similarly, the $d_{min}$ variable is set to 0 or 1 for \pwm{} and \am{} dimming, respectively. For hybrid dimming both variables are set to zero.
  1966. \begin{figure}[!ht]% HybridVSPWM
  1967.         \centering
  1968.         \includegraphics{graphics/System/Optimizations/HybridDimming/Fhybrid-Fpwm.pdf}
  1969.         \caption{Increase of the luminous flux in hybrid dimmed luminaire compared to \pwm{} dimmed luminaire.}\label{fig:HybridVSPWM}
  1970. \end{figure}
  1972. \begin{figure}[!ht]% HybridVSAM
  1973.         \centering
  1974.         \includegraphics{graphics/System/Optimizations/HybridDimming/Fhybrid-Fam.pdf}
  1975.         \caption{Increase of the luminous flux in hybrid dimmed luminaire compared to \am{} dimmed luminaire.}\label{fig:HybridVSAM}
  1976. \end{figure}
  1978. Results (fig.~\ref{fig:HybridVSPWM} and~\ref{fig:HybridVSAM}) are plotted as an increase of luminous flux in hybrid dimmed luminaire compared to \pwm{} and \am{} dimmed luminaire. An increased luminaire gamut the blue-cyan colours is observed in both cases.
  1980. Hybrid \pwm{}/\am{} dimming technique presented in this dissertation is only one variation of current shape that may be used to drive \led{}s. Tse~et~al. \cite{Tse2009a,Tse2009b,Tse2009c,Tse2009d} proposed using two levels of \dc{} current and alternating between them. Colorimetric benefits of using such a current shape has not yet been studied.
  1981. % chapter luminaire_control (end)
  1983. %%%%%%%%%%%%%%%
  1984. %% Chapter 5 %%
  1985. %%%%%%%%%%%%%%%
  1986. \chapter{Power converter} % (fold)
  1987. \label{cha:power_converter}
  1988. Existing \ac{} diodes that can be directly connected to the electric grid are not powerful enough for projectors and moving head applications. Power electronics circuitry is necessary to convert \ac{} grid power into the form applicable to high power \dc{} \led{}s. The structure and complexity of the power electronics converter depends on the power level of installed \led{}s and type of the colour control.
  1990. All electronic circuitry connected to the utility grid have to meet strict demand on the drawn and injected current. \textsc{en}~61000-3-2 standard classifies the device by its application and based on this appropriate limits are imposed. \textsc{en}~61000-3-3 puts a limits on observable light flicker and voltage fluctuations.
  1992. In simplest low power applications all diodes can be driven with a single-stage power converter. This ensures low cost of the electronics and simplicity of the design. The drawback of a single stage converter is that the diode control is not fully decoupled from the grid current control. Any disturbance in the grid may be observable in the output light. Therefore this kind of converters are limited to applications where precise light control is not important e.g. street lighting.
  1994. Power converters for driving high power, polychromatic luminaires are typically composed of a multi-stage converter, each specialised in its own function. In the simplest form one converter (called a power factor correction stage) is controlling the grid current and keeping the input current spectra within defined standards while a second converter drives the light-emitting diodes.
  1996. \begin{figure}[!ht]% Three stages of power conversion
  1997.         \centering
  1998.                 \includegraphics{graphics/Converter/PowerConverterLuminaire.pdf}
  1999.         \caption{Three stages of power conversion from \ac{} grid source to \led{} load. The \dc{}/\dc{} converter provides galvanic isolation between the grid and the load. Three \led{} drivers provides power to three: red, green and blue diode strings.}
  2000.         \label{fig:PowerConverterLuminaire}
  2001. \end{figure}
  2003. \noindent{}The approach chosen for this project is a three stage power converter (fig.~\ref{fig:PowerConverterLuminaire}) where the first \textsc{pfc} stage is an interleaved buck converter with output voltage of 395\thinspace{}\textsc{v}. The role of the second stage, a phase shifted full bridge converter, is to lower the high voltage from the previous stage to a safe 30\thinspace{}\textsc{v} level and to provide galvanic isolation between its input and output. The third stage, described in this chapter, drives the \led{}s.
  2005. In design process of the \led{} driver, the two first stages of the power converter are modelled as an ideal voltage source.
  2007. %% Converter requirements
  2008. \section[Converter requirements]{converter requirements}
  2010. The output voltages and currents of the converter are determined by the load. In case of \led{}s the current-voltage characteristics is well defined and does not vary much. Also, when using a \pwm{} dimming method, there are only two operating points: when the diode is turned off and when the diode is supplied with nominal current value.
  2012. \begin{figure}[!ht]% Luminaire dimming methods
  2013.         \centering
  2014.                 \includegraphics{graphics/Converter/DimmingMethods.pdf}
  2015.         \caption{Three approaches to dimming of individual diodes in the luminaire. Each colour has a separate driver which has a fixed current value and is controlled by an external \pwm{} intensity control (left). One driver supplies power to all diodes and individual diodes are shorted to create individual intensity control (middle). Each colour of the diodes is driven by a separate driver with an adjustable current (right).}
  2016.         \label{fig:DimmingMethods}
  2017. \end{figure}
  2019. \noindent{}Different dimming schemes may influence the driver design (fig.~\ref{fig:DimmingMethods}). The driver may operate with a fixed current and nearly constant voltage if one driver supplies only one colour string. This if very beneficial because the driver may be optimised only for a single operating point.
  2021. If one driver supplies power to all, serial connected diodes in the luminaire, and the diodes are controlled individually by shorting their terminals, the current is fixed to a single value but the voltage changes significantly from zero to the sum of all forward voltages.
  2023. Finally, each of the colours in the luminaire can be controlled by a separate \dc{} current source with variable current magnitude. This limits the operating points of the converter to the current-voltage curves of the connected diodes (e.g. fig.~\ref{fig:IV_n_Rs}). This approach gives the biggest control opportunities as the current in each diode can be adjusted individually allowing the use of \am{}, \pwm{} and hybrid dimming schemes.
  2025. The last scheme was chosen for implementing in a test luminaire. The highest degree of current control gives the opportunity to test the behaviour of hybrid dimming as well as current-voltage \led{} model.
  2027. % figure: current voltage curves of 4xCBT-90
  2029. The test luminaire was chosen to have three: red, green and blue colour strings, each consisting of four \textsc{cbt}-90 diodes. The typical voltage of a string, at 13.5\thinspace{}\textsc{a} nominal current, had a value of 9.6, 17.2 and 15.6\thinspace{}\textsc{v} for red, green and blue string, respectively. The forward voltage of a single colour diode may vary significantly therefore the converter should be designed with higher maximum power than predicted with typical forward voltages. The total nominal power of installed \led{}s is 572.4\thinspace{}\textsc{w} with string powers: 129.6\thinspace{}\textsc{w} for red, 210.6\thinspace{}\textsc{w} for blue to 232.2\thinspace{}\textsc{w} for green string.
  2031. %% Converter topologies
  2032. \section[Converter topologies]{converter topologies}
  2033. Different topologies of power converters may be used to power light-emitting diodes. The difference between topologies include the complexity, control methods, input and output voltages and the shape of input and output waveforms. \cite{Broeck2007,EricksonBook}
  2035. \begin{figure}[!ht]% Buck, synchronous buck, buck-boost, boost
  2036.         \centering
  2037.                 \includegraphics{graphics/Converter/BasicDCDCConverters.pdf}
  2038.         \caption{Basic \dcdc{} converters. Buck (top left), synchronous buck (top right), boost (bottom left) and buck-boost (bottom right).}
  2039.         \label{fig:BasicDCDCConverters}
  2040. \end{figure}
  2042. % Buck, buck-boost, boost, cuk, sepic, interleaved buck
  2043. \noindent{}Three basic topologies of \dcdc{} converters (fig.~\ref{fig:BasicDCDCConverters}) can be created by manipulating a switch cell composed of a transistor, inductor and recirculation diode: buck, boost and buck-boost. The application of one of these converters depends on the input to output voltage ratio of the converter. If output voltage is always lower than input voltage a buck topology can be used, if otherwise, a boost converter can be implemented. Buck-boost topology can be applied to both of the above cases.
  2045. For increased efficiency at high current levels, the recirculation diode present in a converter can be replaced with a transistor which is driven opposite to the other transistor. This reduces conduction power losses in the recirculation loop at the expense on increased complexity of the driving circuitry.
  2047. % Discuss isolated and non-isolated
  2049. Boost topologies require lower voltage on input port than on the load. If the lowest string voltage in the luminaire is 9.6\thinspace{}V, the boost-derived \led{} drivers would require more than 61\thinspace{}A of current when all diodes are driven with nominal current, not including driver losses.
  2051. One of the most important factors when designing a converter for solid-state light sources is the lifetime. Current state of the art diodes can operate for dozen thousands hours maintaining more than 70\thinspace{}\% of initial flux. Power converter driving these light sources should have similar lifetime. Electrolytic capacitors are often the cause of failures in switch-mode power converters~\cite{Chen2008}. It is beneficial therefore, to implement a topology that minimises their use.
  2053. \noindent{}High current electrolytic capacitors are used in power converters as filters and energy storage. Typical \dcdc{} power converter has at least two: input and output capacitors limiting current ripples and providing energy during transients.
  2055. %% Dual interleaved buck topology
  2056. \section[Dual interleaved buck topology]{dual interleaved buck topology}
  2058. \begin{figure}[!ht]% interleaved buck schematic
  2059.         \centering
  2060.                 \includegraphics{graphics/Converter/Schematic_interleaved_buck.pdf}
  2061.         \caption{Schematic of dual interleaved buck topology. 180° phase shift in control signal creates a ripple cancellation effect.}
  2062.         \label{fig:DualInterleavedBuckSchematic}
  2063. \end{figure}
  2065. \noindent{}An interleaved buck topology (fig.~\ref{fig:DualInterleavedBuckSchematic}) was investigated due to its ripple cancelling properties. When introduced, interleaved buck topology was used extensively in voltage regulation modules (\textsc{vrm}s) for driving \textsc{cpu}s. Benefits of using this topology were: increased dynamic performance, cancellation of output ripple current and spreading the current between multiple phases for easier thermal design~\cite{Consoli2001}. \textsc{cpu}s change the current demand often and require stable supply voltage for stable operation and this requires big values of output capacitances. \led{}s on the other hand operate in the working point set by the driver and do not require fast dynamics. It is therefore possible to reduce the output capacitance value to the point where film capacitors can be used instead of electrolytic capacitors.
  2067. \begin{figure}[!ht]% interleaved buck ripple attenuation
  2068.         \centering
  2069.                 \includegraphics{graphics/Converter/RippleAttenuation.pdf}
  2070.         \caption{Current ripple magnitude with respect to output voltage ($V_{g}=30\textrm{ V}$) superimposed on current-voltage characteristics of red, green and blue strings of four, series connected \textsc{cbt-90} diodes.}
  2071.         \label{fig:DualInterleavedBuckRippleAttenuation}
  2072. \end{figure}
  2074. \noindent{}Ripple cancellation effect is strongest at duty cycle equal to 0.5 as shown in figure~\ref{fig:DualInterleavedBuckSchematic} \cite{Wei2001}. The input voltage for the converter was fixed to 30\thinspace{}\textsc{v} so that the working point for the converter lays close to the half of the input voltage.
  2076. %% Small-signal model
  2077. \section[Small-signal model]{small-signal model}
  2078. In order to design a controller, first a small-signal model of the converter is created. Small-signal model averages the signals over the switching period and models the change of electrical quantities disregarding switching ripple \cite{EricksonBook}. This way the dynamics of signal variations can be easily modelled.
  2080. Two buck structures can be modelled independently and, using a superposition principle, connected together. When the upper switch is turned on, the inductor current is equal to the current drawn from the input source. When the top switch of the converter phase is off, the input source current is equal to zero.
  2081. \begin{align}
  2082.         i_{g}(t) &= i_{1} \textrm{ during }d_{1}\\
  2083.         i_{g}(t) &= 0 \textrm{ during }1-d_{1}
  2084. \end{align}
  2085. Applying small-ripple approximation we obtain
  2086. \begin{align}
  2087.         i_{g}(t) &\approx{} \avg{i_{1}}_{T_{s}} \textrm{ during }d_{1}\label{SSMig1d11}\\
  2088.         i_{g}(t) &= 0 \textrm{ during }1-d_{1}\label{SSMig1d12}
  2089. \end{align}
  2090. Similarly, for the other phase
  2091. \begin{align}
  2092.         i_{g}(t) &= i_{2} \textrm{ during }d_{2}\\
  2093.         i_{g}(t) &= 0 \textrm{ during }1-d_{2}
  2094. \end{align}
  2095. Small-ripple approximation yields
  2096. \begin{align}
  2097.         i_{g}(t) &\approx{} \avg{i_{2}}_{T_{s}} \textrm{ during }d_{2}\label{SSMig1d21}\\
  2098.         i_{g}(t) &= 0 \textrm{ during }1-d_{2}\label{SSMig1d22}
  2099. \end{align}
  2100. Equations \ref{SSMig1d11}, \ref{SSMig1d12}, \ref{SSMig1d21} and \ref{SSMig1d22} yield, on average over switching period, input current equal to
  2101. \begin{equation}
  2102.         \avg{i_{g}}_{T_{s}} = d_{1}\avg{i_{1}}_{T_{s}} + d_{2}\avg{i_{2}}_{T_{s}}
  2103. \end{equation}
  2104.  Similar procedure can be applied to output capacitor voltage and inductor currents yielding following equations
  2105. \begin{align}\label{eq:SSMNonlinearConverterEquations}
  2106.  C\deriv{\avg{v_{c}}_{T_{s}}}{t} &= \avg{i_{1}}_{T_{s}} + \avg{i_{2}}_{T_{s}} - \frac{\avg{v_{c}}_{T_{s}} - \avg{v_{f}}_{T_{s}}}{R_{d}} \\
  2107.  L_{1}\deriv{\avg{i_{1}}_{T_{s}}}{t} &= d_{1}(t)\avg{v_{g}}_{T_{s}} - \avg{v_{c}}_{T_{s}} \\
  2108.  L_{2}\deriv{\avg{i_{2}}_{T_{s}}}{t} &= d_{2}(t)\avg{v_{g}}_{T_{s}} - \avg{v_{c}}_{T_{s}}
  2109. \end{align}
  2110. The nonlinear equations describing the converter have to be perturbed and linearised in order to extract the small-signal \ac{} equations. All state variables and inputs to the converter are expressed as quiescent values with superimposed \ac{} variation.
  2111. \begin{align}\label{eq:SSMPerturbation}
  2112. \avg{v_{g}}_{T_{s}} &= V_{g} + \widetilde{v_{1}}(t) \nonumber \\
  2113. \avg{v_{f}}_{T_{s}} &= V_{f} + \widetilde{v_{f}}(t) \nonumber \\
  2114.     d_{1}(t)       &= D_{1} + \widetilde{d_{1}}(t) \nonumber \\
  2115.     d_{2}(t)       &= D_{2} + \widetilde{d_{2}}(t) \\
  2116. \avg{v_{c}}_{T_{s}} &= V_{c} + \widetilde{v_{c}}(t) \nonumber \\       
  2117. \avg{i_{1}}_{T_{s}} &= I_{1} + \widetilde{i_{1}}(t) \nonumber \\
  2118. \avg{i_{2}}_{T_{s}} &= I_{2} + \widetilde{i_{2}}(t) \nonumber
  2119. \end{align}
  2120. Substituting above equations into input gate equation yields
  2121. \begin{align} % *****
  2122.         I_{g} + \widetilde{i_{g}}(t) &= \Big(D_{1} + \widetilde{d_{1}}(t)\Big) \Big(I_{1} + \widetilde{i_{1}}(t)\Big) + \Big(D_{2} + \widetilde{d_{2}}(t)\Big) \Big(I_{2} + \widetilde{i_{2}}(t)\Big) \nonumber\\
  2123.         I_{g} + \widetilde{i_{g}}(t) &= \underbrace{D_{1}I_{1} + D_{2}I_{2}}_{\textrm{DC terms}} + \underbrace{D_{1}\widetilde{i_{1}}(t) + D_{2}\widetilde{i_{2}}(t) + I_{i}\widetilde{d_{1}}(t) + I_{2}\widetilde{d_{2}}(t)}_{\textrm{I-order linear terms}} + \nonumber\\ &\qquad\qquad\qquad + \underbrace{\widetilde{d_{1}}(t)\widetilde{i_{1}}(t) + \widetilde{d_{2}}(t)\widetilde{i_{2}}(t)}_{\textrm{II-order nonlinear terms}}
  2124. \end{align}
  2125. This equation contains three types of terms: \dc{}, time-invariant terms, first order \ac{} terms and second order \ac{} terms. If small signal assumptions are satisfied then the magnitude of second order terms is much lower than other terms and can be therefore neglected. Separation of \dc{} and \ac{} terms yields two set of equations
  2126. \begin{align}
  2127. I_{g} &= D_{1}I_{1} + D_{2}I_{2} \label{eq:SSMIgDC} \\
  2128. \widetilde{i_{g}}(t) &= D_{1}\widetilde{i_{1}}(t) + D_{2}\widetilde{i_{2}}(t) + I_{i}\widetilde{d_{1}}(t) + I_{2}\widetilde{d_{2}}(t) \label{eq:SSMIgAC}
  2129. \end{align}
  2130. Applying perturbations to $L_{1}$ inductor equation one obtains
  2131. \begin{align}
  2132.         L_{1}\deriv{\Big(I_{1} + \widetilde{i_{1}}(t)\Big)}{t} &= \Big( D_{1} + \widetilde{d_{1}}(t)\Big)\Big(V_{g} + \widetilde{v_{1}}(t)\Big) - \Big(V_{c} + \widetilde{v_{c}}(t)\Big) \nonumber \\
  2133.         L_{1}\bigg(\deriv{I_{1}}{t} + \deriv{\widetilde{i_{1}}(t)}{t}\bigg) &= D_{1}V_{g} + D_{1}\widetilde{v_{g}}(t) + V_{g}\,\widetilde{d_{1}}(t) + \widetilde{d_{1}}(t)\widetilde{v_{g}}(t) - V_{c} - \widetilde{v_{c}}(t) \nonumber \\
  2134.         L_{1}\deriv{I_{1}}{t} + L_{1}\deriv{\widetilde{i_{1}}(t)}{t} &= \underbrace{D_{1}V_{g} - V_{c}}_{\textrm{DC terms}} + \underbrace{D_{1}\widetilde{v_{g}}(t) + V_{g}\,\widetilde{d_{1}}(t) - \widetilde{v_{c}}(t)}_{\textrm{I-order linear terms}} + \!\!\!\!\!\!\!\!\underbrace{\widetilde{d_{1}}(t)\widetilde{v_{g}}(t)}_{\textrm{II-order nonlinear terms}}
  2135. \end{align}
  2136. Again, the second-order nonlinear terms are omitted and the terms are separated
  2137. \begin{align}
  2138.  0 &= D_{1}V_{g} - V_{c} \label{eq:SSML1DC}\\
  2139.  L_{1}\deriv{\widetilde{i_{1}}}{t} &= D_{1}\widetilde{v_{g}}(t) + V_{g}\,\widetilde{d_{1}}(t) - \widetilde{v_{c}}(t)\label{eq:SSML1AC}
  2140. \end{align}
  2141. Similar procedure can be applied to $L_{2}$ inductor equation yielding
  2142. \begin{align}
  2143.  0 &= D_{2}V_{g} - V_{c} \label{eq:SSML2DC}\\
  2144.  L_{2}\deriv{\widetilde{i_{2}}}{t} &= D_{2}\widetilde{v_{g}}(t) + V_{g}\,\widetilde{d_{2}}(t) - \widetilde{v_{c}}(t)\label{eq:SSML2AC}
  2145. \end{align}
  2146. Finally, the output capacitor equation is perturbed.
  2147. \begin{align}
  2148.         C\deriv{\Big(V_{c} + \widetilde{v_{c}}(t)\Big)}{t} &= \Big( I_{1} + \widetilde{i_{1}}(t)\Big) + \Big(I_{2} + \widetilde{i_{2}}(t)\Big) - \frac{\Big(V_{c} + \widetilde{v_{c}}(t)\Big) - \Big(V_{f} + \widetilde{v_{f}}(t)\Big)}{R_{d}} \nonumber \\
  2149.         C\deriv{V_{c}}{t} + C\deriv{\widetilde{v_{c}}(t)}{t} &= \underbrace{ I_{1} + I_{2} - \frac{V_{c}-V_{f}}{R_{d}}}_{\textrm{DC terms}} + \underbrace{\widetilde{i_{1}}(t) + \widetilde{i_{2}}(t) -  \frac{\widetilde{v_{c}}-\widetilde{v_{f}}}{R_{d}}}_{\textrm{I-order linear terms}}
  2150. \end{align}
  2151. The \dc{} and \ac{} terms are separated with the omission of nonlinear terms.
  2152. \begin{align}
  2153.  0 &= I_{1} + I_{2} - \frac{V_{c}-V_{f}}{R_{d}} \label{eq:SSMVcDC}\\
  2154.  C\deriv{\widetilde{v_{c}}}{t} &= \widetilde{i_{1}}(t) + \widetilde{i_{2}}(t)  - \frac{\widetilde{v_{c}}(t) - \widetilde{v_{f}}(t)}{R_{d}}\label{eq:SSMVcAC}
  2155. \end{align}
  2157. \noindent{}The \dc{} equations extracted in the procedure (\ref{eq:SSMIgDC}, \ref{eq:SSML1DC}, \ref{eq:SSML2DC} and \ref{eq:SSMVcDC}) collected below are used to calculate the quiescent operating point of the converter
  2158. \begin{align}
  2159.         I_{g} &= D_{1}I_{1} + D_{2}I_{2} \nonumber\\
  2160.         0 &= D_{1}V_{g} - V_{c} \\
  2161.         0 &= D_{2}V_{g} - V_{c} \nonumber\\
  2162.         0 &= I_{1} + I_{2} - \frac{V_{c}-V_{f}}{R_{d}}\nonumber
  2163. \end{align}
  2164. Upon finding the quiescent values of all the variables, the results are inserted into small-signal \ac{} equations (\ref{eq:SSMIgAC}, \ref{eq:SSML1AC}, \ref{eq:SSML2AC} and \ref{eq:SSMVcAC}) collected below
  2165. \begin{align}
  2166.         \widetilde{i_{g}}(t) &= D_{1}\widetilde{i_{1}}(t) + D_{2}\widetilde{i_{2}}(t) + I_{i}\widetilde{d_{1}}(t) + I_{2}\widetilde{d_{2}}(t) \nonumber\\
  2167.         L_{1}\deriv{\widetilde{i_{1}}}{t} &= D_{1}\widetilde{v_{g}}(t) + V_{g}\,\widetilde{d_{1}}(t) - \widetilde{v_{c}}(t) \\
  2168.         L_{2}\deriv{\widetilde{i_{2}}}{t} &= D_{2}\widetilde{v_{g}}(t) + V_{g}\,\widetilde{d_{2}}(t) - \widetilde{v_{c}}(t) \nonumber\\
  2169.         C\deriv{\widetilde{v_{c}}}{t} &= \widetilde{i_{1}}(t) + \widetilde{i_{2}}(t)  - \frac{\widetilde{v_{c}}(t) - \widetilde{v_{f}}(t)}{R_{d}}\nonumber
  2170. \end{align}
  2172. \begin{figure}[!ht]% Small-signal model separate equations
  2173.         \centering
  2174.         \includegraphics{graphics/Converter/SmallSignalModel/small_signal_model_separate.pdf}
  2175.         \caption{Circuits equivalent to the small-signal converter equations: (left)~input port, inductor loops, (right)~capacitor node.}\label{fig:SSModelSeparate}
  2176. \end{figure}
  2177. \begin{figure}[!ht]% Small-signal model ciralcuit equivalent
  2178.         \centering
  2179.         \includegraphics{graphics/Converter/SmallSignalModel/small_signal_model_combined.pdf}
  2180.         \caption{Complete small-signal \ac{} equivalent circuit model of an ideal two phase, interleaved buck converter.}\label{fig:SSModelCombined}
  2181. \end{figure}
  2182. \noindent{}The digital average current controller, chosen for this converter will control the phase current, regulating it, by varying the input duty cycle. Appropriate converter transfer functions must be analysed to design the dynamics of the controller.
  2183. \begin{equation}
  2184.         \widetilde{i_{1}}(s) = G_{i_{1}d_{1}}\widetilde{d_{1}}(s) +
  2185.         G_{i_{1}d_{2}}\widetilde{d_{2}}(s) + G_{i_{1}v_{g}}\widetilde{v_{g}}(s)
  2186. \end{equation}
  2187. First summand in the above equation corresponds to the impact of control signal on the phase current. Second term corresponds to the effect of the duty cycle variation in one phase on the inductor current in the other phase. This corresponds to the circulating current in the converter. Last term in the summation is the effect of variation in the input voltage on the inductor current and is treated as an external disturbance.
  2189. Control to inductor current transfer function is calculated assuming all other input variations are equal to zero:
  2190. \begin{equation}
  2191.         G_{i_{1}d_{1}}(s) = \frac{\widetilde{i_{1}}(s)}{\widetilde{d_{1}}(s)}\bigg|_{\widetilde{d_{2}}(s)=0\textrm{, }\widetilde{v_{g}}(s)=0}
  2192. \end{equation}
  2194. \begin{figure}[!ht]% Small-signal model Gi1d1 transfer function
  2195.         \centering
  2196.         \includegraphics{graphics/Converter/SmallSignalModel/small_signal_model_Gi1d1.pdf}
  2197.         \caption{Manipulation of equivalent circuit of the interleaved buck converter to find the $G_{i_{1}d_{1}}(s)$ control to output transfer function. } \label{fig:SSModelGi1d1}
  2198. \end{figure}
  2200. \noindent{}Because of the assumption $\widetilde{v_{g}}(s)=0$, the transformers in the converter model (fig.~\ref{fig:SSModelCombined}) are shorted and the circuit can be manipulated into the final form, shown in figure~\ref{fig:SSModelGi1d1}.
  2202. The inductor current $\widetilde{i_{1}}(s)$ can be calculated as
  2203. \begin{equation}       
  2204.         \widetilde{i_{1}}(s) = \frac{1}{sL_{1}+\Big( R_{d}\big\|sL_{2}\big\|\frac{1}{sC} \Big)}V_{g}\widetilde{d_{1}}(s)
  2205. \end{equation}
  2206. therefore the control to inductor current transfer function is equal to
  2207. \begin{align}
  2208.         G_{i_{1}d_{1}}(s) &= \frac{\widetilde{i_{1}}(s)}{\widetilde{d_{1}}(s)} = V_{g} \frac{1}{sL_{1} + \frac{sR_{d}L_{2}}{s^2 R_{d}L_{2}C + sL_{2} + R_{d}}} =
  2209. V_{g} \frac{1}{\frac{sR_{d}L_{2} + s^3 R_{d}L_{1}L_{2}C + s^2 L_{1}L_{2} + sR_{d}L_{1}}{s^2 R_{d}L_{2}C + sL_{2} + R_{d}}}=\nonumber\\
  2210. &= V_{g} \frac{s^2 R_{d}L_{2}C + sL_{2} + R_{d}}{s^3 R_{d}L_{1}L_{2}C + s^2 L_{1}L_{2} + sR_{d}(L_{1} + L_{2})}=\nonumber\\
  2211. &= V_{g} \frac{s^{2}\big(\frac{L_{2}}{L_{1}+L_{2}}\big)C + s\frac{1}{R_{d}}\big(\frac{L_{2}}{L_{1}+L_{2}}\big) + \big(\frac{1}{L_{1}+L_{2}}\big)}{s^2\big(\frac{L_{1}L_{2}}{L_{1}+L_{2}}\big)C + s\frac{1}{R_{d}}\big(\frac{L_{1}L_{2}}{L_{1}+L_{2}}\big) + s}
  2212. \end{align}
  2213. Noting the parallel connection of the inductors using following form
  2214. \begin{align}
  2215.          \bigg(\frac{L_{1}L_{2}}{L_{1}+L_{2}}\bigg) &= \big(L_{1}\big\|L_{2}\big)
  2216. \end{align}
  2217. we obtain the final transfer function
  2218. \begin{align}
  2219.         G_{i_{1}d_{1}}(s) &= V_{g} \bigg(\frac{L_{2}}{L_{1}+L_{2}}\bigg)
  2220.         \frac{s^{2}C+s\frac{1}{R_{d}}+\frac{1}{L_{2}}}{s^{3}\big(L_{1}\big\|L_{2}\big) C+s^{2}\frac{1}{R_{d}}\big(L_{1}\big\|L_{2}\big)+s}
  2221. \end{align}
  2223. \noindent{}Control to other phase inductor current transfer function is equal to the current variation in one phase to the duty cycle variation in the other phase:
  2225. \begin{equation}
  2226.         G_{i_{1}d_{2}}(s) = \frac{\widetilde{i_{1}}(s)}{\widetilde{d_{2}}(s)}\bigg|_{\widetilde{d_{1}}(s)=0\textrm{, }\widetilde{v_{g}}(s)=0}
  2227. \end{equation}
  2229. \noindent{}Manipulation of the basic model of the converter, similarly as for the previous transfer function, yields a simplified model shown in figure~\ref{fig:SSModelGi1d2}.
  2231. \begin{figure}[!ht]% Small-signal model Gi1d2 transfer function
  2232.         \centering
  2233.         \includegraphics{graphics/Converter/SmallSignalModel/small_signal_model_Gi1d2.pdf}
  2234.         \caption{Manipulation of equivalent circuit of the interleaved buck converter to find the $G_{i_{1}d_{2}}(s)$ control to output transfer function. } \label{fig:SSModelGi1d2}
  2235. \end{figure}
  2237. \noindent{}The capacitor voltage can be expressed as
  2238. \begin{equation}
  2239.         v_{c} = \frac{\Big(R_{d}\big\|sL_{1}\big\|\frac{1}{sC} \Big)}{sL_{2}+\Big( R_{d}\big\|sL_{1}\big\|\frac{1}{sC} \Big)}V_{g}\widetilde{d_{2}}(s)
  2240. \end{equation}
  2241. while the $\widetilde{i_{1}}(s)$ current is equal to
  2242. \begin{equation}       
  2243.         \widetilde{i_{1}}(s) = \frac{-v_{c}\phantom{-}}{sL_{1}}
  2244. \end{equation}
  2245. and the transfer function can be calculated according to following equations
  2246. \begin{align}
  2247.         G_{i_{1}d_{2}}(s) &= \frac{\widetilde{i_{1}}(s)}{\widetilde{d_{2}}(s)} = \frac{-V_{g}\frac{R_{d}}{s^{2}R_{d}L_{1}C+sL_{1}+R_{d}}}{sL_{2}+\frac{sR_{d}L_{1}}{s^{2}R_{d}L_{1}C+sL_{1}+R_{d}}} = \frac{-V_{g}R_{d}}{sL_{2}(s^{2}R_{d}L_{1}C+sL_{1}+R_{d})+sR_{d}L_{1}} =  \nonumber\\
  2248.         &=\frac{-V_{g}}{s^{3}L_{1}L_{2}C+s^{2}\frac{1}{R_{d}}L_{1}L_{2}+s(L_{1}+L_{2})}=\nonumber\\
  2249.   &= -V_{g}\bigg(\frac{1}{L_{1}+L_{2}}\bigg) \cdot \frac{1}{s^{3}\big(L_{1}\big\|L_{2}\big) C+s^{2}\frac{1}{R_{d}}\big(L_{1}\big\|L_{2}\big)+s}
  2250. \end{align}
  2252. \noindent{}The input voltage to inductor current will not be resolved symbolically due to much higher complexity.
  2254. All transfer function shown above behave like integrators because the parasitic resistances were omitted in the modelling. Adding these resistances will change the transfer functions so that they will have finite \dc{} response.  
  2256. \begin{figure}[!ht]% Small-signal model
  2257.         \centering
  2258.         \includegraphics{graphics/Converter/SmallSignalModel/small_signal_model_nonideal_combined.pdf}
  2259.         \caption{Complete small-signal \ac{} equivalent circuit model of an nonideal two phase, interleaved buck converter. Includes conduction losses in transistors ($R_{tu}$ and $R_{td}$—drain-source resistance of upper and lower transistor, respectively), inductors ($R_{L1}$ and $R_{L2}$) and in output capacitor ($R_{esr}$).} \label{fig:SSModelNonidealCombined}
  2260. \end{figure}
  2261. \clearpage % *** manual new page 100/101
  2262. \begin{figure}[!ht]% Small-signal model Gi1d1 transfer function
  2263.         \centering
  2264.         \includegraphics{graphics/Converter/SmallSignalModel/Gi1d1TransferFunction.pdf}
  2265.         \caption{Ideal and detailed (including parasitic resistances) $G_{i_{1}d_{1}}(s)$ transfer function of the \led{} driver.} \label{fig:SSModelGi1d1TransferFunction}
  2266. \end{figure}
  2268. \begin{figure}[!ht]% Small-signal model Gi1d2 transfer function
  2269.         \centering
  2270.         \includegraphics{graphics/Converter/SmallSignalModel/Gi1d2TransferFunction.pdf}
  2271.         \caption{Ideal and detailed (including parasitic resistances) $G_{i_{1}d_{2}}(s)$ transfer function of the \led{} driver.} \label{fig:SSModelGi1d2TransferFunction}
  2272. \end{figure}
  2274. %% Controller design
  2275. \section[Controller design]{controller design}
  2276. The obtained transfer functions $G_{i_{1}d_{1}}(s)$ and $G_{i_{1}d_{2}}(s)$ are further discretised using a Tustin transform, because this method preserves the stability of the system~\cite{ComputerControlledSystems}.
  2277. \begin{figure}[!ht]% Current control loop
  2278.         \centering
  2279.         \includegraphics{graphics/Converter/Control/currentControlLoop.pdf}
  2280.         \caption{Current control loop in $z$ domain.} \label{fig:ControlCurrentControlLoop}
  2281. \end{figure}
  2282. A model of the control system is then created in $z$ domain, as shown in figure~\ref{fig:ControlCurrentControlLoop}. \textsc{pi} controller parameters were tuned to obtain stable response. The colour control loop runs much lower than the current control loop, therefore the current controller does not need very high bandwidth.
  2284. \begin{figure}[!ht]%
  2285.         \centering
  2286.         \includegraphics{graphics/Converter/Control/ControllerOLBode.pdf}
  2287.         \caption{Open loop bode plots of the control system showing the stability of the design.} \label{fig:ControlOpenLoopBode}
  2288. \end{figure}
  2290. %% Current measurement
  2291. \section[Current measurement]{current measurement}
  2292. Current sensing is required  in the feedback loop for both current mode control and over-current protection of diodes. It is important that the measurement technique is fast, accurate, lossless and immune to switching noise. The dynamics of the current measurement affects the dynamics of the current control loop.
  2293. \begin{figure}[!ht]% Current measurement (Rs, Rdson)
  2294.         \centering
  2295.                 \includegraphics{graphics/Converter/CurrentMeasurementRsRdson.pdf}
  2296.         \caption{Current measurement by adding series measurement resistor $R_{s}$ (left) and using transistor $R_{ds(on)}$ resistance (right).}
  2297.         \label{fig:CurrentMeasurementRsRdson}
  2298. \end{figure}
  2300. % Current sensing resistor
  2301. \subsection{Current sensing resistor}
  2302. Most basic current measurement technique is adding a series resistor and measuring the voltage across it (fig.~\ref{fig:CurrentMeasurementRsRdson}a). This method is very accurate as high precision resistor, with small tolerance and low temperature coefficient, can be used. The drawback is, that in order to obtain reasonably high voltage signal (above 100\thinspace{}mV), a resistor with resistance $R_{s}>100\textrm{mV}/i$ has to be used. The higher the resistance the bigger the losses
  2303. \begin{equation} \label{eq:PowerLossInSenseResistor}
  2304.         p_{loss} = R_{s}i^{2}.
  2305. \end{equation}
  2306. For example, a 1\thinspace{}\textsc{a} current sensor will need minimum 100\thinspace{}mΩ resistor dissipating 100\thinspace{}mW. A 10\thinspace{}\textsc{a} current sensor needs only 10\thinspace{}mΩ resistor that will dissipate 1\thinspace{}\textsc{w} of power. Low value sensing resistors can be manufactured as a fixed length track on a converter \pcb{}.
  2308. \clearpage % *** manual new page 101/102
  2309. \noindent{}\marginpar{sense resistor placement}Placement of the sense resistor influences both conditioning circuit and protection against short circuit. Placing the circuit in the return path of the current produces a ground referenced measurement but also increases the load potential with respect to the ground. \led{} die is typically electrically isolated from its heatsink therefore this method can be used.
  2311. Sensor connected in series with the inductor provides constant current monitoring and load short circuit protection but it is the least efficient solution. It will not detect any short circuit through both transistors. The measurement has high common mode voltage dependent on the output voltage.
  2313. Connecting the sensor in series with high side transistor provides very good short circuit protection and is typically used in peak current controlled converters. Losses depend on both output current and duty cycle of the converter therefore this placement is best used in low duty cycle buck converters. Differential voltage measurement with high common mode voltage is necessary.
  2315. % MOSFET Rds(on)
  2316. \subsection{\textsc{mosfet} drain-source resistance}
  2317. Instead of using an added resistor, a parasitic resistance already present in the circuit can be used (fig.~\ref{fig:CurrentMeasurementRsRdson}b)~\cite{Forghani-zadeh2005}. The drain-source resistance of a open \mosfet{} can be used. This technique is considered as lossless as it does not introduce any additional losses apart from conduction losses of the transistor.
  2318. When transistor is on, the value of the drain current is given by
  2319. \begin{equation} \label{eq:DrainCurent}
  2320.         i_{D} = \frac{v_{sense}}{R_{DS(on)}}.
  2321. \end{equation}
  2322. Again, the minimum $R_{DS(on)}$ resistance depends on the measured current level. By using this technique one may be forced to use worse transistor than available in the current state of the art market. Another drawback is the $R_{DS(on)}$ tolerance reaching 30–40\%~\cite{Hua2006}. Datasheets typically provide the typical and maximum value of the resistance. An additional calibration circuit can be introduced to overcome this problem.
  2323. The on resistance depends also strongly on the temperature and is sensitive to low $V_{GS}$ voltage~\cite{Forghani-zadeh2005}.
  2325. \begin{figure}[!ht]% Current measurement (inductor DCR)
  2326.         \centering
  2327.         \includegraphics{graphics/Converter/CurrentMeasurementInductorDCR.pdf}
  2328.         \caption{Current measurement by inductor \textsc{dcr} current sensing (left) and improved inductor \textsc{dcr} current sensing (right).}
  2329.         \label{fig:CurrentMeasurementInductorDCR}
  2330. \end{figure}
  2332. % inductor DCR
  2333. \subsection{Inductor \textsc{dcr}}
  2334. Instead of using \mosfet{}'s on resistance a scheme based on measuring inductor \dc{} resistance can be implemented~\cite{Hua2006}. Additional \textsc{rc} network has to be added to estimate the current in the inductor. The current is given by the Ohm law in s-domain
  2335. \begin{equation}\label{eq:DCRInductorCurrent}
  2336.         I_{L} = \frac{V_{L}}{R_{L}+sL}
  2337. \end{equation}
  2338. where $V_{L}$ is the voltage on the terminals of the inductor and $R_{L}$ is the inductors \dc{} resistance. The auxiliary $R_{s}, C_{s}$ network forms a voltage divider.
  2339. \begin{align}
  2340. V_{C} &=  \frac{1/sC_{s}}{R_{s}+1/sC_{s}}V_{L} \\
  2341. V_{C} &=  \frac{1}{1+sR_{s}C_{s}}V_{L} \\
  2342. V_{L} &= (1+sR_{s}C_{s})V_{C}
  2343. \end{align}
  2344. Substituting the result to equation \ref{eq:DCRInductorCurrent} yields
  2345. \begin{align}
  2346. I_{L} &= \frac{(1+sR_{s}C_{s})}{R_{L}+sL}V_{C} \\
  2347. I_{L} &= \frac{1}{R_{L}}\left( \frac{1+sR_{s}C_{s}}{1+s\frac{L}{R_{L}}}\right) V_{C} \\
  2348. I_{L} &= \frac{1}{R_{L}}\left( \frac{1+sτ_{RC}}{1+sτ_{RL}}\right) V_{C}
  2349. \end{align}
  2350. When both time constants are equal
  2351. \begin{equation}\label{eq:DCREqualTimeConstants}
  2352.         τ_{RC} = τ_{RL} \;\Rightarrow\; R_{s}C_{s} = \frac{L}{R_{L}}
  2353. \end{equation}
  2354. voltage on the capacitor $C_{s}$ is proportional to inductor current
  2355. \begin{equation}\label{eq:DCRFinal}
  2356.         I_{L} = \frac{V_{C}}{R_{L}}\,.
  2357. \end{equation}
  2358. However, the above equation is only true when the elements are well matched (eq.~\ref{eq:DCREqualTimeConstants}). Because of the high tolerance of inductive and capacitive components and dependence of $L$ on the \dc{} bias this condition is hard to satisfy.
  2360. % Inductor dcr with improved SNR
  2361. \subsection{Inductor \textsc{dcr} with improved signal-to-noise ratio}
  2362. Signal voltage in the inductor \textsc{dcr} method depends on the series resistance of an inductor. If it is too small, the output signal will be susceptible to the interruption of noise. To increase signal to noise ratio (\textsc{snr}) a change in the measuring circuit (fig.~\ref{fig:CurrentMeasurementInductorDCR}b) has been proposed~\cite{Lethellier2002}. Additional transistors drive the measuring \textsc{rc} network. Measured signal includes the resistance of the inductor and $R_{DS(on)}$ of the transistor. This increases the signal by a factor of $(R_{DS(on)}+R_{L})/R_{L}$.
  2364. % Inductor dcr 2
  2365. \begin{figure}[!ht]% Current measurement (inductor DCR 2)
  2366.         \centering
  2367.         \includegraphics{graphics/Converter/CurrentMeasurementInductorDCR2.pdf}
  2368.         \caption{Schematic of current sensing circuit using a differential amplifier. For proper operation $R_{1}=R_{1}'$, $R_{2}=R_{2}'$ and $C_{1}=C_{1}'$ must be used.
  2369.         \label{fig:CurrentMeasurementInductorDCR2}}
  2370. \end{figure}
  2372. \subsection{Improved inductor \textsc{dcr}}\label{ssec:ImprovedInductorDCR}
  2373. Two previous inductor \textsc{dcr} measurement methods require high value of inductor resistance so that the measured signal has sufficient magnitude. On the other hand efficiency requirement calls for the lowest possible resistance. It is however possible to create a measurement circuit that does not have this drawback using a single differential amplifier with impedances matched to inductor's time constant~\cite{Dallago2000}. The schematic of this current measurement scheme is shown in figure~\ref{fig:CurrentMeasurementInductorDCR2}.
  2375. \noindent{}Impedance of paralleled $R_{1}$ and $C_{1}$ is equal to
  2376. \begin{equation} \label{eq:inductorDCR2Zrc}
  2377.         Z_{RC} = \frac{R_{1}}{1+sR_{1}C_{1}}.
  2378. \end{equation}
  2379. Currents flowing in and out $V_{a}$ and $V_{b}$ nodes are calculated as follows
  2380. \begin{equation} \label{eq:inductorDCR2I1I2}
  2381.         I_{1}=\frac{V_{1}-V_{a}}{R_{2}}\qquad{}
  2382.         I_{1}\:\!\!'=\frac{V_{a}-V_{out}}{Z_{RC}}\qquad{}
  2383.         I_{2}=\frac{V_{2}-V_{b}}{R_{2}}\qquad{}
  2384.         I_{2}\:\!\!'=\frac{V_{b}}{Z_{RC}}
  2385. \end{equation}
  2386. Assuming that amplifier input currents are negligible one can write $I_{1}=I_{1}\:\!\!'$ and $I_{2}=I_{2}\:\!\!'$. Equations~\ref{eq:inductorDCR2I1I2} can be rewritten as
  2387. \begin{align}\label{eq:inductorDCR2tf}
  2388. Z_{RC}\left(V_{1}-V_{a}\right) &= R_{2}\left(V_{a}-V_{out}\right)  \nonumber \\
  2389. Z_{RC}\left(V_{2}-V_{b}\right) &= R_{2}V_{b} \nonumber \\
  2390. Z_{RC}\left(V_{1}-V_{a}-V_{2}+V_{b}\right) &= R_{2}\left(V_{a}-V_{b}-V_{out}\right)
  2391. \end{align}
  2392. Assuming that the gain of amplifier \textit{A}~~∞ voltages $V_{a}$ and $V_{b}$ can be treated as equal $V_{a} = V_{b}$ and transfer function of the differential amplifier circuit can be calculated as
  2393. \begin{align}\label{eq:inductorDCR2tf2}
  2394. Z_{RC}\left(V_{1}-V_{2}\right) &= -R_{2}V_{out}\nonumber \\
  2395. Z_{RC}V_{in} &= R_{2}V_{out} \nonumber \\
  2396. \frac{V_{out}}{V_{in}} &= \frac{Z}{R_{2}}
  2397. \end{align}
  2398. Voltage on the inductor is equal to
  2399. \begin{align}\label{eq:inductorDCR2IL}
  2400. V_{L} &= I_{L}\left(R_{L}+sL\right) \nonumber \\       
  2401. V_{L} &= I_{L}R_{L}\left(1+sL/R_{L}\right)
  2402. \end{align}
  2403. Current sensing circuit is connected to the terminals of the inductor, therefore substituting eq.~\ref{eq:inductorDCR2Zrc} and \ref{eq:inductorDCR2IL} into equation~\ref{eq:inductorDCR2tf2} yields
  2404. \begin{align}\label{eq:inductorDCR2tfVoutIL1}
  2405. \frac{V_{out}}{I_{L}} &= \frac{R_{L}R_{1}}{R_{2}}\cdot\frac{1+sL/R_{L}}{1+sR_{1}C_{1}}
  2406. \end{align}
  2407. When time constants $R_{1}C_{1}$ and $L/R_{L}$ are equal, eq.~\ref{eq:inductorDCR2tfVoutIL1} is simplified to
  2408. \begin{equation}\label{eq:inductorDCR2tfVoutIL2}
  2409. \frac{V_{out}}{I_{L}} = \frac{R_{L}R_{1}}{R_{2}}
  2410. \end{equation}
  2411. Current signal is converted into voltage signal with gain set by the ration of $R_{1}$ and $R_{2}$ resistors and the series resistance of the inductor. This circuit requires high voltage, high bandwidth amplifier in order to accurately follow inductor current.
  2413. \begin{figure}[!ht]% Current measurement (Observer, Average)
  2414.         \centering
  2415.         \includegraphics{graphics/Converter/CurrentMeasurementObserverAverage.pdf}
  2416.         \caption{Current measurement by observer technique (left) and average current sensing (right).}
  2417.         \label{fig:CurrentMeasurementObserverAverage}
  2418. \end{figure}
  2420. % observer technique
  2421. \subsection{Observer technique}
  2422. This technique uses the inductor voltage to calculate the inductor current~\cite{Midya1997}. As the voltage-current relation of an inductor is given by
  2423. \begin{equation} \label{eq:VIinductor}
  2424.         v_{L} = L\frac{\ud{}i_{L}}{\ud{}t},
  2425. \end{equation}
  2426. the inductor current can be calculated by integrating the inductor voltage
  2427. \begin{equation} \label{eq:IVinductor}
  2428.         i_{L} = \frac{1}{L}\int{}v_{L}\ud{}t.
  2429. \end{equation}
  2430. The inductor voltage is typically much larger signal than the output of a current sensing resistor and depends only on the input voltage.
  2432. % average current sensing
  2433. \subsection{Average current sensing}
  2434. Figure~\ref{fig:CurrentMeasurementObserverAverage}b shows a simple method of measuring the average current in the inductor~\cite{Xunwei1999}. An additional \textsc{rc} network is added parallel to the recirculation transistor. Under steady state conditions, the average voltage on resistor $R_{s}$ is zero therefore average voltage on the inductor can be written as $V_{out}-\avg{V_{C}}$. The average current in the inductor is derived as
  2435. \begin{equation}
  2436.         \avg{I_{L}} = \frac{\avg{V_{L}}}{R_{L}} = \frac{V_{out}-\avg{V_{C}}}{R_{L}}\,.
  2437. \end{equation}
  2438. For accurate current sensing, the \dc{} resistance of the inductor has to be known. Values of the additional \textsc{rc} filter will affect the measurement bandwidth~\cite{Patel2007} therefore this technique is mostly applicable to control current sharing between phases in multiphase converters.
  2440. % overview
  2441. \subsection{Overview}
  2442. \noindent{}An overview of current measurement schemes is shown in table~\ref{tab:current_measurement_comparison}. An improved inductor \textsc{dcr}, described in chapter~\ref{ssec:ImprovedInductorDCR}, was chosen for the converter. It uses a single operational amplifier and a tuned \textsc{rc} circuit therefore it is not expensive. Moreover, manufacturers provide inductors with specific tolerances for the \dc{} resistance. The accuracy of the method can therefore be calculated, based on the resistance tolerance.
  2444. \input{tables/current_measurement_comparison.tex}
  2446. %% Hardware implementation
  2447. \section[Hardware implementation]{hardware implementation}
  2448. \begin{align}
  2449.         V_{in} &= 30\ \textrm{V}\nonumber \\
  2450.         V_{out} &= \left\{9.6, 17.2, 15.6\right\}\textrm{V\quad{}for red, green and blue diode  string, respectively} \nonumber \\
  2451.         I_{f} &= 13.5\ \textrm{A}\nonumber \\
  2452.         I_{phase} &= 0.5 \cdot I_{f} = 6.75\ \textrm{A}\nonumber \\
  2453.         I_{pp} &=  25\ \% \cdot{} I_{phase} = 1.6875\ \textrm{A}\nonumber \\
  2454.         f_{sw} &= 200\ \textrm{kHz}\nonumber
  2455. \end{align}
  2457. \subsection{Inductor}
  2458. In order to calculate the necessary inductance value, the inductor voltage equation can be used.
  2459. \begin{equation}\label{eq:ConverterDesignInductorVoltageEquation}
  2460.         v_{L} = L \frac{\ud{}i_{L}}{\ud{}t}
  2461. \end{equation}
  2462. During on phase of the switching period ($d\cdot{}t_{sw}$) inductor current increases by the value of peak-to-peak ripple current $I_{pp}$ and the voltage on the inductor is equal to $v_{in}-v_{out}$. Rearranging equation \ref{eq:ConverterDesignInductorVoltageEquation} and substituting appropriate variables yields
  2463. \begin{equation}\label{eq:ConverterDesignInductanceEquation}
  2464.         L = \left(V_{in}-V_{out}\right) \cdot{} \frac{V_{out}}{V_{in}\cdot{f_{sw}}} \cdot{} \frac{1}{I_{pp}}
  2465. \end{equation}
  2466. Calculating above equation yields 19.3\thinspace{}µH for red, 21.7\thinspace{}µH for green and 22.2\thinspace{}µH for blue diode string. A~closest 22\thinspace{}µH standard value for inductors was chosen to be used in the converter.
  2468. Basing on the \dc{} bias current and required inductance value a 77350-\textsc{a}7 Magnetics core has been chosen for the inductor. Magnetics design calculator was used to design and estimate the power losses of the inductor. 21 turns of 16 \textsc{awg} wire yields 21.71\thinspace{}µH at full \dc{} current bias. Core losses are estimated at 90\thinspace{}mW. \dc{} resistance of approximately 9\thinspace{}mΩ generates 410\thinspace{}mW copper losses. The total dissipated power is equal to 500\thinspace{}mW per inductor.
  2470. \subsection{Output capacitor}
  2471. The value of the output capacitor typically is chosen based on the maximum output voltage ripple magnitude and the voltage overshoot during load transient. As the \led{} is driven by the current, the voltage transients are not important. The output capacitor limits the 400\thinspace{}kHz (twice the switching frequency) ripple current flowing through the load. Therefore the only reason for the output capacitor is the \textsc{emi} issues creased by the wires connecting the driver to the diodes.
  2473. \subsection{Transistors}
  2474. \marginpar{mosfet}A \mosfet{} (metal-oxide-semiconductor field-effect transistor) transistors are ideal switches for low voltage applications. When turned on, a channel of n-type or p-type semiconductor is formed between the drain and source terminals. Because of the small on resistance (typically few mΩ) the conduction losses are much lower than in \textsc{bjt} (bipolar junction transistor) experiencing a high voltage drop between output terminals.
  2476. \begin{figure}[!ht]% MOSFET with parasitic components
  2477.         \centering
  2478.         \includegraphics{graphics/MOSFET/MOSFETwParasitic.pdf}
  2479.         \caption{Power \mosfet{} with parasitic components \cite{BaloghGateDrivers}. Nonlinearities of input capacitance $C_{iss}$, reverse transfer capacitance $C_{rss}$ and output capacitance $C_{oss}$ of \abbr{ipb80n04s3-03} transistor shown in the function of drain-source voltage.} \label{fig:MOSFETwithParasiticComponents}
  2480. \end{figure}
  2482. \noindent{}\mosfet{} with basic parasitic components is shown in the figure~\ref{fig:MOSFETwithParasiticComponents}. Drain and source inductance values depend on transistor package (typically few nH). Internal gate resistance has to be included in the driving loss calculations as its value is typically between 0.5–5\thinspace{}Ω.
  2484. \begin{figure}[!ht]% Switching Losses
  2485.         \centering
  2486.         \includegraphics{graphics/MOSFET/SwitchingLosses.pdf}
  2487.         \caption{Switching waveforms of a power \mosfet{} during turn on (left) and turn off (right). $V_{GS}$ gate-source voltage, $I_{D}$ drain current, $V_{DS}$ drain-source voltage and $P$ power losses. $V_{GS}$ waveform can be obtained from specific transistor datasheet~($V_{GS}$ vs. $Q_{G}$ figure).} \label{fig:MOSFETSwitchingLosses}
  2488. \end{figure}
  2490. % Losses
  2491. \marginpar{losses}There are three power loss mechanisms in the \mosfet{}. One is conduction losses, when transistor is fully on (fig.~\ref{fig:MOSFETSwitchingLosses}, period $t_{4}$$t_{7}$), and power is dissipated in the ohmic channel according to $R_{DS(on)}\cdot{}i_{D}^{2}$ equation. Blocking losses are typically neglected because of very low drain leakage current (for example \abbr{ipb80n04s3-03} has a 100\thinspace{}µA leakage current at $T_{j}=125\mathrm{^\circ{}C}$ which gives 4\thinspace{}mW at $v_{DS}=40 \mathrm{ V}$). The last mechanism is the switching losses. Due to high complexity, the analysis is broken into specific periods~(fig.~\ref{fig:MOSFETSwitchingLosses}).
  2493. % t1-t2
  2494. \marginpar{period $t_{1}$$t_{2}$}Before $t_{2}$, the current in the drain is zero and some voltage is present between drain and source. Only negligible blocking losses are present. \mosfet{} driver is charging the input capacitance $C_{iss}$. Transistor is in the off state. Gate-source voltage raises linearly (assuming constant current flowing into the gate). This period is defined as a turn-on delay $t_{d(on)}$.
  2496. % t2-t3
  2497. \marginpar{period $t_{2}$$t_{3}$}Gate-source voltage has reached the threshold voltage $V_{th}$ and transistor is starting to conduct current. $i_{D}$ raises but blocked voltage still remains on the output terminals. Power is dissipated due to non-zero product of $i_{D}$ and $v_{DS}$. Gate-source voltage continues to raise until it reaches plateau voltage $V_{sp}$.
  2499. % t3-t4
  2500. \marginpar{period $t_{3}$$t_{4}$}$C_{GS}$ is charged and gate current starts to charge $C_{GD}$. Due to a Miller effect, $v_{GS}$ is clamped to plateau voltage $V_{sp}$. Drain-source voltage is decreasing. The rate of $\ud{}v_{DS}/\ud{}t$ is dictated by the gate current charging the gate-drain capacitance $C_{rss}$. Assuming constant gate current, the drain-source voltage slope is not constant due to highly non-linear $C_{rss}$ (fig.~\ref{fig:MOSFETwithParasiticComponents}). Gate charge $Q_{GS} + Q_{GD}$ is the minimum charge needed to turn the transistor on. Gate charge value is used to calculate required gate drive current.
  2502. \begin{table}[b]
  2503. \footnotesize
  2504. \caption{Losses in the top \mosfet{} at 6.75\thinspace{}\textsc{a} phase current (half of the nominal 13.5\thinspace{}\textsc{a} diode current) and 16\thinspace{}\textsc{v} output voltage. \label{tab:MosfetLossesComparison}}
  2505. \centering
  2506. \begin{tabular}{lcccc}
  2507.         & \textsc{turn on} & \textsc{turn off} & \textsc{conduction} & \textsc{total} \\
  2508.         \textsc{symbol} & \textsc{losses} [\textsc{w}] & \textsc{losses} [\textsc{w}] & \textsc{losses} [\textsc{w}] & \textsc{losses} [\textsc{w}] \\
  2509.         \hline
  2510.         \textsc{irfp4321}  & 0.552 & 0.071 & 0.409 & 1.032\\
  2511.         \textsc{sud23n06}  & 0.103 & 0.146 & 1.319 & 1.568\\
  2512.         \textsc{irls3036}  & 0.616 & 1.164 & 0.067 & 1.848\\
  2513.         \textsc{irlr3705z} & 0.199 & 0.360 & 0.293 & 0.853\\
  2514.         \textsc{irlr3636}  & 0.245 & 0.400 & 0.161 & 0.806\\
  2515.         \hline
  2516. \end{tabular}
  2517. \end{table}
  2519. \marginpar{mosfet comparison}A loss model, based on above description, has been used to estimate the losses in the transistor. The model uses datasheet parameters ($R_{DS(on)}$, $C_{oss}$, $C_{rss}$, $C_{gs}$, $R_{g}$, $Q_{g}$, $Q_{gd}$, $V_{Gth}$ and $gfs$) to calculate turn on, turn off and conduction losses in the transistors. Table~\ref{tab:MosfetLossesComparison} shows the loss distribution and the total loss in the device for the top transistor. A \textsc{irlr3636} \mosfet{} was chosen based on this loss estimation. Losses in the bottom transistor are dominated by the conduction losses, as the transistor is turning on and off with a diode forward voltage drop across it. The same \mosfet{} device was chosen due to its low cost and a low $R_{DS(on)}$ value.
  2521. \subsection{Gate circuit}
  2522. A \textsc{lm5106} 100\thinspace{}V half bridge gate driver with programmable dead-time has been chosen for driving the \mosfet{}s. The driver needs only one \pwm{} signal to drive two: high and low \mosfet{}s so that one driver containing two phases requires only two \pwm{} lines.
  2524. \subsection{Current measurement}% Current measurement
  2525. The improved inductor \textsc{dcr} current measurement method, described in detail in chapter \ref{ssec:ImprovedInductorDCR}, was chosen for the converter. In the luminaire, each colour diode string is driven with a driver consisting of two buck converters. Each of these converters requires a separate current measurement circuit, therefore the current measurement circuit should be cheap.
  2527. \begin{figure}[!ht]%
  2528.         \centering
  2529.         \includegraphics{graphics/Converter/LedDriverCurrentMeasurement.pdf}
  2530.         \caption{Current measurement scheme implemented in the driver. Inductor \textsc{dcr} measurement circuit (left) and improved circuit with biased input and output (right) allows the use of cheap operational amplifiers. Potentials $V_{1}$ and $V_{2}$ measured across the inductor. } \label{fig:LedDriverCurrentMeasurement}
  2531. \end{figure}
  2533. \noindent{}A low cost \textsc{lm837} operational amplifier was used due to its price, high operating voltage ±18\thinspace{}\textsc{v} and high unity-gain bandwidth of 25\thinspace{}MHz. This amplifier does not have rail-to-rail inputs or outputs, therefore the voltage levels were shifted to an appropriate level using bias voltage (for the output) and additional bias resistors (to bias both inputs) as shown in figure~\ref{fig:LedDriverCurrentMeasurement}.
  2535. \begin{figure}[!ht]%
  2536.         \centering
  2537.         \includegraphics{graphics/Converter/CurrentMeasurement/CurrentMeasurementDynamics.pdf}
  2538.         \caption{Measured effect of difference in $R_{1}C_{1}$ and $L/R_{L}$ time constants. Measured current $I_{Lx,meas}$ is controlled by the controller. Shape of the actual current depends on the relation of time constants. $R_{1}C_{1}$ < $L/R_{L}$ (left), $R_{1}C_{1}$ equal to $L/R_{L}$ (middle) and $R_{1}C_{1}$ > $L/R_{L}$ (right)} \label{fig:LedDriverCurrentMeasurementDynamics}
  2539. \end{figure}
  2541. Time constant of the $R_{1}C_{1}$ must be equal to the time constant of the inductor for the current measurement circuit to track the current accurately, as described in chapter \ref{ssec:ImprovedInductorDCR}. The effects of time constant mismatch can be visible in figure~\ref{fig:LedDriverCurrentMeasurementDynamics}. Although, the feedback current signal experiences normal step response, the actual current response may be very different.
  2543. \begin{figure}[!ht] % Current measurement linearity
  2544.         \centering
  2545.         \includegraphics{graphics/Converter/CurrentMeasurement/CurrentMeasurementLinearity.pdf}
  2546.         \caption{Measured linearity of current sensing circuits for each phase of the converter. \textsc{adc}x values refer to the value of analog to digital conversion in the \textsc{dsp}.} \label{fig:LedDriverCurrentMeasurementLinearity}
  2547. \end{figure}
  2549. \noindent{}The linearity of the current sensing circuit has been measured in the laboratory (fig.~\ref{fig:LedDriverCurrentMeasurementLinearity}) showing a linear relationship between actual and measured current in the inductors. For the amplifier gain and bias voltage used in the experiments the achieved current measurement resolution is around 5\thinspace{}mA.
  2551. \section[Experimental results]{experimental results}
  2553. \begin{figure}[!ht] % Converter
  2554.         \centering
  2555.         \includegraphics[scale=0.6]{graphics/Converter/Converter.jpg}
  2556.         \caption{Dual interleaved buck converter built in the laboratory. Top side, show in the left picture, contains all power components: \mosfet{}s, inductors, gate drivers and filter capacitors. Bottom side of the \pcb{} contains the current measurement circuit and a connector for the control board.} \label{fig:LedDriverPicture}
  2557. \end{figure}
  2559. The design has been built in the laboratory (fig.~\ref{fig:LedDriverPicture}). The converter is controlled using a \textsc{tms320f28027} Piccolo microcontroller board. \dsp{} contains four independent \pwm{} modules, each capable of generating two \pwm{} signals with the same carrier waveform. This means that one \dsp{} can control up to four \led{} drivers. The converter has been tested using a high power \led{} load shown in figure~\ref{fig:LedDriverLEDLoad}.
  2561. \begin{figure}[!ht] % LED load
  2562.         \centering
  2563.         \includegraphics[scale=0.75]{graphics/Converter/LedLoad.jpg}
  2564.         \caption{LED load used in the experiments composed of four series connected \textsc{cbt-90} diodes.} \label{fig:LedDriverLEDLoad}
  2565. \end{figure}
  2567. \begin{figure}[!ht] % Current step
  2568.         \centering
  2569.         \includegraphics{graphics/Converter/Control/CurrentStepSimulationMeasurement.pdf}
  2570.         \caption{Measured (left) and simulated (right) current step in one of the phases of the converter. Gate driver losses are included in the measurement. } \label{fig:LedDriverCurrentStepSimulationMeasurementComparison}
  2571. \end{figure}
  2573. The dynamic response was tested by performing a step in current command. Results, shown in figure~\ref{fig:LedDriverCurrentStepSimulationMeasurementComparison} show small difference between the model and the real system. Real system shows more damped response with little overshot. This phenomenon may be explained by lack of detailed modelling of all parasitic effects e.g. inductance change with the \dc{} current value or switching losses in the transistors.
  2575. \begin{figure}[!ht] % Efficiency
  2576.         \centering
  2577.         \includegraphics{graphics/Converter/LedDriverEfficiency.pdf}
  2578.         \caption{Measured efficiency of the \led{} driver supplying four, series connected green \textsc{cbt-90} \led{}s.} \label{fig:LedDriverEfficiency}
  2579. \end{figure}
  2581. \noindent{}Efficiency was measured with the previously shown load. Power needed to supply the \textsc{dsp} are not included in the efficiency calculations as one processor can be used for controlling multiple \led{} drivers and, at the same time, act as a colour controller.
  2583. \begin{figure}[!ht] % Thermal image
  2584.         \centering
  2585.         \includegraphics[scale=0.6]{graphics/Converter/ThermalImage.png}
  2586.         \caption{Thermal image of power converter driving four, series connected \textsc{cbt-90} \led{}s with 2\thinspace{}\textsc{a}, 20\thinspace{}\textsc{w} (left) and 13.5\thinspace{}\textsc{a}, 190\thinspace{}\textsc{w} (right). No forced cooling was used.} \label{fig:LedDriverThermalImage}
  2587. \end{figure}
  2589. % chapter power_converter (end)
  2591. %%%%%%%%%%%%%%%
  2592. %% Chapter 6 %%
  2593. %%%%%%%%%%%%%%%
  2594. \chapter{Conclusions} % (fold)
  2595. \label{cha:conclusions}
  2596. The aim of this work was to research \led{} driving solutions and to design a power converter for driving high power \led{} light sources. Thorough investigation of the previous scientific work showed three distinct research areas: \led{} properties (colorimetric, thermal, dimming, lifetime, etc.), luminaire control and \led{} drivers. Luminaire control work used the results of studies on diodes' properties together with the knowledge on colour theory to create colour control engines. However, little work bridges the gap between the colour control and \led{} drivers. While designing a \led{} driver, most of the focus is on topologies and converter lifetime and almost no focus is given to the driven diode. Similarly, researchers in the field of luminaire control treat \led{} drivers as black boxes where the only design choice is the dimming method. Therefore, in order to create an intelligent \led{} power converter, all three major research areas have been examined.
  2598. \marginpar{hybrid pwm/am modulation}Previous research on pulse-width and amplitude modulation control show different effects, these methods have on diode performance and characteristics. It is therefore possible to drive the diode with any current shape to obtain different efficacy and colour points than with the two classical driving methods. To the hypothesis a diode was driven with a hybrid \pwm{}/\am{} dimming mechanism. The choice of the driving current was made based on the observation that most commercial \led{} drivers offer both methods of control, so the implementation of this method in existing or future luminaires would not require much change in hardware design. Tests conducted in laboratory conditions proved that by using the hybrid dimming mechanism, many new control opportunities became possible.
  2600. An observation that \pwm{} and \am{} methods yield an opposite colour shift while dimming white phosphor-converted \led{} led to the discovery that by using the hybrid dimming mechanism one is able to control the position of the colour point and intensity of the diode. This phenomenon was used to stabilise the colour point that shifted with the heatsink temperature changes.
  2602. Similar effect on colour point was observed with green and blue InGaN diodes. The hybrid modulation was used to minimise the peak wavelength shift yielding a colour point moving inside the MacAdam ellipse throughout its dimming range. Future investigation should verify if this behaviour can be used to increase the accuracy of the colour sensors, as the spectrum shifts are the main cause of the measurement error.
  2604. Future research on this topic should include different driving currents, as the hybrid \pwm{}/\am{} dimming is only a single example of the possible current shapes. Existing current control methods, like peak current control, are capable of controlling the current with very high bandwidth, therefore very complex current shapes can easily be used to drive light-emitting diodes.
  2606. \marginpar{luminaire control}Review of luminaire control show various feedback mechanisms that stabilise the output light colour. Methods include measuring the colour or the intensity of the light and different indirect junction temperature measurement schemes. Junction temperature information is used, together with a model of the parameters change, to estimate diode colorimetric properties. Review of colour spaces and corresponding colour distance metrics shows that the $ΔE_{ab}^{*}$ colour distance is a good measure of colour control loop accuracy.
  2608. \marginpar{luminaire control optimisation}Polychromatic luminaires, consisting of four or more \led{}s, have been previously shown to have the possibility to optimise the control of primary diodes. The optimisation procedure can maximise various lamp parameters such as: efficacy, luminous flux or colour quality. Trichromatic luminaires did not have this ability using \pwm{} or \am{} dimming methods. Research presented in this dissertation proved that by using the hybrid dimming mechanism, an increase in control allows the luminaire to optimise the same lamp parameters as in polychromatic luminaire consisting of more than three basic diodes. An increase of efficacy has been found, especially at lower intensity levels. An increased device gamut has been also shown, particularly in the cyan area of trichromatic \rgb{} luminaire.
  2610. \marginpar{current-voltage model}In this dissertation, a current-voltage diode model has been presented which utilises the fact that the diode's parameters depend on instantaneous values of junction temperature and forward current. Diode's voltage is dependent on its current and the junction temperature and is very easy to measure. It is therefore possible to create a model of diode's colorimetric and power properties based solely on instantaneous values of diode's current and voltage. The model was proved to generate accurate colorimetric feedback much under just noticeable $ΔE_{ab}^{*}$ colour distance even under pulsed current conditions.
  2612. The current-voltage model can be applied to string connected light-emitting diodes but the accuracy of the method has not yet been verified in laboratory.
  2614. The model can be easily used in trichromatic luminaire colour control, where the only needed feedback value is the diode voltages. The instantaneous current information is taken from current command as the two values are different only during very short current settling transients.
  2616. \marginpar{detailed model of a luminaire}Together with detailed thermal and electrical model of a diode, the model can provide a good platform for simulating colorimetric and power properties of the diode under different current shapes mentioned before. Current-voltage model can provide information about dissipated power to the detailed thermal model. The resulting junction temperature and driving current can be converted into forward voltage used in the current-voltage model. The accuracy of this complete luminaire description has not yet been proven and is a part of suggested future work. The detailed model has been used to show the effect of heastink thermal resistance on maximal achievable flux.
  2618. \marginpar{model generation procedure}The data acquisition procedure used to create the current-voltage model is far from perfect as it includes the delays for thermally stabilising the luminaire system. As part of future work, the fact that the thermal time constant of \led{} structure and heatsink are much different should be utilised to increase the speed of the model generation. Also the means of generating the model of light-emitting diodes without the use of an actively controlled heatsink should be investigated.
  2620. \marginpar{power converter}The review of diodes' behaviour under different dimming mechanisms led to the conclusion that driving the diode with \dc{} current yielded highest efficiency compared to other pulsed current methods. Therefore the converter should supply the diode with constant current of regulated value. Previous research also shows importance of lifetime analysis of the converter. The key component, limiting the lifetime of the converter, is assumed to be the electrolytic capacitor. As electrolytic capacitors are typically used in input and output of power converters as part of the filters, a topology minimising the need for these capacitors was investigated. The dual interleaved buck converter was built and tested driving four series connected high current \textsc{cbt-90} diodes. This converter minimises the output current ripple using interleaving technique, where the current ripple from both phases cancel each other. Around 0.5 duty cycle this effect is the highest. Because of low series resistance of the light-emitting diodes, their voltage does not change much with the change of the driving current. The dual interleaved converter can be used with input voltage close to twice the forward voltage of the diode, therefore operating close to the 0.5 duty cycle in the whole diode current dimming range and thus minimising the need for output capacitor.
  2622. % chapter conclusions (end)
  2624. %%%%%%%%%%%%%%%%
  2625. % Bibliography %
  2626. %%%%%%%%%%%%%%%%
  2627. \cleardoublepage
  2628. \phantomsection
  2629. \addcontentsline{toc}{chapter}{\numberline{ }Bibliography}
  2630. \bibliographystyle{unsrt}
  2631. \bibliography{thesis}
  2632. \end{document}