%in the name of God hold on r =7; turn = 1000; tetha = 0:0.03:2*turn*pi; speed = s1 = 1.0005; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'color', [2/255, 140/255, 207/255]); turn = randi(70) tetha = 0:0.005:2*turn*pi; r1 = randi(15); r2 = randi(200); speed = s2 = r1 + 1 / r2; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'color', [143/255, 210/255, 158/255]); turn = randi(70) tetha = 0:0.005:2*turn*pi; r1 = randi(15); r2 = randi(200); speed = s3 = r1 + 1 / r2; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'w'); ylabel( ['https://nooh124.blog.ir ' ,num2str(s2),' ' , num2str(s3)] ) axis equal; axis([-14.5 14.5 -14.5 14.5]) print -dpng figure.png
clc format short disp (' **********************************************************') disp (' *** Ajustamento de uma linha de nivelamento geometrico ***') disp (' **********************************************************') disp (' **********************************************************') disp (' ******************percurso 1 DINI22***********************') disp (' **********************************************************') pontos=cellstr(['T1';'P7';'C';'P';'B';'A';'I1';'I2';'I3';'I4';'I5';'I6';'I7';'I8';'I9';'I10';'M3']); %cota_T1+desnivel_T1P7+residuo_T1P7-cota_P7=0 %cota_P7+desnivel_P7C+residuo_P7C-cota_C=0 %cota_C+desnivel_CP+residuo_CP-cota_P=0 %cota_P+desnivel_PB+residuo_PB-cota_B=0 %cota_B+desnivel_BA+residuo_BA-cota_A=0 %cota_A+desnivel_AI1+residuo_AI1-cota_I1=0 % ... %cota_I10+desnivel_I10M3+residuo_I10M3-cota_M3=0 % % %Os pontos começados com I são pontos intermédios A=zeros(16,15); for i=1:15 A(i,i) = -1; A(i + 1,i) = 1; end W(1,1)= 80.108-1.8486; W(2,1)=-0.0569; W(3,1)=-0.2131; W(4,1)=0.5814; W(5,1)=-0.1922; W(6,1)=-0.3585; W(7,1)=-0.5723; W(8,1)=-0.2508; W(9,1)=-0.278; W(10,1)=-0.5356; W(11,1)=-0.3216; W(12,1)=-0.1517; W(13,1)=-0.3064; W(14,1)=-0.2759; W(15,1)=-0.2911; W(16,1)=0.2109-75.23; % % % dist(1,1)=53.85; dist(2,1)=49.69; dist(3,1)=30.23; dist(4,1)=37.62; dist(5,1)=57.45; dist(6,1)=70.28; dist(7,1)=54.59; dist(8,1)=56.10; dist(9,1)=49.03; dist(10,1)=60.47; dist(11,1)=65.17; dist(12,1)=57.42; dist(13,1)=59.88; dist(14,1)=60.22; dist(15,1)=56.47; dist(16,1)=9.73; % %Dini22 "Standard deviation on 1 km of double levelling" + "foldable bar code staff" %sigma_Dini22 = 1.3mm % %Dini22 "Accuracy of distance measuremen" + "foldable bar code staff" %sigma_Dini22 = 30mm % sigma_dini22=30/1000; sigma=zeros(16,16); for i=1:16 sigma(i,i)=sigma_dini22*sqrt(2*dist(i,1)/1000); end sigma2=sigma^2; sigma02=1; QL=sigma2/sigma02; P=inv(QL); % %*************AJUSTAMENTO************************** % At= A.'; N=At*P*A; U=At*P*(-W); % % parametros ajustados % X=inv(N)*U cotas=X; % % matriz cofactor dos parametros % Qx=inv(N); % % PRECISAO DAS COTAS % for i=1:15 precisao_cotas(i,1)=sqrt(Qx(i,i)); end % % RESIDUOS % v=A*X+W; % % % Qv=(inv(P)-A*inv(N)*At); % % OBSERVACOES AJUSTADAS % La(1,1)=-1.8486+v(1,1); La(2,1)=-0.0569+v(2,1); La(3,1)=-0.2131+v(3,1); La(4,1)=0.5814+v(4,1); La(5,1)=-0.1922+v(5,1); La(6,1)=-0.3585+v(6,1); La(7,1)=-0.5723+v(7,1); La(8,1)=-0.2508+v(8,1); La(9,1)=-0.278+v(9,1); La(10,1)=-0.5356+v(10,1); La(11,1)=-0.3216+v(11,1); La(12,1)=-0.1517+v(12,1); La(13,1)=-0.3064+v(13,1); La(14,1)=-0.2759+v(14,1); La(15,1)=-0.2911+v(15,1); La(16,1)=0.2109+v(16,1); La % % % QLa=inv(P)-Qv; % % variancia a posteriori por unidade de peso % s02=v.'*P*v/(16-15); % distribuicao do qui-quadrado a 0.05% qui(1,1)=3.841; qui(2,1)=5.991; qui(3,1)=7.815; qui(4,1)=9.488; qui(5,1)=11.071; qui(6,1)=12.592; qui(7,1)=14.067; qui(8,1)=15.507; qui(9,1)=16.919; qui(10,1)=18.307; qui(11,1)=19.675; qui(12,1)=21.026; qui(13,1)=22.362; qui(14,1)=23.685; qui(15,1)=24.996; qui(16,1)=26.296; qui(17,1)=27.587; qui(18,1)=28.869; qui(19,1)=30.144; qui(20,1)=31.410; qui(21,1)=32.671; qui(22,1)=33.924; qui(23,1)=35.172; qui(24,1)=36.415; qui(25,1)=37.652; qui(26,1)=38.885; qui(27,1)=40.113; qui(28,1)=41.337; qui(29,1)=42.557; qui(30,1)=43.773; qui(31,1)=44.985; qui(32,1)=46.194; qui(33,1)=47.400; qui(34,1)=48.602; qui(35,1)=49.802; qui(36,1)=50.998; qui(37,1)=52.192; qui(38,1)=53.384; qui(39,1)=54.572; qui(40,1)=55.758; qui(41,1)=56.942; qui(42,1)=58.124; qui(43,1)=59.304; qui(44,1)=60.481; qui(45,1)=61.656; % teste da razao de variancias y=s02*(4-3)/sigma02; if (y < qui(4-3,1)) % sigma02 e s02 estatisticamente semelhantes 'Aceitar o ajustamento.' else 'Rejeitar o ajustamento.' end
%in the name of God hold on r =7; turn = 1000; tetha = 0:0.03:2*turn*pi; speed = s1 = 1.0005; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'color', [210/255, 30/255, 47/255]); turn = 70; tetha = 0:0.005:2*turn*pi; r1 = randi(15); r2 = randi(200); speed = s2 = r1 + 1 / r2; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'color', [143/255, 210/255, 158/255]); turn = 70; tetha = 0:0.005:2*turn*pi; r1 = randi(15); r2 = randi(200); speed = s3 = r1 + 1 / r2; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'w'); ylabel( ['https://nooh124.blog.ir ' ,num2str(s2),' ' , num2str(s3)] ) axis equal; axis([-14.5 14.5 -14.5 14.5]) print -dpng figure.png
%in the name of God hold on r =7; turn = 1000; tetha = 0:0.03:2*turn*pi; speed = s1 = 1.0005; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'color', [210/255, 30/255, 47/255]); turn = 70; tetha = 0:0.005:2*turn*pi; r1 = randi(15); r2 = randi(200); speed = s2 = r1 + 1 / r2; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'color', [143/255, 210/255, 158/255]); turn = 70; tetha = 0:0.005:2*turn*pi; r1 = randi(15); r2 = randi(200); speed = s3 = s2 + 5; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'w'); ylabel( ['https://nooh124.blog.ir ' ,num2str(s2),' ' , num2str(s3)] ) axis equal; axis([-14.5 14.5 -14.5 14.5]) print -dpng figure.png
%in the name of God hold on r =7; turn = 1000; tetha = 0:0.03:2*turn*pi; speed = s1 = 1.0005 x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'color', [210/255, 30/255, 47/255]); turn = 70; tetha = 0:0.005:2*turn*pi; r1 = randi(15); r2 = randi(200); speed = s2 = r1 + 1 / r2; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'color', [143/255, 210/255, 158/255]); turn = 70; tetha = 0:0.005:2*turn*pi; r1 = randi(15); r2 = randi(200); speed = s3 = s2 / r1; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'w'); ylabel( ['https://nooh124.blog.ir ' ,num2str(s2),' ' , num2str(s3)] ) axis equal; axis([-14.5 14.5 -14.5 14.5]) print -dpng figure.png
%in the name of God hold on r =7; turn = 1000; tetha = 0:0.03:2*turn*pi; speed = s1 = 1.0005 x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'color', [210/255, 30/255, 47/255]); turn = 70; tetha = 0:0.005:2*turn*pi; r1 = randi(15); r2 = randi(200); speed = s2 = r1 + 1 / r2; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'color', [143/255, 210/255, 158/255]); turn = 70; tetha = 0:0.005:2*turn*pi; r1 = randi(15); r2 = randi(200); speed = s3 = r1 + 1 / r2; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'w'); ylabel( ['https://nooh124.blog.ir ' ,num2str(s2),' ' , num2str(s3)] ) axis equal; axis([-14.5 14.5 -14.5 14.5]) print -dpng figure.png
x = [1 2 3 4 5 6 7 8 9 10]; y1 = [.16 .08 .04 .02 .013 .007 .004 .002 .001 .0008 ]; y2 = [.16 .07 .03 .01 .008 .003 .0008 .0003 .00007 .00002 ]; semilogy(x,y1,'-bo;y1;',x,y2,'-kx;y2;'); title('Plot title'); xlabel('X Axis'); ylabel('Y Axis'); print -dpng figure.png
AM = input("\nEnter the amplitude of modulating signal"); fm = input("\nEnter the frequency of modulating signal"); Tm = 1/fm; t = 0: Tm/100: 3*Tm; M = AM * cos(2*pi*fm*t); %Formula for modulating signal subplot(3, 1, 1); plot(t, M); %plot message signal title('Message Signal'); xlabel('Time'); ylabel('Amplitiude'); AC = input("\nEnter the amplitude of carrier signal"); fc = input("\nEnter the frequency of carrier signal"); Tc = 1/fc; C = AC * cos(2*pi*fc*t); %Formula for carrier signal subplot(3, 1, 2); plot(t, C); %plot carrier signal title('Carrier Signal'); xlabel('Time'); ylabel('Amplitiude'); k = AM / AC; % Formula for Modulation Index disp("\nThe modulation index is"); disp(k); %Ampitude Modulated(AM) Signal AM_signal = AC * [ 1 + ( k.*cos(2*pi*fm*t ))].* cos( 2*pi*fc*t ); %Formulfor AmplitudeModulated(AM) signal subplot(3, 1, 3); plot(t, AM_signal); %plot AM signal title('AM Signal'); xlabel('Time'); ylabel('Amplitiude'); print -dpng figure.png
%in the name of God hold on r =7; turn = 1000; tetha = 0:0.03:2*turn*pi; speed = s1 = 1.0005 x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'color', [210/255, 30/255, 47/255]); turn = 80; tetha = 0:0.01:2*turn*pi; r1 = randi(15); r2 = randi(200); speed = s2 = r1 + 1 / r2; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'color', [216/255, 219/255, 35/255]); turn = 80; tetha = 0:0.01:2*turn*pi; r1 = randi(15); r2 = randi(200); speed = s3 = r1 + 1 / r2; x = cos(tetha) * r + cos(speed*tetha) * 7; y = sin(tetha) * r + sin(speed*tetha) * 7; plot(x, y, 'w'); ylabel( ['https://nooh124.blog.ir ' ,num2str(s2),' ' , num2str(s3)] ) axis equal; axis([-16.1 16.1 -16.1 16.1]) print -dpng figure.png
We use cookies to provide and improve our services. By using our site, you consent to our Cookies Policy. Accept Learn more