正交配置法 matlab

正交分解法 物体受到多个力作用时求其合力,可将各个力沿两个相互垂直的方向直行正交分解,然后再分别沿这两个方向求出合力,

正交分解法是处理多个力作用用问题的基本方法,值得注意的是,对、方向选择时,尽可能使落在、轴上的力多；

被分解的力尽可能是已知力.步骤为：

①正确选择直角坐标系,一般选共点力的作用点为原点,水平方向或物体运动的加速度方向为X轴,使尽量多的力在坐标轴上.

②正交分解各力,即分别将各力投影在坐标轴上,分别求出坐标轴上各力投影的合力.

clc

clear

close all

Tb = 1

Ns = 10

fs = Ns/Tb% 采样频率

A = 2

t0 = (0:1/fs:Tb-1/fs)

%% 正交信号

s0 = A*ones(1,10) % 正交信号s0

s1 = [A*ones(1,5) -1*A*ones(1,5)] % 正交信号s1

% figure

% subplot(211)

% plot(t0,s0)

% subplot(212)

% plot(t0,s1)

%% 编码

N = 50% 采样点数

code = [1,-1,1,1,-1]% 5个code

for kc = 1:5

if code(kc)==1

x(Ns*(kc-1)+1:kc*Ns) = s0

else

x(Ns*(kc-1)+1:kc*Ns) = s1

end

end

%% 加噪声

% t = 0:1/fs:1/fs*N-1/fs

Nt = 1:N

index2 = 10:10:50

index2 = index2+2

Nt2 = 1:5

deta = [0,0.1,1.0,2.0]

for kd=1:4

noise = sqrt(deta(kd))*randn(1,N) % 噪声

y = x+noise % 输出信号

% y=x

%% 相关

for ks = 1:5

z0(ks) = y(Ns*(ks-1)+1:ks*Ns)*s0'/10

z1(ks) = y(Ns*(ks-1)+1:ks*Ns)*s1'/10

end

% z0 = xcorr(y,s0)

% z1 = xcorr(y,s1)

str0 = ['噪声方差：' num2str(kd),'与s0相关结果']

str1 = ['噪声方差：' num2str(kd),'与s1相关结果']

figure

subplot(211)

stem(Nt2,z0)

title(str0)

subplot(212)

stem(Nt2,z1)

title(str1)

end

你可以看看这个，你怎么和之前别人提的问题一样啊

---