【文章內容簡介】
如果操作者希望停止這種判斷,則可以退出程序。③如果所輸入的自然數是一個合數,除了給出其不是素數的結論外,還應給出至少一種其因數分解形式。例:輸入 6, 因為6不是素數。則程序中除了有“It is not one prime”的結論外,還應有:“6=2*3”的說明。 【程序】:a=input(39。請輸入一個整數(end in 0):39。)。while(a~=0) if a==1 disp(39。1 is not a prime or a posite number.39。) elseif isprime(a)==1 fprintf(39。%d is a prime.\n39。,a) elseif isprime(a)==0 f=factor(a)。m=length(f)。 fprintf(39。%d is not a prime,%d=%d39。,a,a,f(1)) for i=2:m fprintf(39。*%d39。,f(i))。 end fprintf(39。\n39。)。 end a=input(39。請輸入一個整數(end in 0):39。)。end【運行結果】:請輸入一個整數(end in 0):11 is not a prime or a posite number.請輸入一個整數(end in 0):22 is a prime.請輸入一個整數(end in 0):3535 is not a prime,35=5*7請輸入一個整數(end in 0):0實驗四 MATLAB的符號計算與SIMULINK的使用實驗任務1.掌握MATLAB符號計算的特點和常用基本命令;2.掌握SIMULINK的使用?!境绦?】syms a b c dA=[a b。c d]。B=[1,4。2,9]。 det(A),det(B),eig(A),eig(B)【運行結果1】ans = a*db*cans = 1ans = 1/2*a+1/2*d+1/2*(a^22*a*d+d^2+4*b*c)^(1/2) 1/2*a+1/2*d1/2*(a^22*a*d+d^2+4*b*c)^(1/2)ans = 【程序2】r=solve(39。x^3+5*x^29*x+1739。)r4=vpa(r,4)r10=vpa(r,10)【運行結果2】r =1/3*(557+3*18849^(1/2))^(1/3)52/3/(557+3*18849^(1/2))^(1/3)5/3 1/6*(557+3*18849^(1/2))^(1/3)+26/3/(557+3*18849^(1/2))^(1/3)5/3+1/2*i*3^(1/2)*(1/3*(557+3*18849^(1/2))^(1/3)+52/3/(557+3*18849^(1/2))^(1/3)) 1/6*(557+3*18849^(1/2))^(1/3)+26/3/(557+3*18849^(1/2))^(1/3)5/31/2*i*3^(1/2)*(1/3*(557+3*18849^(1/2))^(1/3)+52/3/(557+3*18849^(1/2))^(1/3)) r4 = .*i .858+*i r10 = .*i .858375306+*i【simulink仿真】【仿真結果】實驗五 MATLAB在控制系統(tǒng)分析中的應用一、實驗任務1.掌握MATLAB在控制系統(tǒng)時間響應分析中的應用;2.掌握MATLAB在系統(tǒng)根軌跡分析中的應用; 3. 掌握MATLAB控制系統(tǒng)頻率分析中的應用; 4. 掌握MATLAB在控制系統(tǒng)穩(wěn)定性分析中的應用1. 求下面系統(tǒng)的單位階躍響應【程序】num=[4] 。 den=[1 , 1 , 4] 。step(num , den)[y , x , t]=step(num , den) 。tp=spline(y , t , max(y)) %計算峰值時間ymaxmax(y) %計算峰值【結果】tp =ymax =2. 求下面系統(tǒng)的單位脈沖響應:【程序】num=[4] 。 den=[1 , 1 ,4] 。impulse(num,den)【結果如右圖】:求系統(tǒng)的零輸入響應和脈沖響應。【程序如下】: 【結果如下圖】a=[0 , 1 。 10 , 2] 。b=[0 。 1] 。c=[1 , 0] 。 d=[0] 。x0=[1 ,0]。subplot(1 , 2 , 1) 。 initial(a , b , c ,d,x0)subplot(1 , 2 , 2) 。 impulse(a , b , c , d):輸入正弦信號時,觀察輸出信號的相位差?!境绦蛉缦隆浚? 【結果如下圖】:num=[1] 。 den=[1 , 1] 。t=0 : : 10 。u=sin(2*t) 。hold onplot(t , u , ‘r’)lsim(num , den , u , t),求出周期為4秒的方波的輸出響應 【程序如下】: 【結果如下圖】:num=[2 5 1]。den=[1 2 3]。t=(0:.1:10)。period=4。u=(rem(t,period)=period./2)。lsim(num,den,u,t)。,繪制系統(tǒng)的根軌跡,并分析其穩(wěn)定性 【程序如下】: