<?xml version="1.0" encoding="UTF-8"?>
<!--
 * Scilab ( https://www.scilab.org/ ) - This file is part of Scilab
 * Copyright (C) 2008 - INRIA
 *
 * Copyright (C) 2012 - 2016 - Scilab Enterprises
 *
 * This file is hereby licensed under the terms of the GNU GPL v2.0,
 * pursuant to article 5.3.4 of the CeCILL v.2.1.
 * This file was originally licensed under the terms of the CeCILL v2.1,
 * and continues to be available under such terms.
 * For more information, see the COPYING file which you should have received
 * along with this program.
 *
 -->
<refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink"
          xmlns:svg="http://www.w3.org/2000/svg" xmlns:ns5="http://www.w3.org/1999/xhtml"
          xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook"
          xmlns:scilab="http://www.scilab.org" xml:id="diag" xml:lang="fr">
    <refnamediv>
        <refname>diag</refname>
        <refpurpose>extraction de diagonale ou construction d'une matrice avec un vecteur en diagonal
        </refpurpose>
    </refnamediv>
    <refsynopsisdiv>
        <title>Séquence d'appel</title>
        <synopsis>
           y = diag(vm)
           y = diag(vm, k)
        </synopsis>
    </refsynopsisdiv>
    <refsection>
        <title>Paramètres</title>
        <variablelist>
            <varlistentry>
                <term>vm</term>
                <listitem>
                    <para>vecteur ou matrice (stockage plein ou creux)</para>
                </listitem>
            </varlistentry>
            <varlistentry>
                <term>k</term>
                <listitem>
                    <para>entier (valeur par défaut 0)</para>
                </listitem>
            </varlistentry>
            <varlistentry>
                <term>y</term>
                <listitem>
                    <para>vecteur ou matrice</para>
                </listitem>
            </varlistentry>
        </variablelist>
    </refsection>
    <refsection>
        <title>Description</title>
        <para>
            Pour <literal>vm</literal> un n-vecteur (ligne ou colonne)
            <literal>diag(vm)</literal> renvoie une matrice diagonale avec les
            composantes de <literal>vm</literal> sur la diagonale principale.
        </para>
        <para>
            <literal>diag(vm,k)</literal> est une matrice
            <literal>(n+abs(k))x(n+abs(k))</literal> avec les composantes de
            <literal>vm</literal> sur la kième diagonale. <literal>k=0</literal>
            correspond à la diagonale principale <literal>k&gt;0</literal> correspond
            aux diagonales supérieures et <literal>k&lt;0</literal> aux diagonales
            inférieures.
        </para>
        <para>
            Pour une matrice <literal>vm</literal>,
            <literal>diag(vm,k)</literal> est le vecteur colonne dont les composantes
            sont les termes de la kième diagonale de <literal>vm</literal>.
            <literal>diag(vm)</literal> est la diagonale principale de
            <literal>vm</literal>. <literal>diag(diag(x))</literal> est une matrice
            diagonale.
        </para>
        <para>
            Si <literal>vm</literal> est une matrice creuse
            <literal>diag(vm,k)</literal> renvoie une matrice creuse.
        </para>
        <para>Pour construire un système linéaire diagonal, utilisez
            <literal>blockdiag</literal>.
        </para>
        <para>
            Notez que <literal>eye(A).*A</literal> renvoie une matrice diagonale
            avec les termes diagonaux de <literal>A</literal>. Ceci est valide pour
            tout type de matrice (constante, polynomiale, rationnelle, système
            dynamique linéaire, ...).
        </para>
    </refsection>
    <refsection>
        <title>Exemples</title>
        <programlisting role="example"><![CDATA[
diag([1,2])

A=[1,2;3,4];
diag(A)  // diagonale principale
diag(A,1)

diag(sparse(1:10))  // matrice diagonale creuse

// Construction d'une matrice tridiagonale de taille 2*m+1
m=5;diag(-m:m) +  diag(ones(2*m,1),1) +diag(ones(2*m,1),-1)
 ]]></programlisting>
    </refsection>
    <refsection role="see also">
        <title>Voir aussi</title>
        <simplelist type="inline">
            <member>
                <link linkend="toeplitz">toeplitz</link>
            </member>
            <member>
                <link linkend="eye">eye</link>
            </member>
            <member>
                <link linkend="blockdiag">blockdiag</link>
            </member>
            <member>
                <link linkend="spec">spec</link>
            </member>
            <member>
                <link linkend="bdiag">bdiag</link>
            </member>
        </simplelist>
    </refsection>
</refentry>
