Please note, this is a STATIC archive of website www.tutorialspoint.com from 11 May 2019, cach3.com does not collect or store any user information, there is no "phishing" involved.
Tutorialspoint

beauty1397-12-22

%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

amon

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

DINI22

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

rasm-beautiful-971218-iiii

%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

rasm-beautiful-971218iii

%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

rasm-beautiful-971218ii

%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

Rasm-beautiful-971218i

%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

Execute MATLAB/Octave Online

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

Execute MATLAB/Octave Online

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

cuteShape 2

%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

Advertisements
Loading...

We use cookies to provide and improve our services. By using our site, you consent to our Cookies Policy.