用matlab程式設計解數學題(急)

時間 2021-10-14 22:41:45

1樓:天上的一條龍

檔案1jisi.m

%creat the function for jisi

%let x(1)=x, x(2)=y

function dx=jisi(t,x,a,b,c)

s=sqrt((c-x(1))^2+(a*t-x(2))^2);

dx=[b*(c-x(1))/s;b*(a*t-x(2))/s];

檔案2seajisi.m

clear,clf,shg

%set the definied time

%ts=0:0.05:.5;

ts=0:0.1:1.6;

n=length(ts);

x0=[0 0];

%a=20;b=40;c=15;

a=35;b=40;c=15;

opt=odeset('reltol',1e-6,'abstol',1e-9);

%[t,x]=ode45(@jisi,ts,x0,,a,b,c);

[t,x]=ode45(@jisi,ts,x0,opt,a,b,c);

%exact solution x1=c

y1=a*t;

%output t,x(t),y(t) and draw x(t),y(t)

[t,x,y1]

plot(t,x),grid,gtext('x(t)','fontsize',16),

gtext('y(t)','fontsize',16),pause

%draw y(x): the position of tatch jisi

plot(x(:,1),x(:,2),'r*'),grid

xlabel('x','fontsize',16),ylabel('y','fontsize',16)

2樓:丿happy灬殤丨

...。。。。//有關緝私艇用_p表示,有關走私船用_c表示,r為兩艇距離

//以緝私艇開始位置為座標原點

double x_p,y_p,x_c,y_c,t;

double dt,r,dx_p,dy_p,dx_c,dy_c,v_p,v_c;

dt=0.001;

v_p=40.0; v_c=20.0;

x_p=0.0; y_p=0.0;

x_c=15.0; y_c=0.0;

t=0.0;

dx_p=v_p; dy_p=0.0;

dx_c=0.0; dy_c=v_c;

r=sprt((x_c-x_p)*(x_c-x_p)+(y_c-y_p)*(y_c-y_p));

for(t=0.0;r<0.0001;t+=dt)cout<<"相遇時間 t="<

有人會用matlab 做數學題嗎?!!**等!急啊!請幫忙!!!

3樓:

十分明晰,趙本山啊,一定正確有``

4樓:

我郵箱[email protected]

樂意幫你解答

5樓:

要看是什麼題啊,不是所有的數學題都會做的

6樓:匿名使用者

會呢,是建模用嗎?我用過的嘍,還有lindo,lingo等等的

7樓:匿名使用者

涉及公式 程式設計類都行 但是 邏輯問題解決不了

求一個可以用matlab程式設計解決的物理或數學問題

8樓:

處理物理實驗資料不就是嗎。。。。我現在就在處理我苦逼的實驗資料。。。。

9樓:魂月夜

同學我是你的老師,你的ip地址已被擷取,自己動手做

用matlab解決一道數學問題

10樓:匿名使用者

v(t)=v-(v-v0)*exp(-t/tau)移項,後兩邊同取對數

log(v-v(t))=log(v-v0)-t/tau上面v,v(t),t已知,v0,tau未知。

令y=log(v-v(t))

x=t那麼就有kx+b=y

其中k=-1/tau,b=log(v-v0)上面的過程就是將非線性化成線性。

將每個資料代入會得到二元一次線性方程組

方程數6,變數數2,可以由最小二乘法求出。

以下是matlab**:

v=14;

t=[0.3,0.5,1.0,2.0,4.0,7.0];

vt=[5.6873,6.1434,7.1633,8.8626,11.0328,12.6962];

y=log(v-vt);

a=[t(:),ones(size(t(:)))];

kb=a\y(:);

tau=-1/kb(1)

v0=v-exp(kb(2))

%%下面是顯示擬合的結果

vt_fit=v-(v-v0)*exp(-t/tau)plot(vt,vt_fit,'*',vt,vt,'r')xlabel('vt原始資料')

ylabel('vt擬合資料')

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

填空:方法:化非線性問題成線性問題

v0=5.0001

tau=3.6165

11樓:匿名使用者

最簡單的方法,列方程求解。裡面只有兩個未知數,v0和時常數t7.1633=14-(14-v0)*exp(-1/t) (1)8.

