推 fishooooo:抱歉～i是複數的i 03/05 20:49
> -------------------------------------------------------------------------- <
作者: zupo (幫MATLAB問題) 看板: MATLAB
標題: Re: [問題] 二次方程式（？）猜其中一值解另一值
時間: Thu Mar 6 04:07:35 2008
※ 引述《fishooooo (x8 x8 )》之銘言：
: 我有一個問題要求解 只有初步的想法 但是寫的過程卡卡的 想請各位大大給點意見
: 問題如下
: ｜-C11+C44*alpha^2-pho*v^2 (C13+C44)*alpha*i (e15+e31)*alpha*i |
: ｜ |
: ｜(C13+C44)*alpha*i C44+C33*alpha^2-pho*v^2 e15+e33*alpha^2 |=0
: ｜ |
: ｜(e15+e31)*alpha*i e15+e33*alpha^2 -(es11+es33*alpha^2) |
: ｜ |
: C11=226 e15=0.08 es11=0.19
: C13=121 e31=-2.1 es33=5.1
: C33=218 e33=6.9 pho=7.5
: C44=48
: 給定 v=2000:10:6000
: 計算每個給定v值的alpha值 如果給定的v值不能使det| |=0就換下一個v值求alpha
: 希望最後能畫出橫軸v,縱軸alpha的圖型
: 目前我在中間的處理有點不知道該如何寫
: 我是新手＠＠ 還請大家幫我解答一下 給我意見
: 感恩感恩<(_ _)>
不是很好解.
不過在超強力函數 fsolve 威能全開+先行宣告矩陣空間下
似乎得到不錯的解答
==============================================================================
function pttex64
global v
clc
format long g
v1 = (2000:10:6000)';
n = length(v1);
alpha = zeros(n,1);
for a=1:n
v = v1(a);
alpha(a) = fsolve(@nonldet,2);
end
anser = [v1 alpha]
function f = nonldet(alpha)
global v
i = (-1)^0.5;
C11=226;e15=0.08;es11=0.19;
C13=121;e31=-2.1;es33=5.1;
C33=218;e33=6.9;pho=7.5;C44=48;
A = [-C11+C44*alpha^2-pho*v^2,(C13+C44)*alpha*i,(e15+e31)*alpha*i
(C13+C44)*alpha*i,C44+C33*alpha^2-pho*v^2,e15+e33*alpha^2
(e15+e31)*alpha*i,e15+e33*alpha^2,-(es11+es33*alpha^2)];
f = det(A);
==============================================================================
所有係數原汁原味複製貼上
最後解答是 anser . 有對應的 v 跟 alpha
你先看看吧
--
數值方法 : 求根值 . ODE系統 . DDE系統. DAE系統. PDE系統
非線性聯立系統. 化工程序控制系統 . 力學系統
RLC 系統 . MCK 系統 . Lalpace 轉換 系統
其他 : 工程數學 .微積分.化工程序.輸送現象.化工程序線性控制
Simulator . MATLAB tool box . Simulink . FORTRAN
不會影像處理.....囧>
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.118.160.160
→ hilosi:1個v可能對應到多個alpha喔... 03/06 11:42
→ zupo:嗯...的確很麻煩.不過這就不知道該怎麼解了. 03/06 20:19
→ zupo:但是我給的應該有滿足 :P 03/06 20:20
推 alamabarry:polyeig...科科....壓電晶體波傳是吧 03/06 22:14
→ zupo:那是?...^^ 03/06 22:44
推 fishooooo:謝謝～～另外四樓＠＠ 03/07 16:08