
clear all;
g=[1 0 1 1 1 0 0 1];%基带信号
f=100; %载波频率
t=0:2pi/99:2pi;
cp=[];sp=[];
mod=[];mod1=[];bit=[];
for n=1:length(g);
if g(n)==0;
die=-ones(1,100); %Modulante
se=zeros(1,100); %
else g(n)==1;
die=ones(1,100); %Modulante
se=ones(1,100); %
end
c=sin(ft);
cp=[cp die];
mod=[mod c];
bit=[bit se];
end
bpsk=cpmod;
subplot(2,1,1);plot(bit,'LineWidth',15);grid on;
title('Binary Signal');
axis([0 100length(g) -25 25]);
subplot(2,1,2);plot(bpsk,'LineWidth',15);grid on;
title('ASK modulation');
axis([0 100length(g) -25 25]);
第一种方法:利用dsolve函数求微分方程的符号解(通解):对于一些不是很难,要求出通解的微分方程,用dsolve函数求解。
1、 打开Matlab软件-->点击新建脚本菜单,新建一个脚本文件用于编写微分方程求解程序。
2、 输入微分方程求解程序-->点击保存-->点击运行。
3、在matlab的命令窗口即可看到求解结果,是一个关于参数a,b的表达式
第二种方法:利用Matlab中的solver函数(包括ode45、ode23、ode15s等)来求解微分方程的数值解,这种方法是最常用的方法,对于dsolve函数难以求解的方程就可以利用这种方法求解方程的数值解。
1、 打开Matlab-->新建一个脚本文件用于编写求解程序。
2、 在脚本文件中输入求解程序-->点击保存-->点击运行。
3、 这里需要先编写一个函数文件用于表示方程-->点击保存-->编写求解程序-->点击保存-->点击运行。
4、在figure页面就可以看到所求解的微分方程的图形。
看您的采样周期是1/5,直接采样就好了(不要纠结于冲击函数的数学表达在matlab里怎么对应,这是硬要用连续函数表示离散信号的结果,matlab里只能表示离散信号,所以直接忽略之即可)
t
=
[-10:1/5:10];
k
=
2;
x
=
exp(1ipikt^2);
得到x是一个向量,是t中各个采样时刻对应的采样值。
以上就是关于求MATLAB仿真程序全部的内容,包括:求MATLAB仿真程序、如何利用matlab解微分方程、想知道一个信号被采样后如何用matlab程序写出来!求程序!等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
微信扫一扫
支付宝扫一扫
评论列表(0条)