| 7 | 1/1 | 返回列表 |
| 查看: 2202 | 回復(fù): 6 | ||
| 【懸賞金幣】回答本帖問題,作者小蟲飛更高將贈(zèng)送您 5 個(gè)金幣 | ||
小蟲飛更高新蟲 (初入文壇)
|
[求助]
matlab光強(qiáng)編程,提示矩陣維度 必須一致,不知是哪里出現(xiàn)的問題,代碼在下面 已有1人參與
|
|
|
代碼如下: clear all; k=10; f=1004.8; x=-1.5:0.01:1.5; y=-1.5:0.01:1.5; [x1,y1]=meshgrid(x,y); r=sqrt(x1.^2+y1.^2); Phi=atan(y1./x1); I=(i*k*exp(i*(pi/2+k*f)).*sin(2.*Phi).*quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r.*sin(theta)),2.23,pi)).^2; subplot(2,1,1); mesh(x1,y1,I); axis([-1.5 1.5 -1.5 1.5 0 1]); subplot(2,1,2); subimage(I*255); axis off; 出現(xiàn)錯(cuò)誤如下: 矩陣維度必須一致。 出錯(cuò) formula>@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r.*sin(theta)) 出錯(cuò) quadl (line 62) y = feval(f,x,varargin{:}); y = y( .';出錯(cuò) formula (line 9) I=(i*k*exp(i*(pi/2+k*f)).*sin(2.*Phi).*quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r.*sin(theta)),2.23,pi)).^2; |
至尊木蟲 (著名寫手)

新蟲 (初入文壇)
送紅花一朵 |
你好,我按照你的編程輸入了下,還是有提示錯(cuò)誤,請(qǐng)你在幫我解答下, clear all; k=10; f=1004.8; x=-0.5:0.01:0.5; y=-0.5:0.01;0.5 [x1,y1]=meshgrid(x,y); r=sqrt(x1.^2+y1.^2); [K,L]=size(r); for m=1:k for n=1:L Phi=atan(y1(n)./x1(m)); r0=r(m,n); eq=quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r0.*sin(theta)),2.23,pi); I0=(i*k*exp(i*(pi/2+k*f))*sin(2*Phi)*eq).^2; I(m,n)=I0; end end Re=real(I);Im=imag(I); s.FaceColor='flat'; %subplot(2,1,1); C=rand(K,L); mesh(x1,y1,real(I).C);hold on 提示錯(cuò)誤如下: formula ans = 0.5000 未定義與 'double' 類型的輸入?yún)?shù)相對(duì)應(yīng)的函數(shù) 'Besselk'。 出錯(cuò) formula>@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r0.*sin(theta)) 出錯(cuò) quadl (line 62) y = feval(f,x,varargin{:}); y = y( .';出錯(cuò) formula (line 13) eq=quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r0.*sin(theta)),2.23,pi); >> |
新蟲 (初入文壇)
|
剛我又調(diào)試了下,出現(xiàn)如下錯(cuò)誤,您再幫我看下: k=10; f=1004.8; x=-0.5:0.01:0.5; y=-0.5:0.01:0.5; [x1,y1]=meshgrid(x,y); r=sqrt(x1.^2+y1.^2); [K,L]=size(r); for m=1:k for n=1:L Phi=atan(y1(n)./x1(m)); r0=r(m,n); eq=quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*besselk(2,k.*r0.*sin(theta)),2.23,pi); I0=(i*k*exp(i*(pi/2+k*f))*sin(2*Phi)*eq)^2; I(m,n)=I0; end end Re=real(I);Im=imag(I); s.FaceColor='flat'; %subplot(2,1,1); C=rand(K,L); mesh(x1,y1,real(I).C);hold on axis([-1.5 1.5 -1.5 1.5 0 1]); subplot(2,1,2); subimage(I*255); axis off 提示錯(cuò)誤: 定義變量 "real" 或類 "real"。 出錯(cuò) para (line 22) mesh(x1,y1,real(I).C);hold on |
至尊木蟲 (著名寫手)

新蟲 (初入文壇)
送紅花一朵 |
我試了試其他筆記本的電腦上的版本matlab,還是提示未定義變量‘real ’或類‘real’,請(qǐng)您運(yùn)行下我的程序,看看是什么問題呢? |
新蟲 (初入文壇)
|
定義變量的問題解決了,但是又遇到了一個(gè)索引函數(shù)編制的問題,問題在代碼后面: clear all; close all; clc; k=10; f=1004.8; x=-0.5:0.01:0.5; y=-0.5:0.01:0.5; [x1,y1]=meshgrid(x,y);%直角坐標(biāo)轉(zhuǎn)化為極坐標(biāo) r=sqrt(x1.^2+y1.^2); [K,L]=size(r); for m=1:k for n=1:L Phi=atan(y1(n)./x1(m)); r0=r(m,n); eq=quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*besselk(2,k.*r0.*sin(theta)),2.23,pi); I0=(i*k*exp(i*(pi/2+k*f))*sin(2.*Phi)*eq).^2; I(m,n)=I0; end end Re=real(I); Im=imag(I); s.FaceColor='flat'; %subplot(2,1,1); C=rand(K,L); mesh(x1,y1,real(I).C); hold on ; axis([-1.5 1.5 -1.5 1.5 0 1]); subplot(2,1,2); subimage(I*255); axis off 錯(cuò)誤: 不能使用 {} 或 . 索引為函數(shù)編制索引。 出錯(cuò) para (line 25) mesh(x1,y1,real(I).C); |
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703求調(diào)劑 +7 | c297914 2026-03-06 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)工程求調(diào)劑 +12 | 化工人999 2026-03-04 | 12/600 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)0856材料與化工,前三科206,總分283,求調(diào)劑 +7 | 26考研求調(diào)劑 2026-03-06 | 7/350 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|
[考研] 287求調(diào)劑 +3 | 看看我. 2026-03-05 | 6/300 |
|
|
[考研] 085600,一志愿鄭州大學(xué),280分求調(diào)劑 +7 | Wuqi725 2026-03-05 | 7/350 |
|
|
[考研] 282求調(diào)劑 +7 | 夕~日 2026-03-05 | 8/400 |
|
|
[考研] 274求調(diào)劑 +9 | 一個(gè)學(xué)習(xí)者 2026-03-04 | 9/450 |
|
|
[考研] 沒上岸的看過來 +3 | tangxiaotian 2026-03-01 | 5/250 |
|
|
[考研] 274環(huán)境工程求調(diào)劑 +6 | 扶柳盈江 2026-03-05 | 6/300 |
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
[考研] 266求調(diào)劑 +7 | 哇塞王帥 2026-03-03 | 7/350 |
|
|
[考研] 能動(dòng)297求調(diào)劑,本科川大 +4 | 邵11 2026-03-04 | 4/200 |
|
|
[考研] 一志愿西交大材料學(xué)碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-02 | 9/450 |
|
|
[考研] 0854復(fù)試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
[考研] 【2026 碩士調(diào)劑】課題組 招收調(diào)劑生 +3 | 考研版棒棒 2026-03-02 | 5/250 |
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|