【正文】 2. 一物理系統(tǒng)可用下列方程組來表示：從鍵盤輸入mm2和θ的值，求aaN1和N2的值。輸入m2=39。B=[0。解：M文件：function [p] = prime(p)% 輸入p的范圍，找出其中的素數(shù)m=p(length(p))。 %找出10到99內(nèi)的所有素數(shù)p=10*rem(p,10)+(prem(p,10))/10。)。命令文件：clc。n3=39。n1=input(39。)。解：M文件如下：clc。y1=x.^2。k39。b39。)。title(39。subplot(2,2,1)。subplot(2,2,2)。bar(x,y2,39。k39。y=(x+sqrt(pi))/(exp(2)).*(x=0)+*log(x+sqrt(1+x.^2)).*(x0)。b=input(39。m39。x=linspace(5,5,21)。)。s=0:pi/100:pi/2。title(39。subplot(2,2,3)。shading interp。,[1 0 0],... 39。)。,[,],39。\leftarrow x^2exp(2x)39。v=10*exp(*x).*sin(2000*pi**x+pi)。,39。yaxis39。pause %按任意鍵繼續(xù)set(hs,39。運行結(jié)果：按任意鍵繼續(xù)：4. 以任意位置子圖形式繪制出正弦、余弦、正切和余切函數(shù)曲線。Xaxis39。axis([5,5,5,5,0,1])grid off。style39。axis off。x=rand(1,30000)。t=45+50*rand(100,5)。T1=interp1(h,t1,39。P=polyfit(x,y,5)y1=polyval(P,x)。(3) 當x取矩陣A的每一元素時，求P(x)的值。p2=[0,0,0,p2]。y=polyval(p,A) %x取矩陣A的每一元素時的p(x)值 運行結(jié)果：p = 0 0 0 0 1 3 8 7 11x = + + y = +003 * 實驗九 數(shù)值微積分與方程數(shù)值求解一、實驗目的1. 掌握求數(shù)值導數(shù)和數(shù)值積分的方法。f=inline(39。 %，可再縮小步長提高精度endg。sqrt(cos(t.^2)+4*sin(2*t).^2+1)39。A=[6 5 2 5。y=[ ]。 x=A\b。 endelse %齊次方程組 disp(39。r39。[x,y]=line_solution(A,b) 運行結(jié)果：有無窮個解，特解x，基礎(chǔ)解系yWarning: Rank deficient, rank = 2, tol = . In line_solution at 11x = 2/11 10/11 0 0 y = 1/11 9/11 5/11 1/11 1 0 0 1 所以原方程組的通解是：，其中為任意常數(shù)。f39。myfun39。(2) 在[0,0]附近的最小值點和最小值。)。x0=。y(1)*y(3)。,[t0,tf],[0,1,1])plot(x,y) 運行結(jié)果：x = 0 y = 0 圖形：實驗十 符號計算基礎(chǔ)與符號微積分一、實驗目的1. 掌握定義符號對象的方法。639。syms x y。syms beta1 beta2 x。 (4) B的行列式值。0 1 0。syms x t a y z。 %(4)Ax1=diff(A,x,1)At2=diff(A,t,2)Axt=diff(Ax1,t)f=(x^22*x)*exp(x^2z^2x*z)。eval(dfxz) %符號運算返回數(shù)值運行結(jié)果：ans = 1/2 ans = Inf y1 = (2*sin(2*x))/x + (cos(2*x) 1)/x^2 y2 = (4*cos(2*x))/x (4*sin(2*x))/x^2 (2*(cos(2*x) 1))/x^3 Ax1 = [ a^x*log(a), 0][ t*sin(x), 1/x] At2 = [ 0, 6*t][ 0, 0] Axt = [ 0, 0][ sin(x), 0] Zx = (exp(x^2 + x*z + z^2)*((2*x 2)/exp(x^2 + x*z + z^2) + ((2*x x^2)*(2*x + z))/exp(x^2 + x*z + z^2)))/((2*x x^2)*(x + 2*z)) ans = 4/exp(1) 6. 用符號方法求下列積分。 %(2)f3=(x^2+1)/(x^4+1)。(2) 求級數(shù)的和函數(shù)，并求之和。)。taylor(f,x,6,1)運行結(jié)果：ans = x (x 1)^2/2 + (x 1)^3/3 (x 1)^4/4 + (x 1)^5/5 13. 求下列方程的符號解。3*x*exp(x)+5*sin(x)39。dsolve(39。) 運行結(jié)果：ans = (3*sin(5*x))/exp(2*x)5. 求微分方程組的通解。Dz=4*x4*y+2*z39。拼一個春夏秋冬！贏一個無悔人生！早安！—————獻給所有努力的人. 學習好幫手。t39。clc。,39。sqrt(x^2+y^2)10039。clc。解：M文件：clear all。clc。 %(4)F1=int(f1)F2=int(f2)F3=int(f3,0,inf)F4=int(f4,0,log(2)) 運行結(jié)果： F1 = (3^(1/2)*log(x^2 + 3^(1/2)*x + 1))/12 (3^(1/2)*(2*atan( (3^(1/2)*x^3)/3 (2*3^(1/2)*x)/3) 2*atan((3^(1/2)*x)/3)))/12 (3^(1/2)*log(x^2 3^(1/2)*x + 1))/12 F2 = 1/asin(x) F3 = (pi*2^(1/2))/2 F4 = exp(18729944304496077/9007199254740992)/3 + exp(6243314768165359/9007199254740992) + exp(6243314768165359/4503599627370496) 7/3實驗十一 級數(shù)與方程符號求解一、實驗目的1. 掌握級數(shù)求和的方法。clc。x=sym(39。 %(1)limit(f1)f2=(sqrt(pi)sqrt(acos(x)))/sqrt(x+1)。A=[a b c。clc。simplify(f1) %（1）問f2=(4*x^2+8*x+3)/(2*x+1)。513539。y=sym(39。3. 掌握求符號函數(shù)極限及導數(shù)的方法。 clc。[x,y]=ode45(39。g39。 y=u(2)。,optimset(39。y=X(2)。(1) 3x+sinxex=0在x0=。 endend clc。)。r39。有唯一解x39。3 4 2 2。I1=quad(f,0,2*pi) g=inline(39。dx=diff(g)/。1 2*x 3*x^2。3. 掌握常微分方程數(shù)值求解的方法。p4=conv(p2,p3)。其中A的值與第(3)題相同。:o39。)%室內(nèi)的3次樣條插值溫度T2=interp1(h,t2,39。 %生成100個學生5門功課成績[x,l]=max(P) %x為每門課最高分行向量,l為相應學生序號[y,k]=min(P)%y為每門課最低分行向列,k為相應學生序號mu=mean(P) %每門課的平均值行向量sig=std(P) %每門課的標準差行向量s=sum(P,2) %5門課總分的列向量[X,m]=max(s)%5門課總分的最高分X與相應學生序號m[Y,n]=min(s)%5門課總分的最低分Y與相應學生序號n[zcj,xsxh]=sort(s) %zcj為5門課總分從大到小排序，相應學生序號xsxh 運行結(jié)果：3. 某氣象觀測得某日6:00~18:00之間每隔2h的室內(nèi)外溫度（0C）如實驗表1所示。 %p=y/30000 %運行結(jié)果：mu = sig = p = 2. 將100個學生5門功課的成績存入矩陣P中，進行如下處理：(1) 分別求每門課的最高分、最低分及相應學生序號。2. 掌握數(shù)值插值與曲線擬合的方法及其應用。infinite39。Color39。Ylabel(39。 解：M文件：[x,y,z]=cylinder(3,500)。,39。zlabel(39。)。view39。fontsize39。,39。解：M文件：clc。,39。shading interp39。shading flat。)。[s,t]=meshgrid(s,t)。surfc(x,y,z)。[x,y]=meshgrid(x,y)。分析結(jié)果：由這8個圖知道， 當a,n固定時，圖形的形狀也就固定了，b只影響圖形的旋轉(zhuǎn)的角度。)。4. 繪制極坐標曲線ρ=asin(b+nθ)，并分析參數(shù)a、b、n對曲線形狀的影響。title(39。)。b39。b39。)。y2=cos(2x)39。title(39。x=pi:pi/100:pi。y3=y1.*y2。y=(+3*sin(x)./(1+x.^2))。y2=g(n2)。)。y1=f(n1)。n1=39。輸入矩陣x=39。解：：funct