在數值方法中，需要Matlab來做這道題（全程需要，急用，需要master！！ ）

1.建立乙個新的 m 檔案並儲存被積函式：

function f=fun(x)

f=1./x;

2.建立另乙個 M 檔案並將其放入 Romberg 整合程式中：

function [i,t]=romberg(fun,a,b,eps,varargin)

Romberg 公式求解數值積分。

輸入引數：-fun： integrand。

a、b：積分間隔的端點。

EPS：精度要求，預設值為1E-6

p1,p2,..樂趣。

輸出引數：-i：得到的積分值。

t：Romberg積分過程得到的下三角矩陣。

if nargin<4|isempty(eps);

eps=1e-6;

endn=1;h=b-a;

t(1,1)=h/2*sum(feval(fun,[a,b],varargin));

tol=1;

while tol>eps

h=h/2;n=2*n;k=log2(n);

x=a+(2*(1:n/2)-1)*h;

fx = feval(fun,x,varargin);計算函式的值。

t(k+1,1)=1/2*t(k,1)+h*sum(fx);

for j=1:k

t(k+1,j+1)=(4^j*t(k+1,j)-t(k,j))/(4^j-1);

endtol=abs(t(k+1,k+1)-t(k,k));

endi=t(k+1,k+1);

3.在命令視窗中呼叫計算：

y=vpa(romberg(@fun,1,3,1e-5),6)y =

尋找異構陣列的問題...

由於要計算剪下頻率和相位角裕量，因此可以考慮直接使用裕量函式繪製波特圖

s=tf('s');

g=10/(s*(2*s+1)*(s^2+;

margin(g)

從圖中可以看出，剪下頻率 wc= 和相角裕度為 -136 度。 顯然，閉環系統是不穩定的。

由於閉環系統的不穩定性，時域指標（過衝和調整時間）不可用，如果要用階躍響應曲線進行驗證，可以使用：

step(feedback(g,1))

不難看出，閉環系統的階躍響應確實是發散的。

呃：可惜我不明白你的專業，MATLAB還可以，你能詳細描述一下問題嗎？

1.求解單個引數的解析形式。

a=dsolve('dp1=40*p1/v','dp2=40*p1/v-40*p2/v')

答案是：>ans =

2*c6*exp((40*t)/v)

c5/exp((40*t)/v) +c6*exp((40*t)/v)

引入源 P1（0）=10C，P2（0）=20C，求解 C5=5C，C6=15C

2.取濃度資料中的兩個點，p2（20）=，p2（40）=，把它們帶入p2的解析公式中，用fsolve（）求解c和v，希望能幫到你。

題主的問題屬於數理統計繪圖問題。 其實施過程。

1. 使用 unifrnd 隨機函式並用整數函式四捨五入，生成 50 100 範圍內的等級向量。

2. 使用平均函式計算其平均分數。

3、使用for迴圈語句對分數大於等於60分的資料進行排序，然後找出4個，計算及格率。

5. 使用 for 迴圈語句，然後使用 if 條件語句來區分“90 100 點”。','成績 80 89 分','成績 70 79 分','成績 60 69 積分','得分在60分以下“，並打分。

6. 使用 sum 函式計算條件數。

7. 使用條形函式製作直方圖。

8.使用餅圖功能製作餅圖。

根據以上要求，可以得到以下結果。

:pi/20:pi;

y1=sin(x)

y2=-sin(x)

y3=sin(x)*sin(9*x)

figure(1)

subplot(3,1,1)

plot(x,y1)

axis([0,pi,-1,1])

subplot(3,1,2)

plot(x,y2)

axis([0,pi,-1,1])

subplot(3,1,3)

plot(x,y3)

axis([0,pi,-1,1])

title('連續調製波形')

9y1=sin(t)

y2=-sin(t)

y3=sin(t)*sin(9*t)

7.函式檔案。

function [s,c]=yuan(r)s=r.^2*pi

c=2*pi*r

留下郵箱！！ 不給郵箱怎麼給源程式。。。