Question

麻烦各位大神帮忙解释一下，matlab中seqmult的用法

我知道这个意思是求序列积，但是为什么不能让两个信号求积？或者是两个信号求积应该怎么编写matlab程序？例如：x1=sint1, x2=cost2, t1=1:1/fs1:1, t2=1:1/fs2:1, fs1=8000,fs2=10000. 如何作图：y=x1*x2的时频图？？

Answer 1

>> help seqmult

未找到 seqmult。

请使用帮助浏览器的搜索字段搜索文档，或者
键入 "help help" 获取有关帮助命令选项的信息，例如方法的帮助。
说明：这不是matlab的函数，可能是自写的。但是，两个信号的乘积是可以用点乘计算的，如：
t=0:1/1000:1;
x1=cos(t);x2=sin(t);y=x1.*x2;

追问

这个函数是从书中看到的，经你一说才发现没有这个语句，但是其实我的主要问题是怎么产生如下信号：s(t)=(sin(π*n/N))^2*p(t-n*T) , 其中p(t)=exp(j(2*k*f0*t+π*k*t^2)，其中t是信号采样：t=0:1/10000:1; N=46, n=N-1，T=0.03（脉冲时间），f0=4000（起始频率），k=B/T，B=500（带宽）。

追答

t=0:1/10000:1;
N=46; n=N-1;T=0.03;%（脉冲时间）
f0=4000;%（起始频率）
B=500;%带宽
k=B/T;
p=exp(1j*(2*k*f0*(t-n*T)+π*k*(t-n*T).^2));

s=(sin(π*n/N))^2*p;

按上面的试试

追问

试了，结果明显不对，只画了时域图，但是明显不对，再麻烦你帮我看看怎么改

close all

clear all

clc

t=0:1/10000:1;     

N=46; n=N-1;T=0.03;

f0=5000;

B=500;

k=B/T;

p=exp(1j*(2*k*f0*(t-n*T)+pi*k*(t-n*T).^2));

s=(sin(pi*n/N))^2*p;

figure(1);

plot(s);

追答

可这样：
close all
clear all
clc
t=0:1/10000:1;
N=46; n=N-1;T=0.03;
f0=5000;
B=500;
k=B/T;
p=exp(1j*(2*k*f0*(t-n*T)+pi*k*(t-n*T).^2));
s=(sin(pi*n/N))^2*p;
figure(1);
plot(t,s);