next up previous contents adam-math, MATmatcc, MMKOMP MMkomp Adam WILEY-VCH Prüfungen zum MATLAB Kurs MATLAB-Filme Projekte und M-Files Math Topics Stefan Adam persönlich Mail to stefan.adam@psi.ch
Nächste Seite: SS 07 - Prüfung Aufwärts: Sommersemester 2007 Vorherige Seite: SS07 - Prüfung 1,   Inhalt

SS 07 - Prüfung 1, Lösungen, R-G-B-Y 30. Mai 2007

R   Ingenieurmathematik Prüfung 1 30.Juni2007
Zeit 90 Minuten, Reihenfolge beliebig, 8 Punkte pro Hauptaufgabe, 40 Pt. = N.6.

1)
Verständnisfragen: Es werden nur ganz kurze Antworten erwartet.

1a)
Geben Sie je eine MATLAB-Bibliotheksfunktionen an, welche die Signatur -komplexe Eingabe $\rightarrow$ reelles Resultat bzw. - komplexe Eingabe $\rightarrow$ komplexes Resultat aufweisen!

L1a)
c-r: real(), imag(), abs(), angle(), c-c: conj(), exp()

1b)
Welche Bedingungen, falls überhaupt Bedingungen nötig sind, muss eine Matrix B erfüllen,
im Fall a) dass das Produkt $B \cdot B$ legal ist, und
im Fall b) dass das Produkt $B^T \cdot B$ legal ist?

L1b)
a) B quadratisch; b) keine Bedingung

1c)
Wieviele frei wählbare Zahlenwerte sind in einer nxn Matrix enthalten, falls diese eine oberere Dreicksmatrix ist.

L1c)
n*(n+1)/2

1d)
Wie erreicht man, dass in einer MATLAB-Grafik die Funktionen x=cos(w) und y=sin(w) einen wirklich runden Kreis produzieren?

L1d)
axis equal oder axis ([-1 1 -1 1]) zusammen mit axis square

2)
Suchen Sie die speziellen Permutations/Auswahlmatrizen $Pl$ und $Pr$, so dass die folgende Matrizengleichung für beliebige Werte der Matrix A gilt!
$\displaystyle{
\left(
\begin{array}{rrrrr}
e_4 & e_1 & 0 & e_5 & 0 \\
c...
..._2 & e_3 & e_4 & e_5
\end{array}
\right) \cdot Pr ~=~ Pl \cdot A \cdot Pr
}$

L2)
 RPl = zeros(5) ; RPr = zeros(5); A=iwmat(5)
 RPl(1,5) = 1;   RPl(2,3) = 1;   RPl(4,2) = 1;  
 RPr(4,1) = 1;   RPr(1,2) = 1;   RPr(5,4) = 1;
 As = RPl*A*RPr, RPl, RPr

3)
Suchen sie alle komplexen Lösungen der Gleichung

\begin{displaymath}
z^3 +2\cdot \sqrt{2} = 0
\end{displaymath}

L3)

\begin{displaymath}
z^3 = - 2\cdot \sqrt{2} = \sqrt{2}^3 \cdot \exp{j\cdot pi}
\end{displaymath}


\begin{displaymath}
z_k = \sqrt{2} \cdot \exp{j\cdot pi/3 + k \cdot 2\pi/3 ~~~ k=0 \ldots 2}
\end{displaymath}

4)
Gegeben ist der Punkt $A (8/ 6)$ in der Ebene. Bestimmen Sie die Geradengleichung in der Hesse'schen Normalform für die beiden Geraden:
1. Gerade g durch den Koordinatenursprung O und den Punkt A
2. Gerade h , senkrecht zu g durch den Punkt A
Bestimmen Sie zudem den Schnittpunkt der Geraden h mit der y-Achse (indem Sie in der zu h gehörenden Geradengleichung den x-Koordinatenwert des allgemeinen Punktes Null setzen und dessen y-Koordinatenwert suchen)!

L4)
OA = [8; 6] ;
N = [-6 ; 8]; eng = N/norm(N) , dkg = eng'*OA % dkg muss = 0 sein
% eng = [-0.6 0.8]  g: eng'*OP - 0 = 0
% N orth OA ; enh orth (orth OA) , also parallel OA
enh = OA/norm(OA), dkh = enh'*OA
% enh = [0.8 0.6]  h: enh'*OP - 10 = 0 
yh = dkh/enh(2)   % = 16.66

5)
Bestimmen Sie die folgenden zwei archimedischen Spiralen:
1) Die sich im Gegenurzeigersinn (mathematisch positive Winkeländerung) öffnende Spirale durch die Punkte $Pp(-2/0)$ und $Qp(0/-4)$ und
2) Die sich im Urzeigersinn öffnende Spirale durch die Punkte $Pn(2/0)$ und $Qn(0/-4)$.
Geben Sie die MATLAB-Befehle an zum Zeichnen der beiden Zweige je mit einer halben Drehung zwischen der +y und der -y-Achse!

L5)
% 2 = an*(0 - wn0) ; 4 = an*(-pi/2 - wn0) ;
%  diff:  2 = an*(-pi/2) ; an = -4/pi; wn0 = pi/2
an = -4/pi ; wn0 = pi/2; wn = (0.5:-0.01:-0.5)*pi ;
xn = an*(wn - wn0).*cos(wn); yn = an*(wn - wn0).*sin(wn); 
plot(xn, yn)
axis([-5 5 -5 5]); axis equal; hold on
plot([0 0],[-5 5],'k'); plot([-5 5],[0 0],'k'); 
% 2 = ap*(pi - wp0) ; 4 = ap*(3*pi/2 - wp0) ;
%  diff:  2 = ap*(pi/2) ; ap = 4/pi; wp0 = 2/(4/pi) = pi/2
ap = 4/pi ; wp0 = pi/2; wp = (0.5:0.01:1.5)*pi ;
xp = ap*(wp - wp0).*cos(wp); yp = ap*(wp - wp0).*sin(wp); 
plot(xp, yp,'r')
hold off

6)
Schreiben Sie eine MATLAB-Funktion, welche eine eingegebene Matrix daraufhin prüft, ob sie eine obere Dreiecksmatrix ist, und je nach dem Ausgang der Prüfung true oder false zurückgibt. (Dazu muss vorher getestet werden ob sie quadratisch ist.)

L6
function iftri =  utritestb(M)
% Test auf obere Dreiecksmatrix
% 1. ist M quadratisch?
  [nzei,nspa] = size(M);
  if nzei ~= nspa
    iftri = false;
    return
  else
%  unteres Feld muss Nullen enthalten
    iftri = true;
    for zei = 2:nzei
      for spa = 1:zei-1
        if M(zei,spa) ~= 0
          iftri = false;
          return
        end
      end
    end    
  end


next up previous contents adam-math, MATmatcc, MMKOMP MMkomp Adam WILEY-VCH Prüfungen zum MATLAB Kurs MATLAB-Filme Projekte und M-Files Math Topics Stefan Adam persönlich Mail to stefan.adam@psi.ch
Nächste Seite: SS 07 - Prüfung Aufwärts: Sommersemester 2007 Vorherige Seite: SS07 - Prüfung 1,   Inhalt
2012-03-21