哪位大侠会写matlab求解分段函数FFT的程序 MATLAB分段函数
\u54ea\u4f4d\u5927\u4fa0\u80fd\u5e2e\u5fd9\u5199\u4e00\u4e2a\u6539\u8fdb\u7684FFT\u7b97\u6cd5\uff0c\u8981\u6c42\u6bd4FFT\u8ba1\u7b97\u901f\u5ea6\u5feb\uff0c\u5982\u679c\u80fd\u63d0\u4f9bMatlab\u7684\u4ee3\u7801\u5c31\u66f4\u4e0d\u80dc\u611f\u6fc0\uff01\uff01\uff01\u4f60\u662f\u4e0d\u662f\u8f93\u5165\u6570\u636e\u592a\u591a\u4e86\u3002\u6570\u636e\u8d8a\u5c11\uff0c\u8fd0\u7b97\u4e5f\u5c31\u8d8a\u5feb~
\u6b65\u9aa4\u5982\u4e0b
1\u3001\u6253\u5f00MATLAB\u8f6f\u4ef6\uff0c\u5982\u56fe\u6240\u793a\u3002
2\u3001\u5efa\u7acb\u4e00\u4e2a\u811a\u672c\u6587\u4ef6\uff0c\u5177\u4f53\u65b9\u6cd5\u5982\u56fe\u6240\u793a\u3002
3\u3001\u5b9a\u4e49\u53d8\u91cf\u3002
4\u3001\u5efa\u7acb\u5faa\u73af\uff0c\u6c42\u89e3\u5206\u6bb5\u51fd\u6570\u3002
5\u3001\u91c7\u7528\u4ee5\u4e0b\u6307\u4ee4\u753b\u56fe\u3002
6\u3001\u753b\u51fa\u7684\u56fe\u7247\u5982\u56fe\u6240\u793a\u3002
\u6269\u5c55\u8d44\u6599
20\u4e16\u7eaa70\u5e74\u4ee3\uff0c\u7f8e\u56fd\u65b0\u58a8\u897f\u54e5\u5927\u5b66\u8ba1\u7b97\u673a\u79d1\u5b66\u7cfb\u4e3b\u4efbCleve Moler\u4e3a\u4e86\u51cf\u8f7b\u5b66\u751f\u7f16\u7a0b\u7684\u8d1f\u62c5\uff0c\u7528FORTRAN\u7f16\u5199\u4e86\u6700\u65e9\u7684MATLAB\u30021984\u5e74\u7531Little\u3001Moler\u3001Steve Bangert\u5408\u4f5c\u6210\u7acb\u4e86\u7684MathWorks\u516c\u53f8\u6b63\u5f0f\u628aMATLAB\u63a8\u5411\u5e02\u573a\u3002\u523020\u4e16\u7eaa90\u5e74\u4ee3\uff0cMATLAB\u5df2\u6210\u4e3a\u56fd\u9645\u63a7\u5236\u754c\u7684\u6807\u51c6\u8ba1\u7b97\u8f6f\u4ef6\u3002
MATLAB\uff1a\u7edf\u4e00\u4e86\u7528\u4e8e\u4e00\u7ef4\u3001\u4e8c\u7ef4\u4e0e\u4e09\u7ef4\u6570\u503c\u79ef\u5206\u7684\u51fd\u6570\u5e76\u63d0\u5347\u4e86\u57fa\u672c\u6570\u5b66\u548c\u5185\u63d2\u51fd\u6570\u7684\u6027\u80fd
MATLAB Compiler\uff1a\u53ef\u4ee5\u4e0b\u8f7d MATLAB Compiler Runtime (MCR)\uff0c\u7b80\u5316\u7f16\u8bd1\u540e\u7684\u7a0b\u5e8f\u548c\u7ec4\u4ef6\u7684\u5206\u53d1
Image Processing Toolbox\uff1a\u901a\u8fc7\u4eae\u5ea6\u6307\u6807\u4f18\u5316\u8fdb\u884c\u81ea\u52a8\u56fe\u50cf\u914d\u51c6
Statistics Toolbox\uff1a\u589e\u5f3a\u4e86\u4f7f\u7528\u7ebf\u6027\u3001\u5e7f\u4e49\u7ebf\u6027\u548c\u975e\u7ebf\u6027\u56de\u5f52\u8fdb\u884c\u62df\u5408\u3001\u9884\u6d4b\u548c\u7ed8\u56fe\u7684\u754c\u9762
System Identification Toolbox\uff1a\u8bc6\u522b\u8fde\u7eed\u65f6\u95f4\u4f20\u9012\u51fd\u6570
\u53c2\u8003\u8d44\u6599\uff1a\u767e\u5ea6\u767e\u79d1 MATLAB\u53d1\u5c55\u5386\u7a0b
最大的频率在f=1.15,w=2*pi*f=7.22处
%--------------------------------
Ts=0.2;
fs=1/Ts;
t1=1:Ts:10;
y1=sin(2*pi*t1) %(0<t<10)
t2=10:Ts:20;
y2=1.1*sin(2.1*pi*t2)% (10<t<20)
t3=20:Ts:30;
y3=1.2*sin(2.2*pi*t3) %(20<t<30)
t4=30:Ts:40;
y4=5*sin(2.3*pi*t4) %(30<t<40)
t5=40:Ts:50;
y5=1.5*sin(2.4*pi*t5) %(40<t<50)
t6=50:Ts:60;
y6=1.2*sin(2.5*pi*t6) %(50<t<60)
y=[y1 y2 y3 y4 y5 y6];
figure;plot(y);
yLength=length(y);
Y=fft(y);
Pyy = Y.* conj(Y) / yLength;
haflength=floor(yLength/2);
f=0:fs/2/haflength:fs/2;
bar(f,Pyy(1:haflength+1));
xlabel('Hz');
grid on;
绛旓細--- 鏈澶х殑棰戠巼鍦╢=1.15,w=2*pi*f=7.22澶 --- Ts=0.2;fs=1/Ts;t1=1:Ts:10;y1=sin(2*pi*t1) %(0<t<10)t2=10:Ts:20;y2=1.1*sin(2.1*pi*t2)% (10<t<20)t3=20:Ts:30;y3=1.2*sin(2.2*pi*t3) %(20<t<30)t4=30:Ts:40;y4=5*sin(2.3*pi*t4) %...
绛旓細function eq=nuowei(x,z) % 淇濆瓨涓 nuowei.m 鏂囦欢 Vwz=-62700*z.*x./(sqrt(400^2-x.^2).*(416.5-sqrt(400^2-x.^2)).^2);Vwx=62700./(416.5-sqrt(400^2-x.^2));Vbx=91200./(424-sqrt(480^2-x.^2));eq=400*sqrt((Vwx-Vbx).^2+Vwz.^2);--- 杩愯浠ヤ笅锛>...
绛旓細integfunction=A;%绉垎鍑芥暟涓婄殑涓浜涚偣 time=0:timeg:n*timeg;%绉垎鍖洪棿 Ssp(1)=0; %绉垎鍊肩殑鍒濆 for i=1:(n-1)xx=[time(i),time(i+1)]; %渚濇閫夊彇涓や釜鐐 ff=[integfunction(i),integfunction(i+1)];%閫夊彇鐐规墍瀵瑰簲鐨勭Н鍒嗗嚱鏁板 pp=spline(xx,ff);%涓夋鏍锋潯宸 int_pp=...
绛旓細y=(-log10((0.1*20-0.2*x)./(20+x))).*(0<x&x<10)+7.*(x==10)+(14+log10((0.2*x-0.1*20)./(20+x))).*(10<x&x<20);plot(x,y,10,7,'o')
绛旓細Q = atan2(y, x);% 瑙涓Q11 = Q - a1 Q21 = a1 + a2% 瑙ｄ簩Q12 = Q + a1Q22 = -a1 - a2% 缁樺浘clfset(gcf, 'color', 'w')set(gca, 'vis', 'off')rectangle('Pos', [-20-x0 y0-30 40 22.5], 'LineStyle', '--', 'Linewidth', 4, 'EdgeColor', [1 1 1]*0...
绛旓細s0=1;[T,X]=ode45('lalan',[0 2],[i0;s0]);plot(T,X(:,1),T,X(:,2))--- function dy=lalan(t,x) %淇濆瓨涓簂alan.m鏂囦欢 lam=3;mu=4;dy=[lam*x(2)*x(1)-mu*x(1);-lam*x(2)*x(1)];
绛旓細{ datatype temp;linklist p;p=head->next;temp=(p->data);head->next=p->next;free(p);return (temp);}
绛旓細-0.2097 (-0.508, 0.08866)p6 = 73.6 (-32.86, 180.1)p7 = -1.067e+004 (-2.636e+004, 5025)Goodness of fit:SSE: 0.02709 R-square: 0.9978 Adjusted R-square: 0.9913 RMSE: 0.1164 涔熷彲浠ュ弬鑰冩枃鐚細鏈灏忎簩涔樻洸绾挎嫙鍚堝強Matlab瀹炵幇 涓婂鏍℃暟鎹簱涓嬪幓 ...
绛旓細鍑芥暟鍚嶄负f2 搴曚笅灏变笉鑳戒娇鐢╢2浜 浣犳妸绗竴琛屽簳涓嬬殑f2鍏ㄩ儴鎹㈡垚y灏辫浜
绛旓細3棰 >>x=0:0.1:8*pi;>> y=5*sin(x./2);>> plot(x,y)>> grid on >> xlim('auto')>> ylim('auto')>> title('y=5sin(x/2)')>> xlabel('xlabel')>> ylabel('ylabel')缁撴灉濡傚浘
