 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 4316 | 回復(fù): 10 | ||
[求助]
求助matlab for循環(huán)中的賦值問(wèn)題。 已有2人參與
|
||
|
dm=0.001*2*pi*c; w=-0.5*2*pi*c:dm:0.5*2*pi*c; Kx=zeros(size(w)); Ay=zeros(size(w)); for m=1:length(w); syms ay kx ky; y=(0.0382*8*(pi^3)*(c^4)*sqrt(ainf-(wp.^2)./(w.^2)))./((wp^2).*w); axx=ainf*((1-(wp*(w+1i*y)*wp))./(w.*(((w+1i*y).^2)-((wc)^2)))); %無(wú)單位 axy=ainf*(((wp^2)*wc*1i)./(w.*(((w+y*1i).^2)-((wc)^2)))); %無(wú)單位 azz=ainf.*((1-(wp^2))./(w.*(w+1i*y))); %無(wú)單位 a1=real(axx); a2=imag(axx); a3=real(axy); a4=imag(axy); a5=real(azz); a6=imag(azz); f1=ay-sqrt((kx.^2)-((w./c).^2)); f2=kx-((sqrt(((a4.*ay./2).^2))+(a1.*(w./c).^2)-((a1.*ay.*ky)./(tan(ky.*t))))-(a4.*ay./2)); f3=ky-(sqrt((((w./c).^2).*(a1-((a4.^2)./a1)))-kx.^2)); [kx,ay]=solve('f1','f2','f3'); Kx(m)=kx; Ay(m)=Ay; end 運(yùn)行提示：在賦值 A(  = B 中，A 和 B 中的元素?cái)?shù)目必須相同。出錯(cuò)行數(shù)： Kx(m)=kx; Ay(m)=Ay; 跪求大神解惑。。。  |
回復(fù)此樓» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有5人回復(fù)
-
一志愿華中科技大學(xué)，080502，354分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有6人回復(fù)
-
085410人工智能專碩317求調(diào)劑（0854都可以）
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助MATLAB編程
已經(jīng)有8人回復(fù)
- 求助matlab寫變限積分
已經(jīng)有4人回復(fù)
- 求助Matlab
已經(jīng)有0人回復(fù)
- 求助matlab圖像
已經(jīng)有0人回復(fù)
- 求助matlab求解三元超越方程
已經(jīng)有4人回復(fù)
- 求助matlab三維圖糾正程序
已經(jīng)有2人回復(fù)
- 求助matlab有限差分程序
已經(jīng)有1人回復(fù)
- 求助matlab寫循環(huán)
已經(jīng)有6人回復(fù)
- 求助MATLAB畫圖，二維曲線，總是出錯(cuò)，求助高手
已經(jīng)有3人回復(fù)
- matlab求助
已經(jīng)有4人回復(fù)
- Matlab求助
已經(jīng)有6人回復(fù)
- 求助matlab問(wèn)題
已經(jīng)有5人回復(fù)
- 求助matlab程序
已經(jīng)有3人回復(fù)
- 求助matlab問(wèn)題
已經(jīng)有5人回復(fù)
- 高金求助matlab解微分方程組
已經(jīng)有12人回復(fù)
- 求助matlab程序
已經(jīng)有4人回復(fù)
- 求助matlab積分的問(wèn)題
已經(jīng)有4人回復(fù)
- matlab求助
已經(jīng)有1人回復(fù)
2樓2016-09-16 13:07:16
guokeqin
金蟲 (正式寫手)
- 應(yīng)助: 14 (小學(xué)生)
- 金幣: 5441.2
- 紅花: 6
- 帖子: 823
- 在線: 118.4小時(shí)
- 蟲號(hào): 1938806
- 注冊(cè): 2012-08-14
- 性別: GG
- 專業(yè): 光學(xué)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
solve函數(shù)，不能計(jì)算點(diǎn)乘，建議你采用solve函數(shù)計(jì)算出來(lái)函數(shù)關(guān)系式，然后手動(dòng)輸入方程再操作 發(fā)自小木蟲Android客戶端 |
6樓2016-09-17 06:57:58
3樓2016-09-16 14:43:59
4樓2016-09-16 15:38:26
FMStation
至尊木蟲 (知名作家)
- 應(yīng)助: 591 (博士)
- 貴賓: 0.03
- 金幣: 18333.2
- 紅花: 97
- 帖子: 8891
- 在線: 1489.8小時(shí)
- 蟲號(hào): 2400059
- 注冊(cè): 2013-04-04
- 專業(yè): 計(jì)算機(jī)應(yīng)用技術(shù)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
https://www.mathworks.com/help/symbolic/solve.html Return Numeric Solutions Try solving the following equation. The symbolic solver cannot find an exact symbolic solution for this equation, and therefore issues a warning before calling the numeric solver. Because the equation is not polynomial, an attempt to find all possible solutions can take a long time. The numeric solver does not try to find all numeric solutions for this equation. Instead, it returns only the first solution it finds. |
5樓2016-09-16 21:56:31
|
我把點(diǎn)乘去掉顯示： 在賦值 A( = B 中，A 和 B 中的元素?cái)?shù)目必須相同。我重新編了一下solve的公式： kx=solve(tan(t*(((sqrt(((w(m)/c)^2)*(a1-((a4^2)/a1)))))-(kx^2)))==(a1*sqrt(((w(m)/c)^2)*(a1-((a4^2)/a1))-(kx^2))*sqrt((kx^2)-((w(m)/c)^2)))/((a1*((w(m)/c)^2))-(kx^2)-(a4*kx*sqrt((kx^2)-((w(m)/c)^2)))),kx); 但是在kx賦值到Kx的時(shí)候出現(xiàn)錯(cuò)誤，另外，solve語(yǔ)句單獨(dú)拿出來(lái)進(jìn)行計(jì)算的時(shí)候會(huì)出現(xiàn)警告： 警告: Cannot solve symbolically. Returning a numeric approximation instead. 這個(gè)時(shí)候我能輸出kx。很疑惑。 |
7樓2016-09-17 22:28:39
|
大神，這個(gè)我了解了，所以嘗試把solve語(yǔ)句修改了一下，單獨(dú)計(jì)算kx=solve(tan(t*(((sqrt(((w/c)^2)*(a1-((a4^2)/a1)))))-(kx^2)))==(a1*sqrt(((w/c)^2)*(a1-((a4^2)/a1))-(kx^2))*sqrt((kx^2)-((w/c)^2)))/((a1*((w/c)^2))-(kx^2)-(a4*kx*sqrt((kx^2)-((w/c)^2)))),kx);時(shí)，現(xiàn)出現(xiàn)警告： 警告: Cannot solve symbolically. Returning a numeric approximation instead. 這個(gè)時(shí)候能輸出kx的數(shù)值解，所以我就比較疑惑為什么有數(shù)值解還不能賦值成功。。。 |
8樓2016-09-17 22:31:36
|
我把方程組進(jìn)行處理，整理成一個(gè)公式進(jìn)行計(jì)算，代碼如下： clear; clc; a0=8.85*1e-12; %c2/m2N ainf=12.37; %極限高頻相對(duì)介電常數(shù) 無(wú)單位 %y=pi*1e11; %衰減頻率 hz B=3; %磁場(chǎng)強(qiáng)度 T或kg/As2 %T=185; %溫度 K t=1.361*1e-6; %厚度 m e=1.6*1e-19; %單位電荷量 C me=9.11*1e-31; %電子質(zhì)量 kg m=0.033*me; %載流子質(zhì)量 kg N=8*1e23; c=3e8; %N=(5.76*1e20)*(T^1.5)*exp(-(0.13/((8.625*1e-5)*T))); %載流子濃度 m-3 此處的0.0151為KB*T KB為玻爾茲曼常數(shù)8.625e-5 eV/K,所以0.0151單位為eV wc=e*B/m; %回旋頻率 rad/s wp=sqrt(N*(e^2)/(a0*m)); %等離子頻率 rad/s %w=0:pi*2*(2*(1e12))/1000:pi*2*2*(1e12); %定義頻率范圍 rad/s %w=1.21*2*pi*1e12; dm=0.0001*2*pi*c; w=0.4*2*pi*c:dm:0.5*2*pi*c; Kx=zeros(size(w)); Ay=zeros(size(w)); for m=1:length(w); syms kx; y=(0.0382*8*(pi^3)*(c^4)*sqrt(ainf-(wp.^2)./(w(m).^2)))./((wp^2).*w(m)); axx=ainf*((1-(wp*(w(m)+1i*y)*wp))./(w(m).*(((w(m)+1i*y).^2)-((wc)^2)))); %無(wú)單位 axy=ainf*(((wp^2)*wc*1i)./(w(m).*(((w(m)+y*1i).^2)-((wc)^2)))); %無(wú)單位 azz=ainf.*((1-(wp^2))./(w(m).*(w(m)+1i*y))); %無(wú)單位 a1=real(axx); a2=imag(axx); a3=real(axy); a4=imag(axy); a5=real(azz); a6=imag(azz); kx=solve(tan(t*(((sqrt(((w(m)/c)^2)*(a1-((a4^2)/a1)))))-(kx^2)))==(a1*sqrt(((w(m)/c)^2)*(a1-((a4^2)/a1))-(kx^2))*sqrt((kx^2)-((w(m)/c)^2)))/((a1*((w(m)/c)^2))-(kx^2)-(a4*kx*sqrt((kx^2)-((w(m)/c)^2)))),kx); A=char(kx); Kx(m)=A; end 顯示： 在賦值 A( = B 中，A 和 B 中的元素?cái)?shù)目必須相同。出錯(cuò) yanzheng (line 38) Kx(m)=A; 。。。。。 |
9樓2016-09-17 22:34:08
czcdxmc
木蟲 (正式寫手)
MEMS小學(xué)生
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 5212.9
- 散金: 70
- 紅花: 8
- 帖子: 845
- 在線: 101.7小時(shí)
- 蟲號(hào): 2088420
- 注冊(cè): 2012-10-26
- 性別: GG
- 專業(yè): 半導(dǎo)體微納機(jī)電器件與系統(tǒng)
10樓2016-09-18 08:17:15
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 duguwuhao 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +4 | 蔡大美女 2026-03-13 | 4/200 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +5 | 陌の森林 2026-03-18 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +9 | 秋有木北 2026-03-14 | 9/450 |
|
|
|
[考研] 266求調(diào)劑 +5 | 陽(yáng)陽(yáng)哇塞 2026-03-14 | 9/450 |
|
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
|
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|