Het grote off-topic topick |
bLAM
Kroquet Bastard Registered: 27 Aug 2001 Location: anywhere Posts: 1669 |
midden "Never let your sense of morals prevent you from doing what's right - Isaac Asimov" |
Gorgeous George
Opmerk Bastard
Registered: 13 Jul 2001 Location: BehindMy21 Posts: 4748
|
wie ziet de fout.
een krat bier voor degene die mijn probleem oplost
function delta=ccc(h,pars)
global Versie Alfaoud F Alfa DeltaX Fcorr FiNU Foud
Fvoorsp = pars(1);
Cbuffer = pars(2);
Cextra = pars(3);
Cveer = pars(4);
DeltaFi = pars(5);
Lslag = pars(6);
Loph = pars(7);
Oafst = pars(8);
Ohoek = pars(9);
Rconst = pars(10);
Xtotaal = pars(11);
Rvar = pars(12);
Alfa = atan(h/Loph)*360/2*pi;
A = sin(0.5*DeltaFi*2*pi/360)*Rvar;
C = Rvar-cos(DeltaFi*2*pi/360)*Rvar;
B = sqrt((2*A)^2-C^2);
D = tan(DeltaFi*2*pi/360)*Loph;
E = C + D;
F = E/sin((DeltaFi-Alfa)*2*pi/360);
G = E/tan((DeltaFi-Alfa)*2*pi/360);
I = Loph*cos(DeltaFi*2*pi/360);
H = B-(Loph-I);
Ecorr = sin(Alfaoud*2*pi/360)*H + cos(Alfaoud*2*pi/360) * E;
Fcorr = Ecorr/sin((DeltaFi-Alfa+Alfaoud)*2*pi/360);
Gcorr = Ecorr/tan((DeltaFi-Alfa+Alfaoud)*2*pi/360);
Hcorr = cos(Alfaoud*2*pi/360) * H - sin(Alfaoud*2*pi/360) * E;
if Foud == 0;
DeltaX = F - (G-H);
else
%DeltaX = Fcorr + Foud - (Gcorr-Hcorr) - Foud;
DeltaX = Foud - (Foud- (Gcorr-Hcorr)) + Fcorr;
end
if Versie == 2 %benadering oude versie
% X2 = ((Rconst-(h/2))*FiNU*2*pi)/360 + sqrt(h^2+Loph^2) - Loph;
X2=(2*pi*Rconst*FiNU)/360 + sqrt(h^2+Loph^2) - Loph;
X1 = (2*pi*Rconst*FiNU)/360;
% F2 = (X2*Cveer + Fvoorsp)*cos(atan(h/Loph));
F2 = (X2*Cveer + Fvoorsp)*cos(atan(h/Loph))+((-Lslag/2)*Cbuffer+Cbuffer*X1+Cextra*X1)*(Rconst/(Rconst - h));;
% h1 = (F2 * Rconst - Fvoorsp * Rconst)/F2;
h1 = (Fvoorsp*Rconst-F2*Rconst)/-F2;
elseif Versie == 3 %simpele model
F2 = (Cveer * (Xtotaal + DeltaX) + Fvoorsp);
h1 = (F2 * Rvar - Fvoorsp * Rconst)/F2;
elseif Versie == 4 %uitbreiding hoekcompensatie en toegevoegde stijfheid
X1 = (2*pi*Rconst*FiNU)/360;
F2 = (Cveer * (Xtotaal + DeltaX) + Fvoorsp)*cos(atan(h/Loph))+((-Lslag/2)*Cbuffer+Cbuffer*X1+Cextra*X1)*(Rconst/(Rvar - h));
h1 = (F2 * Rvar - Fvoorsp * Rconst)/F2;
elseif Versie == 5 %uitbreiding bufferindrukkingscompenatie
X1 = (2*pi*Rconst*FiNU)/360;
FvoorspPap = Fvoorsp/cos(atan(((cos(45*2*pi/360)*Oafst)/cos(Ohoek*2*pi/360))/(((cos(45*2*pi/360)*Oafst)/(cos(Ohoek*2*pi/360))+0.5*Lslag))));
F2 = (Cveer * (Xtotaal + DeltaX) + Fvoorsp)*cos(atan(h/Loph))+((-Lslag/2)*Cbuffer+Cbuffer*X1+Cextra*X1)*(Rconst/(Rvar - h));
F0 = FvoorspPap * cos(atan(((cos(45*2*pi/360)*Oafst)/cos(Ohoek*2*pi/360))/(((cos(45*2*pi/360)*Oafst)/(cos(Ohoek*2*pi/360))+0.5*Lslag-X1))));
h1 = (F2 * Rvar - F0 * Rconst)/F2;
FvoorspPap;
end
delta = h1-h;
clear all
close
%VERSIE: 1=benadering oud model 2=afgelegde weg 3=2+hoekcompensatie 4=3+bufferindrukkingcompensatie
Totalehoek = 90;
global Versie Alfaoud F Alfa DeltaX Fcorr FiNU Foud
DeltaFi = 5; %5
hzplot = zeros ((Totalehoek/DeltaFi),3);
for Versie = 2:1:5
Fvoorsp = 10; %1
Cbuffer = 0; %2
Cextra = 0; %3
Cveer = 0.09; %4
Lslag = 30; %6
Loph = 160; %7
Oafst = 190; %8
Ohoek = 25; %9
Rconst = 40; %10
Foud = 0;
Xtotaal = 0; %11
Rvar =(Rconst * Fvoorsp)/(Fvoorsp - 0.5 * Lslag * Cbuffer); %12
K=1;
%bereken snijpunt voor h=f(h) met h
for Fi = DeltaFieltaFi:Totalehoek %beginpunt:stapgrootte:eindpunt
FiNU = Fi;
pars = [Fvoorsp, Cbuffer, Cextra, Cveer, DeltaFi, Lslag, Loph, Oafst, Ohoek, Rconst, Xtotaal, Rvar];
hz = fsolve('ccc', 0, optimset('fsolve'), pars);
hzplot ((K+1),Versie)=(Rvar-hz);
hzplot (1,Versie)=Rvar;
hzplot (1,1)=40;
hzplot (2,1)=25.62;
hzplot (3,1)=18.86;
hzplot (4,1)=15.09;
hzplot (5,1)=12.65;
hzplot (6,1)=10.93;
hzplot (7,1)=9.64;
hzplot (8,1)=8.64;
hzplot (9,1)=7.83;
hzplot (10,1)=7.16;
hzplot (11,1)=6.6;
hzplot (12,1)=6.13;
hzplot (13,1)=5.72;
hzplot (14,1)=5.36;
hzplot (15,1)=5.04;
hzplot (16,1)=4.77;
hzplot (17,1)=4.52;
hzplot (18,1)=4.29;
hzplot (19,1)=4.09;
K = K + 1;
%opnemen in nieuwe loop Foud, Alfaoud, Xtotaal
if Foud == 0;
Foud = F;
else
Foud = Fcorr;
end
Alfaoud = Alfa;
Xtotaal = Xtotaal + DeltaX;
end
Fi = 0eltaFi:Totalehoek;
Versie
Versie = Versie + 1;
end
fh = plot (Fi',hzplot,'b.-');
set(fh(1),'color',[0,0,1])
set(fh(2),'color',[1,0,1])
set(fh(3),'color',[-0,0,0])
set(fh(4),'color',[0,1,0])
set(fh(5),'color',[0,0.5,0.5])
xlabel ('Fi')
ylabel ('Rvar')
legend('Versie 1 oud model 1.35','Versie 2, simulatie oud model','Versie 3, nieuw model, afgelegde weg','Versie 4, nieuw model, hoekcompensatie','Versie 5, nieuw model, bufferindrukking')
| |
rens
Flying Bastard Registered: 18 Oct 2002 Location: Behind-U Posts: 1257 |
Oh, ik zie de fout al, er staat teveel onzin tussen.
You want some, come get some!!!!!
|
bLAM
Kroquet Bastard Registered: 27 Aug 2001 Location: anywhere Posts: 1669 |
... misschien helpt het als je die smileys eruit gooit
"Never let your sense of morals prevent you from doing what's right - Isaac Asimov" |
Caman
3D Bastard Registered: 11 Apr 2001 Location: up yours Posts: 4721 |
|
fraXTC
Bad-hairday Bastard Registered: 28 Jun 2001 Location: Under_your_skin Posts: 2457 |
|
fraXTC
Bad-hairday Bastard Registered: 28 Jun 2001 Location: Under_your_skin Posts: 2457 |
" http($this->addslashes("http://home.wanadoo.nl/~ingrid.bij/dancebabys/gpalacedb14.gif")) . "\" border=\"0\" alt=\"[user image]\" />"
Electrische bastriangel roeleert
" http("http://www.youtube.com/watch?v=6TUzaW0-K1k") . "\" target=\"_blank\">Wat deze man ook gedaan heeft.. ik vergeef het hem"
|
Fox
The Guzzler Registered: 10 Apr 2001 Location: North-Brabant Posts: 9401 |
|
Rainman
BMW Bastard Registered: 11 Apr 2001 Location: Everywhere Posts: 5391 |
Zo moet je geen last hebben van smilies
function delta=ccc(h,pars)
global Versie Alfaoud F Alfa DeltaX Fcorr FiNU Foud
Fvoorsp = pars(1);
Cbuffer = pars(2);
Cextra = pars(3);
Cveer = pars(4);
DeltaFi = pars(5);
Lslag = pars(6);
Loph = pars(7);
Oafst = pars(8);
Ohoek = pars(9);
Rconst = pars(10);
Xtotaal = pars(11);
Rvar = pars(12);
Alfa = atan(h/Loph)*360/2*pi;
A = sin(0.5*DeltaFi*2*pi/360)*Rvar;
C = Rvar-cos(DeltaFi*2*pi/360)*Rvar;
B = sqrt((2*A)^2-C^2);
D = tan(DeltaFi*2*pi/360)*Loph;
E = C + D;
F = E/sin((DeltaFi-Alfa)*2*pi/360);
G = E/tan((DeltaFi-Alfa)*2*pi/360);
I = Loph*cos(DeltaFi*2*pi/360);
H = B-(Loph-I);
Ecorr = sin(Alfaoud*2*pi/360)*H + cos(Alfaoud*2*pi/360) * E;
Fcorr = Ecorr/sin((DeltaFi-Alfa+Alfaoud)*2*pi/360);
Gcorr = Ecorr/tan((DeltaFi-Alfa+Alfaoud)*2*pi/360);
Hcorr = cos(Alfaoud*2*pi/360) * H - sin(Alfaoud*2*pi/360) * E;
if Foud == 0;
DeltaX = F - (G-H);
else
%DeltaX = Fcorr + Foud - (Gcorr-Hcorr) - Foud;
DeltaX = Foud - (Foud- (Gcorr-Hcorr)) + Fcorr;
end
if Versie == 2 %benadering oude versie
% X2 = ((Rconst-(h/2))*FiNU*2*pi)/360 + sqrt(h^2+Loph^2) - Loph;
X2=(2*pi*Rconst*FiNU)/360 + sqrt(h^2+Loph^2) - Loph;
X1 = (2*pi*Rconst*FiNU)/360;
% F2 = (X2*Cveer + Fvoorsp)*cos(atan(h/Loph));
F2 = (X2*Cveer + Fvoorsp)*cos(atan(h/Loph))+((-Lslag/2)*Cbuffer+Cbuffer*X1+Cextra*X1)*(Rconst/(Rconst - h));;
% h1 = (F2 * Rconst - Fvoorsp * Rconst)/F2;
h1 = (Fvoorsp*Rconst-F2*Rconst)/-F2;
elseif Versie == 3 %simpele model
F2 = (Cveer * (Xtotaal + DeltaX) + Fvoorsp);
h1 = (F2 * Rvar - Fvoorsp * Rconst)/F2;
elseif Versie == 4 %uitbreiding hoekcompensatie en toegevoegde stijfheid
X1 = (2*pi*Rconst*FiNU)/360;
F2 = (Cveer * (Xtotaal + DeltaX) + Fvoorsp)*cos(atan(h/Loph))+((-Lslag/2)*Cbuffer+Cbuffer*X1+Cextra*X1)*(Rconst/(Rvar - h));
h1 = (F2 * Rvar - Fvoorsp * Rconst)/F2;
elseif Versie == 5 %uitbreiding bufferindrukkingscompenatie
X1 = (2*pi*Rconst*FiNU)/360;
FvoorspPap = Fvoorsp/cos(atan(((cos(45*2*pi/360)*Oafst)/cos(Ohoek*2*pi/360))/(((cos(45*2*pi/360)*Oafst)/(cos(Ohoek*2*pi/360))+0.5*Lslag))));
F2 = (Cveer * (Xtotaal + DeltaX) + Fvoorsp)*cos(atan(h/Loph))+((-Lslag/2)*Cbuffer+Cbuffer*X1+Cextra*X1)*(Rconst/(Rvar - h));
F0 = FvoorspPap * cos(atan(((cos(45*2*pi/360)*Oafst)/cos(Ohoek*2*pi/360))/(((cos(45*2*pi/360)*Oafst)/(cos(Ohoek*2*pi/360))+0.5*Lslag-X1))));
h1 = (F2 * Rvar - F0 * Rconst)/F2;
FvoorspPap;
end
delta = h1-h;
clear all
close
%VERSIE: 1=benadering oud model 2=afgelegde weg 3=2+hoekcompensatie 4=3+bufferindrukkingcompensatie
Totalehoek = 90;
global Versie Alfaoud F Alfa DeltaX Fcorr FiNU Foud
DeltaFi = 5; %5
hzplot = zeros ((Totalehoek/DeltaFi),3);
for Versie = 2:1:5
Fvoorsp = 10; %1
Cbuffer = 0; %2
Cextra = 0; %3
Cveer = 0.09; %4
Lslag = 30; %6
Loph = 160; %7
Oafst = 190; %8
Ohoek = 25; %9
Rconst = 40; %10
Foud = 0;
Xtotaal = 0; %11
Rvar =(Rconst * Fvoorsp)/(Fvoorsp - 0.5 * Lslag * Cbuffer); %12
K=1;
%bereken snijpunt voor h=f(h) met h
for Fi = DeltaFieltaFi:Totalehoek %beginpunt:stapgrootte:eindpunt
FiNU = Fi;
pars = [Fvoorsp, Cbuffer, Cextra, Cveer, DeltaFi, Lslag, Loph, Oafst, Ohoek, Rconst, Xtotaal, Rvar];
hz = fsolve('ccc', 0, optimset('fsolve'), pars);
hzplot ((K+1),Versie)=(Rvar-hz);
hzplot (1,Versie)=Rvar;
hzplot (1,1)=40;
hzplot (2,1)=25.62;
hzplot (3,1)=18.86;
hzplot (4,1)=15.09;
hzplot (5,1)=12.65;
hzplot (6,1)=10.93;
hzplot (7,1)=9.64;
hzplot (8,1)=8.64;
hzplot (9,1)=7.83;
hzplot (10,1)=7.16;
hzplot (11,1)=6.6;
hzplot (12,1)=6.13;
hzplot (13,1)=5.72;
hzplot (14,1)=5.36;
hzplot (15,1)=5.04;
hzplot (16,1)=4.77;
hzplot (17,1)=4.52;
hzplot (18,1)=4.29;
hzplot (19,1)=4.09;
K = K + 1;
%opnemen in nieuwe loop Foud, Alfaoud, Xtotaal
if Foud == 0;
Foud = F;
else
Foud = Fcorr;
end
Alfaoud = Alfa;
Xtotaal = Xtotaal + DeltaX;
end
Fi = 0eltaFi:Totalehoek;
Versie
Versie = Versie + 1;
end
fh = plot (Fi',hzplot,'b.-');
set(fh(1),'color',[0,0,1])
set(fh(2),'color',[1,0,1])
set(fh(3),'color',[-0,0,0])
set(fh(4),'color',[0,1,0])
set(fh(5),'color',[0,0.5,0.5])
xlabel ('Fi')
ylabel ('Rvar')
legend('Versie 1 oud model 1.35','Versie 2, simulatie oud model','Versie 3, nieuw model, afgelegde weg','Versie 4, nieuw model, hoekcompensatie','Versie 5, nieuw model, bufferindrukking')
The only reason that you're conscious right now is because I don't want to carry you. - Jack Bauer
" http($this->addslashes("http://87.233.202.7/stuff/bauer/sig.php")) . "\" border=\"0\" alt=\"[user image]\" />" |
Rainman
BMW Bastard Registered: 11 Apr 2001 Location: Everywhere Posts: 5391 |
Dus wel, maar het ziet er al overzichterlijker uit. The only reason that you're conscious right now is because I don't want to carry you. - Jack Bauer
" http($this->addslashes("http://87.233.202.7/stuff/bauer/sig.php")) . "\" border=\"0\" alt=\"[user image]\" />" |
Rainman
BMW Bastard Registered: 11 Apr 2001 Location: Everywhere Posts: 5391 |
Wat voor taal is het eigenlijk? The only reason that you're conscious right now is because I don't want to carry you. - Jack Bauer
" http($this->addslashes("http://87.233.202.7/stuff/bauer/sig.php")) . "\" border=\"0\" alt=\"[user image]\" />" |
Caman
3D Bastard Registered: 11 Apr 2001 Location: up yours Posts: 4721 |
honderd Take my hand and lead me to myself...
" http($this->addslashes("http://www.caman.nl/zooi/sig.jpg")) . "\" border=\"0\" alt=\"[user image]\" />"
|
rens
Flying Bastard Registered: 18 Oct 2002 Location: Behind-U Posts: 1257 |
Oh, ik zie de fout al, er staat teveel onzin tussen. |
such a rookie mistake ... |
Your jealous aren't you?
You want some, come get some!!!!!
|
Gorgeous George
Opmerk Bastard
Registered: 13 Jul 2001 Location: BehindMy21 Posts: 4748
|
het is matlabs
|
bLAM
Kroquet Bastard Registered: 27 Aug 2001 Location: anywhere Posts: 1669 |
matlabs komt van origine uit matlabland waar het wiskundig begaafde volk der matlabianen woont. dit unieke volk werd pas in de laatste helft van de vorige eeuw ontdekt en kentmerkt zich door de enorme wiskundeknobbels op hun hoofden. "Never let your sense of morals prevent you from doing what's right - Isaac Asimov" |
Rainman
BMW Bastard Registered: 11 Apr 2001 Location: Everywhere Posts: 5391 |
Hehe matlabland. En waar ligt dat dan? The only reason that you're conscious right now is because I don't want to carry you. - Jack Bauer
" http($this->addslashes("http://87.233.202.7/stuff/bauer/sig.php")) . "\" border=\"0\" alt=\"[user image]\" />" |
Caman
3D Bastard Registered: 11 Apr 2001 Location: up yours Posts: 4721 |
bij de "m" Take my hand and lead me to myself...
" http($this->addslashes("http://www.caman.nl/zooi/sig.jpg")) . "\" border=\"0\" alt=\"[user image]\" />"
|
fraXTC
Bad-hairday Bastard Registered: 28 Jun 2001 Location: Under_your_skin Posts: 2457 |
|
bLAM
Kroquet Bastard Registered: 27 Aug 2001 Location: anywhere Posts: 1669 |
tussen pythagoria en kwadratica, waar ze overigens dialecten van het matlabs als taal hebben. "Never let your sense of morals prevent you from doing what's right - Isaac Asimov" |
rens
Flying Bastard Registered: 18 Oct 2002 Location: Behind-U Posts: 1257 |
Het Labview volk woont daar toch ook? You want some, come get some!!!!!
|
Rainman
BMW Bastard Registered: 11 Apr 2001 Location: Everywhere Posts: 5391 |
Om het maar niet over die Solidworks schoften te hebben... The only reason that you're conscious right now is because I don't want to carry you. - Jack Bauer
" http($this->addslashes("http://87.233.202.7/stuff/bauer/sig.php")) . "\" border=\"0\" alt=\"[user image]\" />" |
bLAM
Kroquet Bastard Registered: 27 Aug 2001 Location: anywhere Posts: 1669 |
al enkele jaren woedt er een hevige oorlog tussen de inwoners van deze landen en het nabijgelegen Bastardia. "Never let your sense of morals prevent you from doing what's right - Isaac Asimov" |
Caman
3D Bastard Registered: 11 Apr 2001 Location: up yours Posts: 4721 |
paraflexo-idische tegencontracties Take my hand and lead me to myself...
" http($this->addslashes("http://www.caman.nl/zooi/sig.jpg")) . "\" border=\"0\" alt=\"[user image]\" />"
|
bLAM
Kroquet Bastard Registered: 27 Aug 2001 Location: anywhere Posts: 1669 |
paraflexo-idische tegencontracties |
oftewel: oorlog ?
"Never let your sense of morals prevent you from doing what's right - Isaac Asimov" |
Rainman
BMW Bastard Registered: 11 Apr 2001 Location: Everywhere Posts: 5391 |
Gisteren gekeken of ik een e46 M3 bij kon houden. Maar dat viel toch niet mee The only reason that you're conscious right now is because I don't want to carry you. - Jack Bauer
" http($this->addslashes("http://87.233.202.7/stuff/bauer/sig.php")) . "\" border=\"0\" alt=\"[user image]\" />" |