can anyone help me run this, or get this running
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
Relly Syam
am 2 Aug. 2021
Beantwortet: Walter Roberson
am 2 Aug. 2021
clear; clc;
format long e
tic %bagian
n=10 ;
syms c0 c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 c11 t r0 r1 r2 r3 r4 r5 r6 r7 r8 r9 r10 r11 c0=0;
c1=1/11; c2=2/11; c3=3/11; c4=4/11; c5=5/11; c6=6/11; c7=7/11; c8=8/11; c9=9/11; c10=10/11; c11=1;
r0=0; r1=1; r2=2; r3=3; r4=4; r5=5; r6=6; r7=7; r8=8; r9=9; r10=10; r11=11;
EvalAt = [c0, c1, c2, c3, c4, c5, c6, c7, c8, c9, c10, c11];
ktemp = arrayfun(@(EA) euler([r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11], EA).', EvalAt, 'uniform', 0);
ptemp=arrayfun(@(EA) int((3^(EA-t))*euler([r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11],t),t,0,EA).', EvalAt, 'uniform', 0);
E = horzcat(ktemp{:}).' ;
K = horzcat(ptemp{:}).' ;
Ek=E-K ;
Inv_Ek= inv(Ek) ;
F=[c0*3^c0;c1*3^c1;c2*3^c2;c3*3^c3;c4*3^c4;c5*3^c5;c6*3^c6;c7*3^c7;c8*3^c8;c9*3^c9;c10*3^c10;c11*3^c11];
C=Ek\F ;
Ua=@(x)(C(1)*euler(0,x)+C(2)*euler(1,x)+C(3)*euler(2,x)+C(4)*euler(3,x)+C(5)*euler(4,x)+C(6)*euler(5,x)+C(7)*euler(6,x)+C(8)*euler(7,x)+C(9)*euler(8,x)+C(10)*euler(9,x)+C(11)*euler(10,x)+C(12)*euler(11,x)) ;
Ue=@(x)((3^EA)*(1-e^(-EA)));
uaa=zeros(12,1) ;
uee=zeros(12,1) ;
xx=zeros(12,1) ;
k=0;
for i=1:12
uaa(i)=Ua(k);
uee(i)=Ue(k);
xx(i)=k;
k=k+.1;
end
Uap= uaa;
Uex= uee;
y=(abs(uaa-uee));
[xx uee uaa y];
uee uaa y
%gambar plot(xx,uee,'k',xx,uaa,':o') grid on legend({'eksak','aproksimasi N=10'},'Location','Northwest')
%plot(xx,y,'k') %legend({'error N=10'},'Location','Northwest') %grid on xlabel('x'),ylabel('Eksak, Aproksimasi'); toc
0 Kommentare
Akzeptierte Antwort
Walter Roberson
am 2 Aug. 2021
format long e
n=10 ;
syms c0 c1 c2 c3 c4 c5 c6 c7 c8 c9 c10 c11 t r0 r1 r2 r3 r4 r5 r6 r7 r8 r9 r10 r11
c0=0;
c1=1/11; c2=2/11; c3=3/11; c4=4/11; c5=5/11; c6=6/11; c7=7/11; c8=8/11; c9=9/11; c10=10/11; c11=1;
r0=0; r1=1; r2=2; r3=3; r4=4; r5=5; r6=6; r7=7; r8=8; r9=9; r10=10; r11=11;
EvalAt = [c0, c1, c2, c3, c4, c5, c6, c7, c8, c9, c10, c11];
ktemp = arrayfun(@(EA) euler([r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11], EA).', EvalAt, 'uniform', 0);
ptemp=arrayfun(@(EA) int((3^(EA-t))*euler([r0, r1, r2, r3, r4, r5, r6, r7, r8, r9, r10, r11],t),t,0,EA).', EvalAt, 'uniform', 0);
E = horzcat(ktemp{:}).' ;
K = horzcat(ptemp{:}).' ;
Ek=E-K ;
%do not calculate inv(Ek) !!
%{
Inv_Ek= inv(Ek) ;
%}
F=[c0*3^c0;c1*3^c1;c2*3^c2;c3*3^c3;c4*3^c4;c5*3^c5;c6*3^c6;c7*3^c7;c8*3^c8;c9*3^c9;c10*3^c10;c11*3^c11];
C = double(Ek)\F ; %<=== important change!!
Ua=@(x)(C(1)*euler(0,x)+C(2)*euler(1,x)+C(3)*euler(2,x)+C(4)*euler(3,x)+C(5)*euler(4,x)+C(6)*euler(5,x)+C(7)*euler(6,x)+C(8)*euler(7,x)+C(9)*euler(8,x)+C(10)*euler(9,x)+C(11)*euler(10,x)+C(12)*euler(11,x)) ;
Ue = @(EA)((3^EA)*(1-exp(-EA))); %<=== recheck this !!!
uaa=zeros(12,1) ;
uee=zeros(12,1) ;
xx=zeros(12,1) ;
k=0;
for i=1:12
uaa(i)=Ua(k);
uee(i)=Ue(k);
xx(i)=k;
k=k+.1;
end
Uap= uaa;
Uex= uee;
y=(abs(uaa-uee));
[xx, uee, uaa, y];
uee, uaa, y
%gambar plot(xx,uee,'k',xx,uaa,':o') grid on legend({'eksak','aproksimasi N=10'},'Location','Northwest')
%plot(xx,y,'k') %legend({'error N=10'},'Location','Northwest') %grid on xlabel('x'),ylabel('Eksak, Aproksimasi'); toc
0 Kommentare
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Time Series finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!