大家好,今天小编关注到一个比较有意思的话题,就是关于gridsum的问题,于是小编就整理了2个相关介绍gridsum的解答,让我们一起看看吧。
余弦函数的功率谱密度?
模型是给出了一个加速度信号,xn2=t*exp(a*t)*(sin(pi*t)+b*sin(8*pi*t)),要通过两次积分求出位移信号,其中a,b为常数,我取的1。其中加速度采样频率是100Hz,位移采样频率是5Hz,并在位移信号上添加10%RMS的高斯白噪声,要求用welch方法计算位移的功率谱密度。
我的程序:
xn1(1)=0; %速度
xn(1)=0; %位移
Fsa=100; %加速度采样频率
Fs=5; %位移采样频率
t=0;
for n=1:1:10000 %进行加速度采样
xn2(n)=t*exp(t)*(sin(pi*t)+sin(8*pi*t));
t=t+1/Fsa;
end
for n=1:1:500 %进行位移采样
xn1(n+1)=xn1(n)+xn2(n*20)/Fs;
xn(n+1)=xn(n)+xn1(n)/Fs;
n=n+1;
end
rmsx=sum(xn.*xn); %加噪声
rmsx=(rmsx/length(xn))^0.5;
nxn=xn+0.1*rmsx*randn(1,length(xn));
Nfft=512;
window=hanning(256);
noverlap=128;
dflag='none';
[Pxx,f]=pwelch(nxn,window,noverlap,Nfft,Fs);
plot(f,10*log10(Pxx));
xlabel('频率/Hz');ylabel('功率谱/dB');
title('PSD--Welch方法');grid on;
在三维空间直角坐标系中的五点任取几个点,能确定几个平面?
%P1=[x1,y1,z1];
%P2=[x2,y2,z2];
%P3=[x3,y3,z4];
%P=[P1;P2;P3];
% | x1 y1 z1 |
%P =| x2 y2 z2 |
% | x3 y3 z3 |
P=rand(3);%三点的xyz坐标正好是合拼成3x3的矩阵,这里用3x3的随机矩阵测试
%利用随机数生成3个点的坐标
%实际使用中,将已知的三点坐标按上述顺序排成P矩阵既可以了
A=diff(P);
n=cross(A(1,:),A(2,:)); %n是三点组成的平面的法向量
if all(n==0),error('not a triangle');end %如果n是0向量,表示三点共线
A=[A;n]; %从三点坐标很容易求得三点平面所在的方程
B=sum(diff(P.^2),2)/2; %而外接圆心都各点距离相等,又可以列几条方程
B=[B;dot(n,P(1,:))]; %将方程联立可以解得圆心坐标cc
cc=(A\B)'; %这里有个复杂的推导过程,
%最终结果cc(1) cc(2) cc(3)就是圆心的x,y,z坐标
r=sqrt(sum((cc-P(1,:)).^2)); %点cc到任意一点的距离都可以求出圆的半径;
C1=P(1,:)-cc; %第一点到圆心的向量C1
C2=cross(C1,n./sqrt(sum(n.^2))); %C2是三点所在平面内与C1垂直的向量
theta=(0:360)';%画圆所需的角度采样间隔,0到360度,间隔可以自己根据需要调整
cirxyz=repmat(cc,[length(theta) 1])+cosd(theta)*C1+sind(theta)*C2;
%cirxyz是361x3的矩阵,3列分别是圆上每点的x,y,z坐标
plot3(P(:,1),P(:,2),P(:,3),'ko');grid on;hold on; %画出三点
plot3(cirxyz(:,1),cirxyz(:,2),cirxyz(:,3),'r-');hold off; %画出圆
以下是某次空间随机三点和他们的外接圆图像
到此,以上就是小编对于gridsum的问题就介绍到这了,希望介绍关于gridsum的2点解答对大家有用。