Tråden för oss som använder Matlab

Inga problem, gl !

Tjenare

Slitit mig i håret i två dagar strax och tänkte att det får fan räcka att försöka själv.
Har till uppgift att skapa en cirkel med hjälp av två inputs, den ena är arean för cirkeln och den andra är avståndet som cirkeln skall rotera runt från centrum på figuren. Alltså (0,0) och inputen är 5 för rotationsradien så skall cirkeln befinna sig 5 steg utåt. Här är den kod jag kommit fram med hittills men får fan inte skiten att rotera runt den radie jag själv sätter utan den snurrar runt sin egen satta cirkelpunkt. Att den bara snurrar runt sin egen axel kan jag förstå någorlunda på grund av cos och sin funktionen men vet inte hur jag ska lösa problemet.

function Lab9a_1(a,rotRad)
close all
[x,y] = meshgrid(-200:200);
circleRadius = a/pi;
radius = sqrt(circleRadius);

for j = 1:30

cx = x + cos(j*2*pi/30);
cy = y + sin(j*2*pi/30);
drawDisk(cx-rotRad,cy-rotRad,radius);
mov(j) = getframe;
end
movie(mov,10);
end

Min hjälpfunktion:
function drawDisk(cx,cy,radius)
[x,y] = meshgrid(-200:200);
circleImage = (cx).^2 + (cy).^2 < radius^2; % the disk
imagesc(x(:,1),y(:,1),circleImage);

Tack på förhand och trevlig eftermiddag!

Tjena! Utgick inte så mycket från ditt lösningsförslag, utan från din problembeskrivning, så jag kan ha missat något. Gjorde den på ett lite annorlunda sätt dessutom. Men du kan ju kolla om det är något som skulle vara av intresse.

Btw, har du inte formulerat din radie lite fel? Om din variabel a var arean, så borde väl radien vara sqrt(a/pi) ?

%--------------

clc, clear variables, clf, close all
d = 4;
a = 12.5664;
r = sqrt(a/pi)

theta = linspace(0,2*pi);

% hold on
for i = 1:360

x = r.*cos(theta)+d.*cosd(i); %cosd och sind är samma, bara att
y = r.*sin(theta)+d.*sind(i); %deras argument är i grader istället för rad

plot(x,y)
axis([-10 10 -10 10])

pause(0.01)
end

%--------------

Ja det är ju så funktionen skall vara, men skall skapa en bild och göra en film utav det med funktionerna imagesc och movie etc.

Löste problemet, här kommer slutgiltiga koden.

function Lab9a_1(a,rotRad)
close all
circleRadius = sqrt(a/pi);

cx = 0;
cy = 0;

for j = 1:30
cx = rotRad*cos(j*2*pi/30);
cy = rotRad*sin(j*2*pi/30);
drawDisk(cx,cy,circleRadius);
mov(j) = getframe;
end
movie(mov,11,24);
end

function drawDisk(cx,cy,radius)
[x,y] = meshgrid(-200:200);
circleImage = (x-cx).^2 + (y-cy).^2 < radius^2; % the disk
imagesc(x(:,1),y(:,1),circleImage);

När man använder [V,D] = eig(A) får man ju som bekant en diagonal matris D och en matris V med normerade egenvektorer, som tillhör A's egenrum.

Kan man få icke-normerade egenvektorer istället, så man slipper decimaltalen?

Jag tror inte man kan göra det, det är väl inte säkert att alla egenvektor kan representeras som heltal?

Är väl bara skala med lämpligt tal? Om egenrummet som svarar mot λ är n-dimensionellt så innebär det ju att lösningen till (A - λI)v = 0 har n stycken fria variabler (t1, ... ,tn) som kan väljas lämpligt och på så viss få heltalsvektorer.

Min fråga är då om det finns något sätt att från dom normerade vektorerna skala om dom till heltalsvektorer.

Just nu är e1 vinkelrät mot e2 och e3, men e2 = -e3. Vad kan vara fel?

Lat T beteckna tetraedern vars horn ar punkterna (0; 0; 0), (a; 0; 0), (0; b; 0) och
(0; 0; c), dar a; b; c > 0.
(a) Bestam ekvationerna for de plan som utgor T:s randytor. (Denna uppgift kan
med fordel goras hemma.)
(b) Skriv en matlabfunktion vars forsta rad ar:
function [d n]=tetradist(x,y,z,a,b,c)
Om punkten (x; y; z) ligger inuti T ska de returnerande variablerna d och
n vara avstandet fran (x; y; z) till randen respektive enhetsnormalvektorn i
narmaste randpunkt som pekar mot (x; y; z). Om (x; y; z) ligger utanfor eller
pa randen till T returneras avstandet till det narmsta av tetraederns horn
(d) samt koordinaterna for det hornet (n).
2 (4)
(c) Visa att din kod verkligen fungerar genom att valja a = 1, b = 2, c = 3 och
minst tre punkter (x; y; z) varav en ligger inuti T, en ligger pa randen till T
och en ligger utanfor T.
Tips: Du kan anvanda foljande kod for att slumpa fram en punkt inuti T.
x=rand(1)*a, y=rand(1)*b*(1-x/a), z=rand(1)*c*(1-x/a-y/b)



Hjälp tack !!

hur långt har du kommit, jag kommer bara ge dig svaret om du inte har en specik fråga.

Liten fråga:

>> u = [8 3 -1 7];
>> x = linspace(-45,45);
>> v = polyval(u,x);
>> plot(x,v)

u är alltså vektorn som beskriver polynomet 8x^2 + 3x -x + 7

Fråga: Hur beskriver man sin(pi*x) i en vektoroperation? Eller för den delen e^x ?
Alltså om man då har sin(pi*x) + e^x hur avbildar man det polynomet som en vektor som i första exemplet?