1樓:匿名使用者
% 個人基本資訊,請自行修改
xuehao = '2009301234';
addr = '北京市海淀區xx大學xx樓***房間';
name = '張三丰';
% 第一題:要求影象橫座標為你的學號,縱座標為你的地址,標題為你的姓名
% (1)對訊號y=sin(314t)+sin(628t)加上白噪聲,畫出影象
fs = 1000; % 取樣率
t = 0:1/fs:10-1/fs; % 時間軸
y = sin(314*t)+sin(628*t); % 訊號
snr = 5; % 訊雜比
% 後面3條語句可以合併為1條:y1 = awgn(y, snr),但需要通訊工具箱的支援
np = 10^(-snr/10); % 噪聲功率
noise = sqrt(np) * randn(size(y)); % 產生白噪聲
y1 = y + noise; % 以訊雜比新增白噪聲
% 繪製訊號波形
figure(1), plot(t,y1);
xlabel(xuehao);
ylabel(addr);
title(name);
% (2)進行離散傅立葉變換,並畫出傅立葉變換後的頻譜圖、功率譜圖
n = length(t); % 訊號取樣點數
freq = fs/n*(-n/2:n/2-1); % 建立頻率軸
y1_spec = fftshift(fft(y1)); % 訊號的頻譜
% 繪製訊號頻譜
figure(2), plot(freq, abs(y1_spec));
xlabel(xuehao);
ylabel(addr);
title(name);
% 繪製訊號功率譜
p = abs(y1_spec).^2; % 訊號的功率譜
figure(3), plot(freq, p);
xlabel(xuehao);
ylabel(addr);
title(name);
% 第二題:用三次多項式擬合,求出表示式的值,畫出原數值圖與擬合示意圖,
% 圖的標題設定為你的學號,並測算x為7.5、12時的y值
x = 1:9;
y = [10.2 14.1 13.9 12.3 10.5 11 12.4 13.5 16];
% 用三次多項式擬合,從效果看差強人意(對本題資料而言,以4-6階為宜)
n = 3; % 多項式階次
p = polyfit(x, y, n); % 多項式擬合
figure(4)
plot(x, y, 'bo-', x, polyval(p,x), 'r:x');
legend('實驗資料', '擬合資料');
title(xuehao);
% 輸出多項式及測算點資訊
poly = vpa(poly2sym(p,'x'), 3);
fprintf('\n擬合多項式為:\n y = %s', char(poly))
y2 = polyval(p, [7.5 12]);
fprintf('\nx為7.5、12時的y測算值分別為%.2f、%.2f\n', y2)
2樓:匿名使用者
第一問subplot(2,1,1)
ezplot('sin(2*x)')
title('你的單位')
xlabel('姓名');
ylabel('學號')
subplot(2,1,2)
ezplot('sin(x)*cos(2*x)')title('你的單位')
xlabel('姓名');
ylabel('學號')
第二問x=1:9;
y=[10.2 14.1 13.9 12.3 10.5 11 12.4 13.5 16];
p=polyfit(x,y,3);
plot(x,y,'r')
y1=polyval(p,x);
hold on
plot(x,y1,'y')
y2=polyval(p,[7.5 12]) %7.5,12時的y值
matlab求兩個曲線的交線
x,y,z sphere 50 mesh 5 x,5 y,5 z 畫球面hold on ezmesh 0 x y 5 5 畫平面 x,y,z meshgrid linspace 5,5 contourslice x,y,z,x y z,5 x,5 y,5 z,0 0 畫交線 axis equal v...
如何用matlab求兩個函式的卷積
雙槍老椰子 function convolution t 4 0.001 4 e rectpuls t 0.25,1.5 plot t,e pause h tripuls t 1,2,1 plot t,h pause s 8 0.001 8 f conv e,h plot s,f 卷積後訊號就變長了...
目標函式有兩個未知變數,怎麼用matlab求最優解
z x 1.8 x 1 400 x 2 2 1.8 x 1 2 146.24 x 1 1270 50 x 2 3.24 x 1 3 262.872 x 1 2 3771 x 1 82375 100 x fmincon x z x 0 0 0 0 得到的結果為 x 35.6887 195.5207 其...