应用matlab程序求函数值~~

x0=0

x=x0

x2=[]

for i=1:n

x2=[x2,x(i).^2]

x=[x,x(i)^2+0.25]

end

plot(0:n-1,x2)

xlabel('n')

ylabel('x_n^2')

%%%%%%while

n=200

x0=0

x=x0

x2=[]

i=1

while i<=n

x2=[x2,x(i).^2]

x=[x,x(i)^2+0.25]

i=i+1

end

plot(0:n-1,x2)

xlabel('n')

ylabel('x_n^2')

我画的是n与xn^2的曲线毕唤,你如果想要n与xn的曲线,

你用下面稿盯三行

plot(0:n-1,x(1:end-1))

xlabel('n')

ylabel('x_n')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%比方你想要0,5,10,20,50,100,200的值(下次别用省略,谁知道你想要哪些啊,我这手敬凯里就随便给了一些数)

nfind=[0,5,10,20,50,100,200]

xfind=x(nfind+1)

（1-1）：

function [b,y] = by(m,n)

b=lcm(m,n)%最小公倍数

y=gcd(m,n)%最大公约数

end

（1-2）：

function t=sushu(n)

m=floor(sqrt(n))

for i=2:m

if mod(n,i)==0

t=0

break

else

t=1

end

end

end

（派正猛2）：可以一次算出每一项，发现他们都是同一个值，这个值的求法：SUM=sum(1:n^2)/n%n阶魔幻阵特有的和

function [A,H,L,DJ,SUM]=mag(n)

A=magic(n)%产生一个n阶魔幻阵

H=sum(A)%每行的和

L=sum(A)%每列的和

DJ=sum(diag(A))%对角线和

SUM=sum(1:n^2)/n%n阶魔幻阵特有的和

end

（3）：这题我给出“该数各位数字之和为奇数”清漏的求法，结合上面有判断素数的代码自己可以求出

该数各位数字之和为奇数：

function funct

for i=2:999

%l=floor(log10(i))+1 %提取位数

n1=floor(i/100) %取出百位数字n1

n2=mod(floor(i/10),10)%取出十尘桥位数字n2

n3=mod(i,10) %取出个位数字n3

nsum=n1+n2+n3

if rem(nsum,2)~=0%判断是否为奇数

disp (i)

end

end

end

自己可以调试看看，祝你成功！

程序：

wlp = 0.35*pi

whp = 0.65*pi

wc = [wlp/pi,whp/pi]

N = ceil(8/0.15)

n=0:N-1 

window= hanning(N)

[h1,w]=freqz(window,1)

subplot(411) 

stem(window)%画窗函数

axis([0 60 0 1.2])

grid

xlabel('n') 

title('Hanning窗函数')

subplot(412)   

plot(w/pi,20*log(abs(h1)/abs(h1(1))))  %画Hanning窗频谱图

axis([0 1 -350 0])

grid

xlabel('w/pi')

ylabel('幅度(dB)'猜配) 姿芹

title('Hanning窗函数的频谱')

hn = fir1(N-1,wc, hanning (N))%基于窗函数的 FIR 滤波器

[h2,w]=freqz(hn,1,512)   %freqz是迹兆毕用离散傅里叶变换的标准公示计算的 fft使用快速傅里

subplot(413)  

stem(n,hn) 

axis([0 60 -0.25 0.25])

grid

xlabel('n')

ylabel('h(n)') 

title('Hanning窗函数的单位脉冲响应')

subplot(414) 

plot(w/pi,20*log(abs(h2)/abs(h2(1))))

grid

xlabel('w/pi')

ylabel('幅度(dB)')

title('滤波器相对幅度响应')

结果：