8626=14-(14-v0)*exp(-2/t) (2)把(1)平方與(2)相比,就可以得到v0=4.9019

下面這道題,求數學建模問題的最優解,用matlab怎麼程式設計? 20

12樓:匿名使用者

function love%主函式,也可以把這一段放在命令窗中執行.那樣得把sub_f(x)函式存放在work中

clear;clc;

% x11,x12,x13,y22,y33,x22,x23 ->x(1),x(2),x(3),x(4),x(5),x(6),x(7)

maxf=inline('0.25*(x(1)+x(2)+x(3))*x(4)+(x(1)+x(2)+x(3))*(1-x(4))*0.2+(-x(2)+x(6)+x(7))*x(5)*0.

23+(-x(2)+x(6)+x(7))*(1-x(5))*0.185','x') ;

%st. x(2)*80+x(3)*100+x(7)*85<=900

%-x(3)-(7)<=0

%x(3)+x(7)<=3.5

a=[0 80 100 0 0 0 85

0 0 -1 0 0 0 -1

0 0 1 0 0 0 1];

b=[900;0;3.5];

%0<=x(2)<=8.2

%0<=x(4)<=1

%0<=x(5)<=1

lb=[-inf 0 -inf 0 0 -inf -inf];

ub=[inf 8.2 inf 1 1 inf inf];

x0=0.5*ones(1,7);

aeq=;beq=;

[x,f]=fmincon(maxf,x0,a,b,aeq,beq,lb,ub,@sub_f);

x=x,f=-f

function [f feq]=sub_f(x)

f(1)=(x(1)+x(2)+x(3))*x(4)*1.4+(x(6)-x(2)+x(7))*x(5)*1.65-7.5;

f(2)=(x(1)+x(2)+x(3))*(x(4)*1.4+6.1)+(x(6)-x(2)+x(7))*(x(5)*1.65+7.35)-76.5;

feq(1)=0;feq(2)=0;

結果:x =

1.0e+015 *

-4.4110 0.0000 -0.1376 0.0000 0.0000 -4.1338 0.1376

f =4.8952e+018

也可以看我的部落格文章的例子;

13樓:匿名使用者

用fmincon函式

你看看help fmincon

主要是建立兩個m檔案,一個myfun放你的f(x),注意,這裡要放-f(x),因為你要求最大值,fmincon是求最小值的

另一個mycon放非線性約束條件,這倆,(x11+x12+x13)y11*1.4+(x22-x12+x23)y22*1.65<=7.5

(x11+x12+x13)(y11*1.4+6.1)+(x22-x12+x23)(y22*1.65+7.35)<=76.5

把線性約束條件寫成矩陣a,b

直接呼叫

x=fmincon(@myfun,x0,a,b,,,lb,ub,@mycon)

lb,ub是自變數的上下限

x0為初值,一般需要多嘗試幾個初值

14樓:

其實數學建模還有很多軟體,這樣的優化問題都是用lingo求解的~

15樓:匿名使用者

對於優化問題,你可以嘗試使用lingo10,做起來比matlab要方便的多

求解,數學題,求解,數學題

72,減少了正方體的四個面,增加了兩個面,就是6 6 2 72平方釐米 其實只有從一角切是減少了72平方釐米 若從一邊切,就表面積不變,若從中間挖,表面積增加了72平方釐米 原表面積為 2x 9x8 9x6 8x6 348 平方釐米 最大的正方體的表面積為 6x6x6 216 平方釐米 減少的表面積...

初二數學題,急求解答,求解數學題。

洪山心雨 我想了會兒,應該還有更簡單的方法,我的方法如下 作ch垂直ao於h,令of x,oh 3,he ae 1,ch 1,co 根號下10 第一步 odf與 och相似,得 of oh od oc,帶入資料,可得od與x的關係 第二步 odf與 ceh相似,得 of ch od ce,其中,ce...

用算術方法解數學題

解 原甲 原乙 66 60 70 66 3 2 36 24後甲 後乙 66.25 60 70 66.25 25 13此時36 24 25 13 12 5 12 36 15 升 甲 15 3 2 10 升 乙 答 原來甲瓶藥水有15升,乙瓶藥水有12升.70 60 70 60 66 60 70 66 ...